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GENERAL INTRODUCTION AND SUMMARY 

This final report discusses the work accomplished over the past lour years on 
segregation effects. The past year, this program was partially supported under NASA 
contract NAS 8-27891, "Segregation Effects During Solidification in Weightless Melts. " 
The contract covered the period from January 1, 1972 to December 31, 1972. 

Two types of melt segregation effects were studied: evaporative segregation, or 
segregation due to surface evaporation; and freezing segregation, or segregation due to 
liqu’d-solid phase transformation. 

These segregation effects are close., related. In fact, evaporative segregation 
always precedes freezing segregation to some degree and must often be studied prior to 
performing meaningful solidification experiments. This is particularly true since 
evaporation may cause the melt composition, at least at the critical surface regions or 
layers to be affected manyfold within seconds so that the surface region or layer melting 
point and ether thermophysical properties, nucleation characteristics, base for under- 
cooling, and critical velocity to avoid constitutional supercooling, may be completely 
unexpected. 

An important objective was, therefore, to develop the necessary "normal evapor- 
ation equations" for predicting the compositional changes within specified times at tem- 
perature and to correlate these equations with actuul experimental data collected from 
the literature. An evaporative congruent temperature (or equi-evaporative temperature) 
was then defined and listed for various binary or ternary alloys. The application of this 
unique temperature in explaining, predicting, or planning "anomalous" evaporative or 
constitutional melting (on cooling) or solidification (on heating) experiments were dis- 
cussed. We then computed the reactive jetting forces due to surface evaporation and, 
in particular, showed that these forces can be very substantial on a differentially heated 
sample and may completely destroy the unique zero-gravity environment in space manu- 
facturing. In addition, these jetting forces may initiate surface deformation and vibration 
or other fluid disturbances, and may even produce new types of important convection 
currents probably never anticipated. These studies also showed which sample materials 
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are preferable and which to avoid, and what impurities are effective as stabilizing 
influences. The relation between normal evaporation and normal freezing was then 
considered. Finally, applications of evaporation to space manufacturing concerning 
material loss and dimensional control, compositional changes, evaporative purification, 
surface cooling, materials standards, and freezing data interpretation were briefly 
described. 

In the area of segregation due to solidification, we started with an explanation of 
the normal freezing process and its successful use in the semiconductor industry. 

Various constitutional diagrams demonstrated the desirability of using non-constant 
segregation coefficient techniques in metallurgical studies. We then stated the basic normal 
freezing differential equation, together with its solutions for cases w fc ere the liquidus and 
solidus are quadratic, cubic, high-degree polynomial, and exponential functions of the melt 
temperature. The meaning of constant segregation coefficient was discussed, together 
with the associated errors due to curvatures of the liquidus and solidus lines and the 
best value of constant segregation coefficient for a given solidification experiment. 
Numerical methods for computing the normal freezing behavior were then given. Finally 
and as an example, the steady state solidification of the Ni-Sn system under conditions of 
limited liquid diffusion and non-constant segregation coefficients was described. 

The report ends with some suggestions for future work. 
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PART I 

EVAPORATIVE SEGREGATION 
Section 1. INTRODUCTION 


In normal evaporation, we assume that the evaporating alloy is always homogeneous 
in composition, that the vapor is instantly removed, and that the alloy follows Raoult's 
law (Dushman 1962). Such conditions exist or are approached in an induction-stirred, 
melt-in-vacuum, or liquid drop-in-space. 

We deal here mainly with the normal evaporation of binary alloys. i< particular, 
we study the evaporative segregation patterns, i.e. , the type and degree of enrichment 
or depletion of solute in the evaporating source at different evaporation temperatures and 
times. 

Normal evaporation is important in space melting and solidification for the following 
reasons: 

• Significant evaporation of alloy components always occurs at high temperatures 
in space vacuum environments 

• High -temperature evaporation of alloys is generally a neglected area of 
systematic research. Yet, unless the complete evaporative segregation 
behavior is understood and analyzed, solidification and its related segregation 
effects may not be properly studied because of ill-defined initial conditions. 

Before the liquid alloy can be controllably solidified or even melted, there is 
invariably some surface evaporation to cause changes in composition, freezing 
temperature, supercooling characteristics, nucleation and growth morphology 
conditions, and the like 

• Controlled space evaporation probably most closely meets the requirements 
of normal evaporation. We may thus be able to obtain material purity or 
evaporation standards, thermal properties, or even basic thermodynamic 
properties such as heat of evaporation, activity coefficients, sticking 


1-1 





1 

« 


t* 

* 


% 


I 

5 


coefficients, etc. , that are difficult or impossible to obtain on earth 

• Normal evaporation is a much simpler process than normal freezing, 
since the former does not involve such complicated phenomena as 
nucleation, phase transformation, and constitutional or nonconstitutional 
supercooling. Thus in normal evaporation we may, for specific geo- 
metries or alloy systems, ideally isolate and investigate such other 
phenomena as heat conduction or radiation, liquid or solid diffusion, 
fluid dynamics, and convection currents. Exact knowledge of these 
phenomena is necessary to understand solidification 

• Evaporation causes surface cooling due to the heat of evaporation. This 
evaporative cooling effect is particularly important in low-melting materials 

• Different rates of evaporation at various surface regions give rise to 
unbalanced forces and momenta that may produce erratic or unwanted 
accelerations, surface distortions and vib-ations, exceedingly large 
"equivalent gravities", and possibly new types of powerful convection 
currents in zero-gravity conditions, as will be discussed later 

• Evaporation may cause the surface composition of certain unwanted or 
unsuspected impurities to be increased a thousandfold or millionfold 
within seconds so that the layer's melting point and other thermophysical 
properties, nucleation characteristics, base for undercooling, and critical 
velocity to avoid constitutional supercooling, may be completely unexpected. 
In fact, anomalous "constitutional" or evaporative melting on cooling, or 
solidification on heating, is possible because of surface evaporation 
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Section 2. DIFFERENTIAL EQUATION OF NORMAL EVAPORATION 
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Zinsmeister (1964) and Li (1966) gave the exact equation for the normal evaporation 
of binary alloys. This equation, relating the concentration m of the solute in the 
evaporating alloy at time t, when the mole fraction of the alloy remaining is F, is 



where N and N q are, respectively, the number of moles of both solvent and solute at 
evaporating times t = t and t = 0; m Q is the initial molar concentration of the evaporating 
alloy; and U and V are, respectively, the evaporation rate of pure solute and solvent. 

2 

The evaporation rates in mol/m /sec for pure solute and solvent are 

A -B /T i 

U = K 10 U U (M u T)" 2 (2) 


and 

A -B /T j 

V = K 10 V V (M T)’ 2 (3) 

-9 

where K = 5.833 x 10 <y, a = 1 for most metals (Dushman 1962); M y and are, 

respectively, the molecular weights of the solute and solvent atoms; T is the evaporating 
temperature in degrees Kelvin; and A^ and B^ or A v and B v are evaporating constants for 
the solute or solvent. Table 1-1 gives the values of the constants A and B for 51 different 
elements, as given by Dushman (1962). Data on the melting point and density of these 
elements were obtained from Smithell (1967). 
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Element 

A 

B/1000 

At. Wt 

Melting Pt 

Density 





K 


Ag 


14.27 

107.9 

1234 

10.5 

A1 

jj^K' rBKTt^Hj 

15.94 

26.98 

933.3 

2.7 

Au 

nrs 

gpHlff 1 W*r^Mf 

17.58 

197.0 

1336 

19.3 

B 


29.62 

10.82 

2300 

2.34 

Ba 


8.76 

137.4 

977.1 

3.5 

Be 

12.01 

16.47 

9.013 

1556 

1.848 

Bi 

11.18 

9.53 

209.0 

544.1 

9.8 

C 

15.73 

40.03 

12. 01 

3973 

2.25 

Ca 

11.22 

8.94 

40.08 

1123 

1.54 

Cb 

14.37 

40.40 

92.91 

2688 

8.57 

Cd 

11.56 

5.72 

112.4 

594.1 

8.64 

Ce 

13.74 

20,10 

140.1 

1077 

6. 75 

Co 

12.70 

21.11 

58.94 

1765 

8.9 

Cr 

12.94 

20.00 

52. 01 

2176 

7.1 

Cs 

9.91 

3.80 

132.9 

301.8 

1.87 

Cu 

11.96 

16.98 

63.54 

1356 

8.96 

Fe 

12.44 

19.97 

55.85 

1812 

7.87 

Ga 

11.41 

13.84 

69.72 

302.8 

5.91 

Ge 

11.71 

18.03 

72.60 

1210 

5.32 

In 

11.23 

12.48 

114-8 

429. 7 

r> o 
t . O 

Ir 

13.07 

31.23 

192.2 

2716 

22.4 

!• 

10.28 

4.48 

39.10 

336.3 

.86 

La 

11.60 

20.85 

138.9 

1103 

6.19 

Li 

10.99 

8.07 

6.94 

459.7 

.534 

Mg 

11.64 

7.65 

24.32 

924 

1.74 

Mn 

12.14 

13.74 

54.94 

1517 

7.4 

Mo 

11.64 

30.85 

95.95 

2883 

10.2 

Na 

10.72 

5.49 

22.99 

371 

.97 

j Ni 

12.75 

20.96 

58.71 

1726 

8.9 

Os 

13.59 

37. 00 

190.2 

2970 

22.5 

t»b 

10.77 

9.71 

207.2 

600.6 

11.68 

Pd 

11.78 

19.71 

106.7 

1825 

12. 

Pt 

12.53 

27.28 

195.1 

2045 

21.45 

Rb 

10.11 

4.08 

85.48 

312.1 

1.53 

Rh 

12.94 

27.72 

102.9 

2239 

12.4 

Ru 

13.50 

33.80 

101.1 

2700 

12.2 

Sb 2 

11.15 

8.63 

243.5 

903. 7 

6.68 

Sc 

11.94 

18.57 

44.96 

1673 

2.99 

Si 

l 

12.72 

21.30 

28.09 

1685 

2.34 










TABLE 1-1 EVAPORATING CONSTANTS FOR 51 CHEMICAL ELEMENTS (Continued) 


Element 

A 

B/1000 

At. Wt 

Melting Pt 
K 

Density J 

Sn 

10.88 

14.87 

118.7 

505 

7.3 

Sr 

10.71 

7.83 

87.63 

1033 

2.6 

Ta 

13.04 

40.21 

180.9 

3270 

16.6 

Th 

12.52 

28.44 

232.0 

1968 

11.85 

Ti 

12.50 

23.23 

47.90 

1941 

4.5 

T1 

11.07 

8.96 

204.4 

576. 8 

11.85 

U 

11.59 

23.31 

238.1 

1406 

19.97 

V 

13.07 

25.72 

50.95 

2290 

6.1 

W 

12.40 

40.68 

183.9 

3653 

19.3 

Y 

12.43 

21.97 

88.91 

1773 

4.47 

Zn 

11.63 

6. 56 

65.38 

692.6 

7.14 j 

Zr 

12.33 

30.26 

91.22 

2130 

6.49 


1 : 

1: 

i; 


Differentiating Eq (1) yields 


but 


dm 


[mU+(l-m)V]N d 


dN = - [mU+(l-m)V] A dt 


where A is the evaporating area of the alloy, assumed constant here. 
Substituting Eq (5) into Eq (4) results in 

dt = Gm y (l-m)°dm 

where N a-m 

C = - — 2 — — 

rv-fl 

A(U-V)m" 

o 

a = (2V-U)/U-V) 

P = (V-2U)/(U-V) 


(4) 


(5) 


( 6 ) 

(7) 

( 8 ) 

(9) 
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2.1 EVAPORATING RATES AND EVAPORATIVE SEGREGATION COEFFICIENT 

The evaporating rates of elements, being exponential functions of the evaporating 

temperature as shown in Eq (2) and (3), are strongly temperature-dependent. At the 

-15 -14 

melting point of nickel, i.e., 1726K, these rates vary from 1. 30 x 10 , 1.05x10 

-13 -12 -13 2 

1.25 x 10 , 1.68 x 10 , and 2. 76x10 Kg/m /sec respectively, for W, Ta, Cb, 

C, and Os to 26,200, 8270, 7670, 7240, 4250, 2330, and 1120 Kg/m 2/ sec for Cd, Cu, 

-15 

Zn, Rb, K, Na, and Mg. The molar evaporating rates vary from 7.11 x 10 , 

-14 -12 -11 -14 2 

5.57 x 10 , 1.35 x 10 , 1.45 x 10 , and 1.40x10 moles/m /sec respectively 

for W, Ta, Cb, Os, and C to 2.34 x 10 5 , 1.31 x 10 5 , 1.17 x 10°, 1.09 x 10 & , 1.00 x 10 5 , 

4 2 

8.5x10 moles/m /sec for Cd, Cu, Zn, K, Na, and Rb. The evaporating rates can also 

be expressed in terms of the velocity of interfacial movement in m/sec, i.e. , from 

-16 -15 -13 -12 -12 

6. 74 x 10 , 6.31 x 10 , 1.36 x 10 , 1.23 x 10 , and 7.45 x 10 respectively 

for W, Ta, Cb, Os, and C to 4.96 x 10 6 , 4. 72 x 10 6 , 3. 03 x 10 6 , 2.40 x 10 (> , 1.07 x 10 6 , 

5 

9.25 x 10 for K, Rb, Cd, Na, Zn, and Cu. For the solvent nickel in nickel alloys, the 

-42 -32 

evaporating rates at 1726K are 4.35 x 10 Kg/m /sec, 7.42 x 10 moles/m /sec, and 
-8 

4.89 x 10 m/sec (or 489 A /sec). Note that the evaporating rates of the various solute 
elements in nickel generally vary by many orders of magnitude. 

When the solute concentration in the evaporating alloy is m mole fraction, the 
solute concentration of the. vapor is mU/[mU+(l-m)V] , where V and U are, as defined 
previously, the molar evaporating rates of the solvent and solute, respectively. The 
evaporative segregation coefficient k^ is, therefore, 

k = U/[mU+(l-m)v] 

which varies with m and U and V and, hence, the evaporating temperature T. 

Since U and V are generally widely different, we have the following two usual cases: 

• Us>V: k v = l/m, i.e., the evaporative segregation coefficient depends only 
on m, but not on T. 

• For the more important case of U v>V, in which the solute may be concentrated 
manyfold on the surface, 
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k = — — 

v (l-m)V 


If, in addition to U « V, m is small, i.e. , m «1, k^ = U/V = constant 
at a given T. 

Thus, for small, particularly trace, amount of low -evaporative impurities or 

solutes in nickel alloys evaporating at 1724K, the evaporating segregation coefficient 

-12 -11 -10 -9 

k = U/V can be as small as 1. 38 x 10 , 1. 29 x 10 , 2. 78 x 10 , 2. 51 x 10 , and 

^ -8 

1.53 x 10 for W, Ta, Cb, Os, and C, respectively. 

For highly evaporative impurities (relative to nickel), the evaporative segregation 
coefficient at 1726K is 

k = 1/m for Us> V 
v 

-6 -9 

which can be very large for small values of m. Thus, if m - 10 or 10 , the instan- 

6 9 

taneous evaporative segregation coefficient is 10 or 10 . Since U »V, m is rapidly 
decreasing. Hence, the evaporative segregation coefficient may rapidly decrease by 
many orders of magnitude within a single experimental run, even for the same impurity, 
when U »V. 

2.2 SPECIAL CASE SOLUTIONS OF THE NORMAL EVAPORATION EQUATION 

Equations (6) - (9) enable us to determine the evaporating time, t, for an evapor- t 

ating alloy to reach a specific solute concentration, m. Unfortunately, these equations 
are not exactly solvable in the general case. All of the following special and important 
cases, (except Case V), however, are solvable in closed forms. 

Case I: The solute is much more evaporative than the solvent, i.e. , U »V; 
or a -1 and 8= -2. In this case 

G dm 

dt, = ^ (10) 

A m(l-m) 4 
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where 


°1=- 


N (1-m ) 
o o 

AU 


and the evaporating time is: 


r a-„ 

*1 = G 1 ['“ rTff: 


)m 

nT 


m-m 


(l-m Q ) (1-m) 


(ID 


In terms of the initial and final weight fractions of the solute in the evaporating 
alloy, w^ and w respectively, the evaporating time can be shown to be: 


r 

h G i [ ; n w Q (i- 


)W 

W) 


w-w 


o 


M 


U 


(1-W) (1-W Q ) • M u 


] 


(11a) 


For given m Q and evaporating temperature T or evaporating rate U, the evaporation 
constant 0, Is proportions! to N/A. For a dilute evaporating alloy in spherical form of 
diameter d, density p, and weight W, is proportional to dp/W. Hence, for dilute 
alloys and given initial and final conditions, the evaporatirg time is directly proportional 
to the diameter d of the evaporating source. 


When m ~ m Q , the second te: m in the bracket in Eq. 11 (or 11a) may be nearly 0. 
The evaporation time for this first case, t^, is then a logarithmic function of m (or w), 
as has been experimentally observed (see e.g. , Tefelske and Chang 1972). 


Also, when m m Q 0, 



N 

o 

AU 



( 12 ) 


When, however, m ~ m Q — ht the logarithmic term in Eq 11 or 11a may be small. 
The evaporating time tj is then approximately linearly related to the final solute con- 
centration (m or *v) in the source. 

Case II: The solvent is much more evaporative than the solute, i.e. , V » U; or 
a= -2 and 0 = -1. In this case 
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P 





dl 2 = 


1 

(l-m)m 


(13) 


where 


m N 
r ° 9 
G 2 AV 


and 


l 2 “ G 2 


(1-m )m m-m 

o o 

f n — + 

m (1-m) m m 

o o 


(14) 


In terms of weight fractions 


‘ 2 ' G 2 


(1-W )W w-w „ 

^ o o_ Mu_ 

,n W (1-W) + W W Mv 
o o 


(14a) 


For the initial stages of evaporation, i.e. , m m , t also becomes a logarithmic 

O c, 

function of m, as has been observed. On the other hand, when m m 0, t may be 

O Ld 

approximately linearly related to 1/m. 

Case 111: a = (2V-U)/(U-V) = 0, i.e. , U = 2V and fl = -3 




(1-mf 2 - (l-m 0 )‘ 2 


1 


where 


N (1-m ) 

G _ _ _2 2 _ 

3 2AVm 


(15) 


Case IV: 0= (V-2U)/(U-V) = 0, i.e., V = 2U and a = -3. Jn this case 


t. = G 

4 4 




(16) 


where 


N m 

r _ - o o 

4 2AU(l-m ) 


a. 

IT 

ii 
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Equations (15) and (16) show that under some conditions, the evaporating time, 

-2 -2 

t, is better represented by linear functions of (1-m) or m than by logarithmic 
functions of m. 


Case V: For relatively dilute alloys, i.e. , m and Bm «1, the following solution 
by series expansion can be obtained from Eq (6): 

.2 


dt,. = Gm ll - pm + 


* 

1 


(Bm) 


+ . . . + (- 1 ) 


i (ftn) 1 
il 


dm 


(17) 


and 


‘5=° 


0+1 0+1 
m - m 
o 

o + l 


o+2 o+2 a+i + l » +1+1 

m - m m - m 

° &+•••+ (-1) - -- -• 7 T -", P ••• (18) 


o 


+ 2 


(o + i + 1 ) i 1 


th 


For computer calculations, it is desirable to know the ratio of the i term to the 


- ,th 

(l-l) term. Thus, 


T. B(o+l)(m a+1+1 - n" +1+1 ) 

( a+i o+i \ 

m -m o j 


(19) 


i-1 i (o + i+1) 


which is generally less than 0m/i or BnWi. 


Because Bm « 1 and i constantly increases with each additional term, this series 
converges rapidly unless 0 is very large, i.e. , unless U = V, which leads to Case VI. 

Case VI: The solute and solvent are evaporating at equal rates, i.e. , U = V. In 
this case, we would expect 

m = m for all t . (20) 

o 


There is, then, no evaporative segregation, i.e., there is neither solute enrichment nor 
depletion in the evaporating source. 

For any pair of solvent and solute, there is a unique temperature, T , at which 

s 

U = V and, hence, the alloy concentration remains stable or constant. From Eq (2) and 
(3), 
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B - B 
u v 


( 21 ) 


T = 


s A - A - 0. 5 log(M /M ) 
u v u v 


This unique "equi -evaporative, or evaporative congruent, temperature' 1 T has 

s 

some important implications and applications, as will be shown. 


Case VII: With extremely dilute alloys, i.e. , m m Q 0, ve have (Li, 1.96") 


F = N/N = (m/m ) ! 
o o 


a+1 


( 22 ) 


Hence, 


t 7*X?[ 1 - (m/n \/ +1 ] 


(23) 


2.3 EXAMPLES OF COMPUTATION 


A generalized computer program has been written to analyze the evaporative seg- 
regation pattern (i.e. , type and degree), based on the equations given above, on various 
potential aerospace materials suitable for space experiments. Computations have been 
made that give the maximum allowable time (often less than a minute or second) for space 
experiments without material loss or compositional changes of a specified amount for 
various solute elements in nickel, iron, or other alloys. 

2. 3. 1 Nickel Alloys 

Several nickel alloys have been selected for the M512 Skylab experiments. As an 
example of the use of the various derived equations, the evaporation behavior of an alloy 
containing 8% by weight of A1 in Ni at the melting point of pure Ni (i.e. , 1726K) is com- 
puted. Here, the solute ele\ient (Al) is comparatively high'y evaporative relative to 
the solvent (Ni). Equation (11) therefore applies, and the time, tj, to reach a final 
solute concentration m from a specified initial concentration m Q is directly pro- 
portional to N A/U (in the G constant). Table 1-2 gives the times to reach various 
° 1 -2 

final Al concentrations for one mole (5.367 x 10 Kg) of the 8% Al in Ni alloy 

(m = 0.159) evaporating at 1726K (the melting point of Ni) from its (supposedly 

0 -3 2 

constant) 10 m surface. 
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TABLE 1-2 NORMAL EVAPORATION OF 8^ BY WEIGHT ALUMINUM 
I IN NICKEL AT 1726K 


Evaporation Time, 
Sec 

Final Solute Concentration 

Alloy Remaining 

Weight Fraction 

Mole Fraction 

Weight, Kg 

Mole 

0 . 

0.080 

0.1591 

5.366 x 10"® 

1.0000 

28.3 

0.078 

0.1555 

5.353 x 10 ® 

0.9955 

57.1 

0.076 

0.1518 

5.341 < 10 ® 

0.9910 

86.5 

0.074 

0.1481 

5.328 xlO " 

0.9865 

117 

0. 072 

0.1444 

5.315 x 10 ® 

0.9820 

147 

0.070 

0.1407 

5.303 x 10 ® 

0.9775 

179 

0.068 

0.1370 

5.290 v i.o “ 

0.9731 

211 

0. 066 

0.1333 

5.2 " O'® 

0.9686 

244 

0. 064 

0.1295 

5.?'' 0 

0. 9642 

277 

0.062 

0.1258 

5.25 

0.9598 

312 

0. 060 

0.1220 

5.240 x 10 ® 

0.9554 

347 

0.058 

0.1182 

5.227 x 10"“ 

0.9510 

384 

0.056 

0.1143 

5.214 x lo"® 

0.9467 

421 

0. 054 

0.1105 

5.202 x 10"® 

0.9423 

460 

0. 052 

0.1066 

5.189 x 10"® 

0.9379 

500 

0. 050 

0.1028 

5.177 x 10~® 

0.9336 

542 

0.048 

0. 0989 

5.164 x 10~l 

0.9292 

584 

0. 046 

0. 0950 

5.152 x 10~® 

0.9249 

629 

0.044 

0.0910 

5.139 x 10'® 

0. 9206 

675 

0.042 

0.0871 

5.126 x 10"® 

0.9163 

723 

0.040 

0. 0831 

5.114 x 10"? 

0.9120 

774 

0. 038 

0. 0792 

5.101 x 10 _ ® 

0.9076 

827 

0.036 

0. 0752 

5.088 x 10'® 

0.9033 

882 

0.034 

0.0711 

5.075 x 10'® 

0.8990 

940 

0. 032 

0. 0671 

5.062 X 10'® 

0.8947 

1000 

0.030 

0.0631 

5.049 x 10"® 

0.8904 

1070 

0.028 

0. 0590 

5.036 x 10 _ ® 

0.8860 

1130 

0.026 

0.0549 

5.023 x 10‘‘ 

0.8817 

1210 

0.024 

0. 0508 

5.009 x 10 ® 

0.8773 

1290 

0.022 

0.0467 

4.996 x 10'® 

0.8729 

1380 

0.020 

0.0425 

4.982 x 10 ® 

0.8685 

1480 

0. 018 

0. 0384 

4.967 x 10'® 

0.8640 

1590 

0.016 

0. 0342 

4.953 x 10 ® 

0. 8594 

1710 

0.014 

0. 0300 

4.938 x 10 ® 

0.8549 

1850 

0.012 

0. 0257 

4.922 x 10"® 

0.8501 

2010 

0.010 

0.0215 

4.905 x 10"® 

0.8453 

2210 

0.008 

0.0172 

4.887 x 10 _ ® 

0.8401 

2470 

0.006 

0.0130 

4.866 x 10“® 

0.8346 

2820 

0.004 

0.0087 

4.841 x 10 _ ® 

0.8284 

3430 

0.002 

0.0043 

4.805 x 10'® 

0.8204 


i 


m 
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Thus, it takes 312 seconds for 5.37 x 10 Kg (1 mole) of an 8 % A1 in Ni to 

-3 2 -3 

change to 6 % in A1 concentration (A = 10 m ), with a decrease of 1. 26 x 10 Kg 

(2.33%) in the weight of the evaporating alloy. The same evaporating Ni-Al alloy takes 

an additional 411 seconds to change its aluminum concentration to 4% with an additional 

loss of 4. 71% in alloy weight. An additional time of 657 seconds is required for the same 

alloy to change in aluminum conct tration by another 2 % (to 2 %) with an additional loss 

of 13. 15% in the alloy weight. If the evaporating alloy is spread out as a thin layer with 

-2 2 -3 2 

an evaporating area A = 10 m instead of 10 m , all the evaporating times given in 
Table 1-2 must be reduced by a factor of 10, in accordance with Eq 11 or 11a. 

The alloy evaporating rate is directly proportional to the evaporating area A 

and gene, ally must, therefore, continuously decrease with shrinking alloy volume, mass, 

or surface. If the geometric shape of the evaporating alloy .mple is fixed and known, 

the time to reach a certain sample . ize or deposit composition can oe computed by 

numerical techniques based on Eq 6-9. Molten materials in zero-gravity environments 

readily assume spherical shapes due to surface tension. An example of such spherical 

-3 

evaporation at the melting point of pure nickel, i.e. , 1726K, for a 6.35 x 10 m 
(£ in. ) sphere with a composition containing 8 weight percent of aluminum in nickel is 
depicted in Figure 1-1. In this figure, curve 1 represents the weight fraction of the 
remaining evaporating sphere (WFR); curve 2 the diameter of the remaining sphere (DM); 
curve 3 the weight fraction of aluminum in the remaining sphere (WFA); and curve 4 the 
weight fraction of aluminum in the vapor (WFV). All four parameters are plotted against 
the evaporation time in seconds. 

Figures 1-2 and 1-3 show the normal evaporation behavior at 1726K of several nickel 
alloys containing 10% by weight of such highly evaporative solutes as Cd, K, Na, and Rb, 
these elements being much more volatile than Al. Figure 2 gives the weight fraction 
remaining in grams upon reaching specified solute concentrations in weight fractions, or 
weight percentages from the original 10 weight percent and decreasing successively by 

0.5%. Figure 3 gives the evaporating times (t. in Eq 11) to reach the same successive 

1 -3 2 

solute concentrations. Again, the evaporating area A is assumed constant at 10 m . 

It can be seen that all these fast-evaporating elements reach the different final solute 

concentrations five orders of magnitude faster than Al. Specifically, they take only about 
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500 microseconds, vs 58.3 seconds for Al, for the solute concentration to change from 
the initial 10 weight percent to 9.5 weight percent. Very small droplets or thin film;- 
having large evaporating areas relative to their masses would, therefore, change com- 
positions significantly, so rapidly that the change can never be prevented. The behavio- 
of any subsequent solidification may thus be completely unexpected, and accurate 
solidification analysis may thus be rendered impossible or meaningless unless the exact 
thermal history of the alloy is first very accurately measured and an evaporation study 
then undertaken to determine the magnitude of the surface composition changes in these 
alloys. 

Comprehensive computations have been made for nickel alloys containing initial 
solute concentrations of from 2 to 10 weight percent (at intervals) of 50 different 
solute elements evaporated at 1726K, the melting point of nickel. For each solute 
element, the evaporating times were computed for 0. 1 % decrements in the solute con- 
centrations until the final concentration of 1 % was reached. The final molar and weight 
fractions of the alloy remaining at the end of each evaporation period were also computed. 
These computations, designed to furnish a handy reference for future experiment planning, 
are available on request. In our original proposal, we recommended for study nickel 
alloys containing 2 or 8 weight percent of tin or aluminum as the specimen materials for 
our proposed experiment. 

2.3.2 Iron Alloys 

The evaporation behavior of binary iron alloys containing 20 different solute 
elements has also been studied. Table 1-3, listing the equi-evaporative temperatures for 
these 20 different alloy systems, gives the ratios of the solute evaporating rate, U, at 
1873K, to that of the solvent iron, V. At 1873K, 10 of these solute elements evaporate 

faster than the solvent (the three left columns) and 10 slower (the three right columns). 

-12 

Moreover, these ratios vary widely over 18 decades, from 4.41 x 10 for the slowest 

0 

evaporating, W, to 3.76 x 10 for the fastest evaporating, Cd. Because of this wide 
variation, and due to the extreme sensitivities of the evaporating surface to unsuspected 
contaminants, predicted or experimental evaporating results cannot generally be very 
accurate. 
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TABLE 1-3 EQUI-EVAPORATIVE TEMPERATURES 
FOR TWENTY BINARY IRON ALLOYS 


Solute 

U/V at 
1873 K 

T s’ K 

Solute 

U/V at 
1873K 

T s , K 

Cd 

3.765 x 10 6 

13,813 

W 

4.412 x 10" 12 

-69,317 

Zn 

2.065 x 10 

15,883 

c 

8.210 x 10~ 8 

5 o'sj 

Mg 

9.074 x 10° 

19,383 

Mo 

1.879 x 10“ 7 

-11,857 

Ca 

5.503 x 10 4 

9,608 

Zr 

1.960 x 10" 6 

-*7,517 

Pb 

3.333 x 10 3 

5,249 

B 

6.83*: x 10‘ 5 

9,783 

Mn 

1.070 x 10 3 

21,013 

V 

3.804 x 10" 3 

8,847 

Al 

45.65 

8,291 

Ti 

2.254 x 10~ 2 

34,923 

Cu 

12.25 

5,886 

Co 

4.362 x 10 _1 

4,591 

Sn 

9.978 

2,959 

Si 

5.238 x 10' 1 

3,098 

Cr 

3.158 

58.1 

Ni 

5.897 x IQ' 1 

3,309 


The equi-evaporative temperatures in iron alloys also v.iry widely. Binary iron 
alloys containing Cr, Zr. Mo, and W have no practical equi- evaporative temperatures. 
One can, therefore, always expect these alloys to change compositions continuously with 
the evaporating time. 

Table 1-4 shows the effect of evaporating temperature on the U/V ratios for four 
different solute elements Mg, Ca, Mn, and Al. In the range of 1773K to 2173K and 
beyond, increasing the evaporating temperatures always decreases the U/V ratios. 

This can also be seen from Table 1-3, as the equi-evaporative temperatures shown for 
these four binary iron alloys are higher than 1873K, at which temperature the four solute 
elements evaporate much faster than the solvent iron. Table 1-4 also shows that, for 
the same temperature variation, the more evaporative the solute element the greater the 
percentage variation in the U/V ratio. In the case of Mg the U/V ratio decreases by 
about 20 times from 1773K > 2173K, whereas for Al the same ratio decreases only by 
less than three times over the same temperature interval. 


1-18 




TABLE 1-4 EFFECT OF SOLUTE ELEMENTS AND EVAPORATING TEMPERATURES 
ON THE VALUES OF U/V (x 1000) IN IRON ALLOYS 


Temp (K) 

Solutes 

Mg 

Ca 

Mn 

A1 

1773 

2132.0 

118.2 

1.649 

0. 06037 

1873 

907.4 

55. 03 

1.070 

0.04565 

1973 

421.2 

27.68 

0. 7261 

0. 03552 

2073 

210.5 

14.88 

0.5110 

0. 02831 

2173 

112.1 

8.466 

0.3719 

0.02304 


Figure 1-4 shows the effect of solute elements and evaporating temperature, T, 
on the evaporating time t, for a given set of initial and final solute concentrations (i.e. , 
m Q = 0. 01 and m = 1 ppm). For given m Q and m, the log t versus 1/T curves for these 
highly evaporative solute elements are approximately linear and have positive slopes. 
This can be predicted from Eq (11) since, for given m Q and m, log t is linearly related 
to B/T - 0.5 log T, and 0.5 log T is small relative to B/T within the evaporating tem- 
perature range studied. Thus, one can use Fig. 1-4 to determine the value of the ele- 
mental evaporating constant B, or heat of evaporation aH = 4.574 B (Dushman, 1962), 
for the solute elements by plotting log t versus 1/T and measuring the slope of the 
resultant, nearly straight lines. 

In Fig. 1-5, the log t versus 1/T relationships also appear nearly linear for all 

-1 -2 -3 -4 

the four different initial concentrations (i.e. , m = 10 ,10 ,10 , and 10 ) of A1 in 

Fe. This can also be seen from Eq (11) and (12). All these nearly straight lines have 
identical slopes, from which the heat of evaporation of pure A1 can be evaluated. 

Figures 1-6 and 1-7 display the effect on the evaporating time, t, of final solute 
concentration, m, and either initial concentration m Q (in Fig. 1-6) or solute elements 
(in Fig. 1-7). 

Figure 1-8 indicates that Mg, Ca, and Mn in Fe alloys are so evaporative at 
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Fig. 1-4 Evaporating Times for Different Dilute Iron Alloys to Change 
Concentration from = 0.01 to m = 1 ppm 
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Fig. 1-5 Evaporating Times for Dilute Iron Alloys Containing A1 and 
Starting at Four Different Initial Concentrations , to Reach 
a Final Concentration of m = 1 ppm 
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Fig. 1-6 Effect of Initial and Final Solute Concentrations on the Evaporating 
Times in Dilute Iron Alloys Containing A1 at 1873K 
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Fig. 1-7 Effect of Solute Elements and Final Concentration on the 
Evaporating Times in Dilute Iron Alloys at 1873 K 
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Fig. 1-8 Evaporation Loss of Dilute Iron Alloys at 1S73K as a 

Function of Solute Elements and Final Solute Concentrations 
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1873K that practically all of these elements are removed by evaporation, without much 

evaporative loss of the solvent Fe atoms. On the other hand, A1 is comparatively less 

evaporative, and much of the solvent Fe atoms are evaporated off together with Al. To 

7 -9 

achieve a purification factor of 10 (i. e. , to m = 10 ) from m Q = 0. 01, for example, 
the initial evaporating alloy must lose over 30% of its material. 

2.4 ACCURACY OF PREDICTED RESULTS 

To check the validity of our derived equations for normal evaporation of binary 
alloys, we searched the current literature. Several sets of alloy evaporation data have 
been found that are amenable to normal evaporation analysis. 

2.4.1 The Fe-Ni and Ni-Cr Systems 

An analysis has been made of the data by Obradovic et al. (1969) for Fe-Ni and 

Ni-Cr alloys evaporated at 1873K for different times under various ambient pressures. 

In the 80% Ni-20% Cr case, evaporation started with N = 0. 7582 moles of the alloy having 

-3 2 ° 

an evaporating area A = 10 m . The solute and solvent evaporating rates are, respec- 

-1 -2 9 

tively, U_ = 3.479 x 10 and V... = 6.3»ti x 10 mol/m /sec. In this case the solute, 
Cr Ni 

Cr, is comparatively highly evaporative relative to the solvent, Ni. Hence, the evapor- 
ating time t^ given in Eq (11) applies, i.e. , 

[ (1-m )m m-m 

O 0 

jfn :: ~ + ~z “ - 

m (1-m) (1-m ) (1-m) 

o o 

where G„ = -N (1-m )/AU, and m and m are the initial and final molar solute con- 
1 o o o 

centrations in the alloy. 

Least-square fits of the Obradovic data give nearly constant values of observed 

4 

G ; -1. 785, -5. 064, and -Z. 064 (x 10 sec) for ambient pressures of 1, 100, and 500 

* _3 

(x 10 torr), respectively. The calculated values of G^, though also nearly constant, 
are, however, one order of magnitude smaller, indicating surface solute depletion. 
Similar analysis for the 60% Fe-40% Ni alloy evaporated with N q = 0. 7478 moles and 
A = 10 3 m 2 gives U Ni = 6.386 x 10 2 and V Fe = 1. 081 x 10 1 mol/m 2 /sec. Here, the 
solvent Fe is much more evaporative than the solute Ni, and Eq (14) applies. The 
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values of observed evaporating coefficients G = m N /AV are 2.469, 1.369, and 
„ 2 o o 

4 2 

8.386 (x 10 sec) for ambient pressures of 0.133, 13.3, and 66.7 Newtons/m , 

respectively. The calculated G 's are again nearly constant but also an order of 

dt 

magnitude smaller than the observed G^s, again indicating solvent depletion (or 
solute enrichment) at the surface. Details of these analyses are given in Appendix A. 

2.4.2 Beryllium Purification 

The kinetics of normal evaporation for beryllium have been quantitatively 

checked with actual results of beryllium purification during vacuum induction melting. 

Details are given in Appendix B. In these tests by Bunshah and Juntz (1966), beryllium 

-4 

was crucible-free, induction-melted under an ambient pressure of 1.33 x 10 
2 

Newton/m . The actual temperature of the melt was not known, and the exact evaporating 
rates for the beryllium solvent and various solute elements cannot be computed. How- 
ever, because all solute elements (i.e., Fe, Cr, Mn, Ni, Si, Al, Mg, Cu, Zn, and Na) 
and the solvent, Be, were evaporating from the same or common liquid-gas interface 
of a fixed area for the same length of time, we can compute the values of P, defined as 
the product of the evaporating time, t, and solvent evaporating rate, V (for V » U), or 
solute evaporating rate, U (for U » V); thus, P = tV or tU. Table 1-5 gives the initial 
concentrations and final concentrations (in ppm) of the various solutes, together with the 
value of P, actual surface concentrations, ratio of actual surface to bulk concentrations, 
and effective times to reach the final concentrations under the assumption that the eva- 
porating temperature was 1523K. The following conclusions can be drawn from Table 5-1: 

• The solute elements can be divided into three groups: Fe, Cr, Ni, Si, 
and Cu evaporate much more slowly than Be; Mg, Zn, and Na evaporate 
much more rapidly than Be; and Al evaporates at about the same rate as 
Be 

• The computed times to reach the final concentrations under normal evapor- 
ation conditions (complete liquid mixing) are fairly constant for the solute 

elements that evaporate much slower than Be (i.e., Fe, Cr, Ni, Si, and 

4 

Cu), being about 4x10 
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• The computed times to reach the nnai concentrations for the highly 
evaporative Mg, Zn, and Na are also fairly constant, but about four 
orders of magnitude smaller 

• After correction for limited liquid mixing (A*w>nrRx B), the effective 
times to reach the final concentrations are much more constant, even 
between the groups of solute elements. In particular, the highly 
evaporative elements Mg, Zn, and Na have their results improved by 
several orders of magnitude 

• Limited liquid diffusion must, therefore, be considered, particularly for 
the highly evaporative elements 

• The surface solute concentrations can be changed by up to six orders of 
magnitude so that the effective evaporating rates can be similarly 
changed. 

2.5 NORMAL EVAPORATION VS NORMAL FREEZING 

Normal evaporation and normal freezing (Pfann 1966) are strikingly sim.lar in 
many respects. Both deal with idealized transition phenomena between two phases. In 
normal evaporation, the two phases are either vapor and liquid or vapor and solid. In 
normal freezing, the two phases are always solid and liquid. In either case, the interface 
between these two phases moves in accordance with the constraints of mass transfer, 
heat transfer, and liberation or absorption of the heat of evaporation or fusion. Either 
phenomenon is idealized under the three assumptions tabulated below. 


TABLE 1-6 COMPARISON BETWEEN ASSUMPTIONS IN NORMAL 
EVAPORATIONS AND NORMAL FREEZING 


Normal Evaporation 

Normal Freezing 

Vapor or solid: complete vapor removal, i.e., D y = ® 
Liquid: complete liquid mixing, Dj = ® 

Phase Equilibrium: Raoult's Law 

Zero solid diffusion, D = 0 

s 

Complete liquid mixing, ® 

Phase diagram 


1 

1 
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In practice, complete vapor removal in normal evaporation has the same effect 
as no solid diffusion in normal freezing. Both assume that the material changed into the 
vapor phase (in normal evaporation) or solid phase (in normal freezing) is permanently 
and completely lost, and no longer has anything to do with the remaining matenal, as 
far as mass transfer is concerned. Complete liquid mixing is assumed in both cases. 
From the phase diagram, the freezing segregation coefficient k can be determined; bastd 
on Raoult's Law, the evaporative segregation coefficient can be obtained. Both segre- 
gation coefficients are generally functions of the solute concentration and phase transition 
temperature, but can be considered constant under certain special conditions. Specifi- 
cally, in extremely dilute solutions with limited solidifcation range, k is nearly constant; 
while if U <s: V and m « 1, k^ is constant at a given evaporation temperature, as in- 
dicated above. 

The value of k in normal freezing generally is within the relatively narrow range 
-3 

of 10 to 10. On the other hand, k in normal evaporation of nickel alloys at 1726K 

-12 V 9 

varies from 1.38 x 10 for W to over 10 for such highly evaporative solute elements 
as K, Rb, and Cd as shown above. The values of k^ are thus many orders of magnitude 
wider in range than those of k. 

The velocity v of interfacial movement has a much broader range in normal 

—6 _3 

evaporation than in normal freezing. Normally, v varies from 10 to 10 m/sec, only 

three orders magnitude from maximum to minimum. In normal evaporation of nickel 

-20 

alloys at 1726K, on the other hand, v varies from 6.74 x 10 m/sec for W to about 
5. 0 m/sec for K and Rb. The spread is almost 20 orders of magnitude. 

Among the three assumptions (Table 1-6) in normal evaporation or freezing, the 
one on complete liquid diffusion is most subject to error. Much work has been done 
to include the effect of limited liquid diffusion in normal freezing. The equa f! ons of 
Tiller et al (1953) consider the initial transient state, the subsequent steady state, and 
the final stage of solidification, because of solute pile-up in front of the solid-liquid inter- 
face due to limited liquid diffusion. These equations are not very accurate for the general 
case (Pohl, 1954), except when k is nearly zero. Fortunately, for the important case of 
surface enrichment of trace impurities of a highly non-volatile impurity in solution, the 
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evaporative segregation coefficient at a given evaporation temperature is k^ = U/V, 
which is indeed close to zero. Hence, Tiller's equation for limited liquid diffusion can 
be bodily transferred to the normal evaporation study with limited liquid diffusion. 

We thus see that with limited liquid diffusion before the steady state is reached, 
there exists a transient state during which the solute concentration and thermophysical 
properties of the surface-enriched layer change constantly in a manner that can be, but 
apparently never has been predicted. 

For the special case of U «V and m « 1 involved in space experimentation, it is 
necessary to know the "characteristic distance" and "decay constant" (across and over 
which the relative solute concentration changes by a factor of e, as defined by Tiller 
et al (1953)). Such determinations allow us to see how much unavoidable evaporation affects 
solidification during a given test period, or indeed to determine if the solidification experi- 
ment is meaningful or definitive and not unduly distorted by surface evaporation. 

The characteristic distance is, as defined by Tiller et al (1953), D /kv, where 

-8 2 ^ 

D is the liquid diffusion coefficient (generally about 10 m /sec). This distance for 

normal freezing thus covers about seven orders of magnitude, from 10 to 10. 0 m. 

The characteristic distance for normal evaporation has a much broader range, covering 

-18 23 

over 41 orders of magnitude and varying from 10 to 10 m. 

2 

The decay constant is, as defined by Tiller et al (1953), D./kv and for normal 

-3 ^ 7 

freezing thus covers about 10 orders of magnitude, from 10 to 10 sec. The decay 

constant for normal evaporation also has a much broader range, covering over 71 orders 

-28 43 

of magnitude, from 10 to 10 sec. 


These much greater evaporative segregation effects would almost certainly conceal 
any minor or subtle zero-gravity effects, especially in the presence of other unknown or un- 
controlled effects. Definitive space solidification work should, therefore, be preceded by 
an evaporative compatibility study of the sample materials and associated impurities. In 
fact, evaporation is almost certain to be very important or so overwhelming that the effect 
of zero-gravity or freezing segregation may be masked or even reversed. A freely sus- 
pended molten drop in space may, for example, have its surface solute concentration 
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greatly enriched, even as much as a millionfold by neglected and undetectable trace 
impurities within a fraction of a second of its deployment. We are then dealing at the 
critical surface layer with a completely new and unanticipated alloy, e.g. , trace impurity 
indium instead of the original zinc, having entirely different composition, melting point, 
surface tension, thermophysical properties, latent heat of fusion, undercooling and 
nucleation characteristics, growth morphology, and the like. 

These and similar analyses suggest that the normal evaporation equations pre- 
sented above work surprisingly well for the several alloy systems studied so far. The 
equations generally give results accurate to within an order of magnitude of the correct 
values, perhaps even for highly concentrated solutions (e.g. , 209? Cr in Ni or 40 9? Ni in 
Fe). Such accuracies are sufficient for many evaporation studies. 




In more critical studies, however, refined evaporation analyses may be needed. 

In particular, the effect of limited diffusion and, hence, surface depletion or enrichment 
of the solute must often be accounted for. This is true especially for the highly evaporative 
solutes such as Mg, Zn, and Na in the solvent Be. We already have results that confirm 
this conclusion. Details of these results will be published in Ref. 8 and elsewhere. 

Section 3. ACCELERATIONS DUE TO DIFFERENTIAL EVAPORATION 


In any zero-gravity experiment, it is important not to subject the test samples to 
accelerations due to gravity or other forces. Further, the molten samples should not 
touch the wall of the test chamber before solidification. Different rates of evaporation 
at various surface regions of a freely suspended sample do, however, give rise to un- 
balanced forces and momenta that may produce erratic or unwanted accelerations and 
undesirable sample impingement onto the test chamber wall. Such differential evaporation 
may result from solute segregation, nozzle detr.rbment, or temperature grad ; rnt. For 
example, one-sided heating by an electron beam, introduces severe e* ^porative segre- 
gation and very steep temperature gradients. Previous work \y Dr. T. Kassal at 
Grumman Research on water jetting into vacuum indj£t..£s that differential evaporation 
introduces unpredictable movements that are,^^ily observed in motion pictures. 


When planning space experiments on material solidification it is, therefore, 
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desirable to have quantitative answers to .e following specific questions: 

• What are the jetting forces on samples of different elemental materials 
under given temperature nonuniformity conditions ? What elements are 
stable (i.e. , subjected to small jetting forces) at their melting points, 
and what elements are unstable ? 

• What is the effect of sample size and density on jetting forces ? 

• For a given sample material, what are the critical or undesirable im- 
purities, and what are the desirable or stabilizing impurities ? 

• What are the equivalent accelerations or gravities on the samples under 
specified test conditions ? 

• Starting from rest, how long does it take for a differentially evaporating 
sample to reach the wall of test chamber of a given size ? 

• What is the minimum test chamber size or maximum degree of super- 
heating for a given sample material and size? 

• How should the sample be heated and released, and with what equipment ? 

3.1 ANALYSIS OF THE JETTING ACTION 

The above questions can be answered by a simple analysis of the jetting action on 
a differentially heated sphere, such as is shown in Fig. 1-9. The right half of the sphere 
is assumed to have a surface temperature T , which is lower than T , the surface 

•L b 

temperature of the left half of the sphere. As a result, the vapor pressure p on the 

(U 

left half is greater than that (p^) on the right half. There is then a net force f which 
produces an acceleration a. This acceleration causes the sphere (diameter d, density p ) 
to travel a distance s to impact the wall of the test chamber or other nearby objects. 

The jeiMng force due to differential evaporation is (see Appendix C) 

'^(f d2 ) lp 2- p 1 l = ('F d3 =) “ 
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Hence, the acceleration due to the jetting action is: 

a = (Pg-P^dp 

The equivalent gravity is then a divided by the gravitational constant g = 9.81 
m/sec^, 

g g = a/g = a/9. 81 

And the time t for the jetting sample to travel the distance s before impacting 
the nearly wall or object is 

t =^2s/a -^8dsp /3 (p 2 - Pl ) 

In the abo'^e analysis, we assumed zero initial velocity and acceleration; the 
equations are easily corrected for nonzero initial velocity and acceleration. The times 
computed are, therefore, the maximum allowables for performing the experiment or pro- 
cessing. Note that these times are proportional to pds raised to the one-half power. In 
addition, they depend even more critically on AT = T -T and on the evaporation behavior 

M -L 

of the sample material, i.e. , on p 9 -p^. This is because the elemental vapor pressures 
vary by over 38 orders of magnitude at the melting point of nickel (1726K), and because 
these vapor pressures are highly temperature sensitive, as previously shown. 

3.2 COMPUTED JETTING ACTION ON SINGLE ELEMENTS 

The 51 metallic elements listed in Table 1-1 have been screened for potential, 

_3 

space-processing sample materials. Here, the sample diameter is fixed at 6 x 10 m, 

while the allowable distance of travel is constant at 0.3 m. Further, half of the spherical 

surface is assumed to be at T^, the melting point of the sample material, while the other 

half is at T , being superheated by 5 to 100% of the absolute melting temperature. Com- 

putations show that the allowable times at 5% superheating vary widely from tens of milli- 

g 

seconds to over 10 years for gallium. Table 1-7 summarizes the data on these highly 
stable sample materials, as compared to nickel. These results further stress the 
critical importance of evaporative segregation. 

Figures 1-10 through 1-12 show the jetting characteristics of the top four, or most 
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stable elements, i.e., Ga, Sn, In, and La. These figures illustrate the superheated 

-3 

vapor pressure, equivalent gravity, and time to travel 0. 30 m, for a 6. 0 x 10 m molten 
metal sphere, heated to its melting point on one side but superheated on the other. The 
amount of superheating varies from 5 to 100% at 5% intervals. Figures 1-13 and 1-14 
give similar jetting data, but for nickel and iron samples respectively. 

Thus, the unwanted accelerations due to differential evaporation, while completely 

negligible in many materials, are serious factors to be considered in others, including 

the specified nickel samples in M512 Skylab experiments. Specifically, it takes only 

-3 

1 0 or 20% superheating on one side of a 6 x 10 m molten nickel sphere to impact the 
experimental chamber wall 0.3 m away within three or one seconds, respectively, of its 
release. Complete solidification in free suspension may, therefore, be very difficult 
under these conditions. 

_3 

In other words, nonuniformly heated, 6x10 m spheres are subjected to equi- 

-39 -3 

valent gravities of from 2.26 x 10 Kg for gallium to 3. 25 x 10 Kg for magnesium, if 
the spheres are superheated by only 5% on one side relative to the other side at its 
absolute melting point. From this, we anticipate new types of powerful evaporative con- 
vection currents to be described later. 

Notice that, in both nickel and iron spheres of the specified size, the computed 
equivalent gravity exceeds the one-g level if one side is superheated by 35 and 267c, 
respectively, relative to the other side at its absolute melting point. Such evaporation 
and jetting action is, of course, merely the familiar rocket action. In the performance of 
some Skylab experiments, however, such differential jetting can have important and 

_3 

serious implications. The 6x10 m spheres in the M512 sphere-forming experiment, 
for example, are known to be very unevenly heated by the electron beam from one side, 
to the point of occasional spitting or explosion, but remain solid on the other side. Such 
samples may be subjected to very large, unwanted accelerations (and hence equivalent 
gravities) of unpredictable magnitude, direction, and duration. Unless corrective action 
is taken, such as changing the sample material, improving the hardware design (now 
considered too late), or optimized heating procedures, any solidification results may be 
difficult to interpret in terms of zero-gravity effects. 
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Fig. 1-13 Jetting Action on 6x10 m Nickel Sphere, One Side at the Absolute 
Melting Point of 1720K while the Other Side at Specified Degree of 
Superheating 
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Figures 1-15 and 1-16 display the jetting characteristics on partially superheated, 
6x10 m spheres of highly volatile materials due to differential evaporation. These 
figures are, respectively, for equivalent gravity and time to travel 0. 3 rn, for seven 
different volatile elements at eight different degrees of superheating from 5 to 100%. 

3.3 COMPUTED JETTING ACTION OF IMPURITIES IN NICKEL 

The undesirable jetting action due to certain evaporating impurities can also be 

replaced by the desirable stabilizing action of other impurities. Specifically, a 
-3 

6. 0 x 16 -m nickel sphere, at its melting point of 1726K, may have its normal jetting 

-3 -15 

action at 5% superheating on one side reduced from 2.15 x 10 g to 7.34 x 10 g by 

-14 

an enriched surface layer of tungsten, or to 6. 75 x 10 g by tantalum. These data on 
the stabilizing impurities on nickel are illustrated in Fig. 1-17. 

Figures 1-18 to 1-23 show the stabilizing actions of surface layers of the six most 
powerful stabilizers in nickel. These stabilizers, as shown in Fig. 1-17, are W, Ta, 

Cb, Os, C, and Mo. These figures show, as ordinates, the superheated vapor pressure, 

_3 

equivalent gravity, and time to travel 0.3 m, for a 6 x 10 m molten nickel sphere. 

This sphere has a surface layer of the specified stabilizing element and is heated to the 
melting point of nickel (1726K) on one side but superheated on the other. The amount of 
superheating in percentage points is plotted as the abscissa. 

Figures 1-24 and 1-25 give the jetting characteristics on partially superheated, 

_3 

6x10 m nickel spheres, surface-contaminated with powerful agitators or highly 
volatile elements. These figures are, respectively, for equivalent gravity and time to 
travel 0.3 m for several volatile elements at different degrees of superheating up to 100%. 

_3 

If we compare Fig. 1-18 with Fig. 1-24 or 1-25 we see that, on a 6 x 10 m 

nickel sphere superheated above its absolute melting point by 5% on one side, the equi- 

7 -15 

valent gravity varies by 19 orders of magnitude from 4. 87 x 10 for K to 7. 34 x 10 
for W. For the same partially superheated nickel sphere, the time to travel 0.3 m varies 
by 10 orders of magnitude froml.84 x 10 4 sec (i.e., 184 usee) for K to 2.89 x 10® sec 
(i.e. , 0.92 years) for W surface layer. 
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Fig. 1-16 Time for a Partially Superheated, 6x10 m Sphere of Highly Volatile 
Material to Travel 0.3m Due to Differential Evaporation and Jetting 
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Fig. 1-24 Eqivalent Gravity on 6x10 m Partially Superheated Nickel Sphere with 
Surface Layer of Highly Volatile Element Due to Differential 
Evaporation and Jetting (Sheet 2 of 2) 
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3.4 FLUID DISTURBANCES DUE TO JETTING FROM DIFFERENTIAL EVAPORATION 


Two of the great advantages of manufacturing in space arc zero-gravity conditions 
and the complete absence of convection currents due to the earth's gravity. With the 
absence of all convection currents, localized fluid disturbances, temperature differences, 
concentration variations, and complex heat and mass transfer conditions do not exist. 
Solidification can, therefore, be more easily controlled, and defect-free crystals are 
then readily obtained. 

However, even in zero gravity, localized or interfacial fluid flow, movement, 
instability, or other disturbance can still occur due to surface tens i on gradients resulting 
from temperature or concentration gradients. This phenomenon, generally called the 
Marangoni effect or instability tBikerman 1958, Scriver, & Sterling 1960, Sterling & 
Scriven 1959) manifests itself in such forms as "interfacial engine", instability, 
twitching, cellular motion, Benard's cells, and tear drop in a wine glass. 

Another type of fluid disturbance not previously reported may result > m jetting 
or rocket action due to differential evaporation as described above. A nonuniformly 
heated liquid drop in space will be non-spherical because both surface tension and jetting 
action are strongly dependent on temperature. At equilibrium, the internal pressure of 
the drop is, of course, constant. The varying jetting pressure therefoie appears as a 
variation in the curvature of the drop so that the algebraic sum of the jetting pressure and 
the internal pressure due to surface tension is constant. Since the jetting pressure and 
surface tension vary according to entirely different functions of the drop temperature, 
such constancy can be maintained for only a selected uoint or ring on the drop surface. 
The overall shape of the liquid drop must, therefore, be non-spherical. 

Under gravity or the equivalent, jetting action may show up as local variations in 
the fluid level. Figure 1 -26(a) shows the heating characteristics of t liquid in a con- 
tainer. The instantaneous profiles for the temperature, evaporation rate, jetting 
pressure, surface tension, and surface curvature are shown. For a liquid being heated 
in a container, the sides are at higher temperatures and, hence, have higher evaporating 
rates and jetting pressure, but generally lower surface tension than the center region. 
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Fig. 1-26 Heating or Cooling Characteristics of Liquid in a Container 
Under s / Tie True or Equivalent Gravity. 
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The liquid surface must therefore convex upward if the jetting pressure is large. 

Large jetting pressures are often possible, as shown previously. Note that the effect 
of lowering surface tension with temperature partly compensates for that of the jetting 
action; the liquid surface is more planar than if surface tension were absent. Similarly, 
Fig. l-26(b) shows the cooling characteristics of the same liquid in the same container. 

Jetting action may also introduce a new type of powerful convection current. 

The surface of liquid drop is usually not clean. Certain spots or regions S on the surface 
may be contaminated with slow -evaporating material (Fig. l-27(a)). These spots or 
regions may be oxidized into relatively non-evaporating oxides. Even under uniform 
heating conditions, such a drop may be evaporating nonuniformly over its surface. 
Localized, differential evaporation or jetting forces may then develop to supply the 
necessary mechanical energy for initiating or maintaining convection currents, as 
shown in the Figure. 

Figure 1 -27(b) shows a similar, uniformly heated liquid drop that is surface- 
contaminated with highly-evaporative material at certain spots or regions S, resulting in 
higher jetting pressures p,. at S compared to neighboring regions. Such differential 

-4 

jetting forces ir. 'v, again, initiate powerful convection currents. It has been shown by 
Prof. Reid of MIT (1971) that even a few atomic layers of selected materials can increase 
the local evaporating rates by orders of magnitudes for sustained periods of time. 

Note that these convection currents, resulting from differential evaporation, are 
completely different in origin and characteristics from those induced by differential 
surface tension (Marangoni effect). In fact, the two types of convection currents often 
have counteracting effects, one tending to compensate for or cancel the other. Also, they 
differ from gravitational convection currents. Gravitational convection currents are 
generally very massive, extending deep into the body of the liquid. Convection currents 
due to differential surface tension, on the other hand, are b" nature surface limited. 

The new convection currents induced •' jetting due to differential evaporation extend to a 
degree intermediate between the ou.ex two types of convection currents. Solidification 
behavior predictions of a specified system based on results on one type of convection 
current may, therefore, not be valid when applied to other types of convection currents. 
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(b) LOCALIZED CONTAMINATION BY HIGHLY VOLATILE MATERIALS 


Fig. 1-27 Convection Currents Induced by Differential Evaporation 
and Jetting 
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Section 4. EQUI-EVAPORATIVE TEMPERATURES AND THEIR APPLICATIONS 


The normal evaporation techniques presented above give results accurate within 
an order of magnitude of the correct values, at least for the several systems studied. 
More refined methods, including those accounting for the limited liquid diffusion in the 
evaporating source, should be considered, as evidenced by the great improvement in the 
predicted results in the purification of beryllium containing such highly evaporative 
solutes as Mg, Zn, and Na. During evaporation, not only the solute concentration, but 
also the shape, size, and evaporating area of the evaporating source may be constantly 

_3 

changing. Figure 1-1 clearly shows the normal evaporation behavior of a 6 x 10 m 
sphere containing 8 weight percent of aluminum in nickel. 

In the most general case, we consider steady-state and transient-state evaporation 
with constant changes in solute concentration, shape, and size of the evaporating source. 
The liquid diffusion is usually far from complete, and Raoult's Law of ideal solution may 
not be valid. Not only mass diffusion, but also heat transfer, fluid dynamics, and other 
related phenomena must be dealt with. 

4.1 EQUI-EVAPORATIVE TEMPERATURES IN ALLOY SYSTEMS 

The evaporation of alloy samples likely to be met in space processing (zone 
melting, crystal growth, thin film formation or, in general, controlled solidification and 
heat treatment) is complicated. Fortunately, for a given solvent and solute there exists 
a unique equi -evaporative temperature at which the solvent and solute are equally evapor- 
ative, i.e., U = V. Hence, at this unique temperature the alloy remains constant in 
composition, independent of the initial solute concentration and sample shape or size. 

This temperature, which can be computed from Eq (21), may be non-existent because 
it is negative or it may be too high or too low to be practical. 

However, for many binary or even ternary alloy systems, the equi-evaporative 
temperatures do exist and are practical. A computer search has been conducted among 
the 2550 binary combinations (i.e. , 51 solvents x 50 solutes) and of 62,475 (i.e. , 

51 solvents x 50 x 49/2 first and second solutes) ternary combinations, of the 51 metallic 
elements listed in Table 1-1. By "practical" we mean that the equi-evaporative 
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temperature is no higher than 200% of the absolute melting point of the solvent. At each 
such equi-evaporative temperature, the relevant binary phase diagram in Hansen (1958) 
was consulted to determine the range of solute concentrations for which the alloy is in 
liquid form. 

The results of these computer searches are summarized in Tables 1-8, 1-9, and 
1-10. Table 1-8 gives the equi-evaporative temperatures for binary alloy systems con- 
taining liquid phases. Table 1-9 lists the same temperatures, but for other binary alloy 
systems which either contain no liquid phases at these unique temperatures or have ill- 
defined phase diagrams according to Hansen (1958). Table 1-10 gives the equi-evaporative 
temperatures for ternary alloy systems. With the exception of alloys of carbon, which 
generally have ill-defined phase diagrams, all other alloys having practical equi- 
evaporative temperatures are listed in these tables. 

Consider, for example, Table 1-8. The Ag-Ga system has a equi-evaporative 
temperature of 1245. 8K or 972. 7°C. This unique temperature is 101% of the absolute 
melting point of Ag (1233.9 K) but over 411% of the absolute melting point of Ga (302. 9 K 
or 29.81°C. Hence, Ag alloys containing Ga having a "practical" equi-evaporative 
temperature of 1246K are listed in Table 1-8; but Ga alloys containing Ag have an 
"impractically" high equi-evaporative temperature (also at 1246K) and are not listed. 
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TABLE 1-8 EQUI-EVAPORATIVE TEMPERATURES FOR TWO-COMPONENT 
SYSTEMS CONTAINING LIQUID PHASES 
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TABLE 1-9 EQUI-EVAPORATIVE TEMPERATURES FOR 
OTHER TWO-COMPONENT SYSTEMS 
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Solvent 

Melting 
Point, K 

Solute and TsK 

Au 

1336 

Sc 2670 

B 

2573 

La 4333; Rh 3069; Pt 2003; Th 970. 7; V 2933 

Ba 

983.1 

Tl 704.8 

Be 

1557 

Ga 2519; In 2994; Sn 946.8 

Ca 

1116 

Sb 671.3; Sr 1633 

Cb 

2741 

Mo 3489; Os 3634; Zr 4980 

Cd 

594.1 

Cs 1139; Rb 1179 

Ce 

1077 

Au 1310; Cu 1940; Ge 1097; Sc 985. 0; Sn 1852 

Co 

1765 

Au 3293; Ge 2975; Pd 1335; Si 1050; Sn 3164 

Cr 

2173 

A1 4030; Au 1807; Cu 2951; Ga 3865; Ge 1513; 
In 3996; Pb 4166; Sc 1477; Sn 2291 

Fe 

1809 

Au 2902; Ge 2465; Ni 3309; Pd 324.8; Sc 3091; 
Si 3099; Sn 2959 

Ge 

1210 

Ce 1097; Cr 1513; Sc 1617 

Ir 

2716 

Mo 297.1; Ru 4513; Th 4722; U 5188; Zr 1678 

La 

1193 

Si 306.7; Ti 2104; Y 1208 

Li 

453.1 

Sr 288.9 

Mg 

923.1 

Cs 1834 

Mn 

1517 

Pb 2430 

Mo 

2893 

Cb 3489; IR 297.1; Os 3414; Ru 1596 

Ni 

1726 

Au 3110; Fe 3309; Ge 2798; Pd 1137; Sc 3178; 
Si 2614; Sn 3111 

Os 

3318 

Cb 3634; Mo 3414; Zr 6125 

Pd 

1825 

Co 1335; Fe 324.7; Ni 1137; Si 10470; Y 3277 

Pt 

2042 

B 2003; Ru 801.5; U 4038 

Rh 

2233 

B 3069; Pt 801.5; U 2878 

Ru 

2583 

Ir 4513; Mo 1596; Th 4619; V 5005; Zr 3085 

Sb 

903.6 

Ca 671.3 

Sc 

1811 

Au 2670; Ce 985. 0; Cr 1477; Fe 3091; Ge 1617; 
Ni 3178; Si 3095; Sn 2912 

Si 

1685 

Au 2969; Co 1010; Fe 3099; Ge 2689; La 306.7; 
Ni 2614; Pd 1293; Sc 2095; Sn 2987 

bn 

505.0 

Be 946.8 

Sr 

1043 

Ca 1633; Li 288.9 

Th 

2027 

B 970. 7 

Ti 

1940 

La 2104 

T1 

576.1 

Ba 704.8 

U 

1403 

V 1328 

V 

2193 

La 2885; Pd 4143; U 1328 

Y 

1793 

La 1208; Pd 3277 

Zr 

2133 

Ir 1678; Ru 3084 
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TABLE 1-10 EQUI-EVAPORATIVE TEMPERATURES FOR 
THREE-COMPONENT SYSTEMS 


Solvent 

Melting 
Point, K 

Solute 

Ts, K 

Liquid 

Concentration 

Range 

Cd 

321. 0 

K 

1180 

0 - 1.00 



Rb 

1179 


Si 

1412.0 

Sc 

3095 




Fe 

3099 


Ni 

1453.0 

Au 

3010 




Sn 

3011 



4.2 SIGNIFICANCE OF THE EQUI-EVAPORATIVE TEMPERATURE 

The binary or ternary alloy systems presented in Tables 1-8 through 1-10 are 
ideal, at their respective unique temperatures, for the study of heat transfer, sphere 
detachment, deformation and oscillation, bubble formation and movement inside the 
liquid, and the like, since their solute concentration is constant regardless of the initial 
concentration and sample shape or size. In addition, by changing the initial solute con- 
centration we can isolate and study the effect on these phenomena of such concentration- 
related properties as surface tension. In all such experiments, the solute concentration 
will remain constant throughout the entire sample during the experiment. There is no 
concentration nonuniformity and no surface concentration gradient. The initial, boundary, 
intermediate, and final conditions of the sample composition are completely defined, 
known, and stable. 

At these unique temperatures, we can accurately determine the equilibrium 
liquid and solid solute concentrations for use, for example, as reference points on 
conventional phase diagrams. These solute concentrations can be obtained with an 
accuracy completely independent of the sample shape and size, and even independent of 
the heating rats if the sample is equilibrated at the equi -evaporative temperature. The 
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resultant sample will display a sharp solid-liquid interface. There will be a com- 
pletely uniform liquid composition, and a completely uniform but different composition 
in the solid. No concentration gradients exist in either the liquid or the solid. Any 
portion, large or small, of the liquid or solid is perfectly representative of the same 
phase and resolution of the analytic instrument or technique becomes immaterial. 

We can, therefore, determine these concentrations by Auger, electron probing, 
spectroanalysis, chemical analysis, etc., with almost the same precision. This is 
true even if the alloy contains one or more highly volatile components. The segregation 
coefficient k at this unique temperature can, therefore, also be accurately determined. 

We can also propose the following unique, isothermal solidification experiment. 

A liquid alloy sphere containing 0. 75 mole fraction of tin in aluminum is evaporated at 
slightly below the equi-evaporative temperature, say at 863K. The solvent aluminum 
now evaporates more slowly than solute tin. This condition results in a solute enrich- 
ment in the remaining liquid. Upon reaching 0. 74 mole fraction of tin, the alloy 
solidifies according to the phase diagram. The solidification proceeds from the sur- 
face, but remains at the constant temperature of 863K. The evaporating time can be 
calculated by means of our normal evaporation equations presented in a previous section 
of this report, after the shape and size of the evaporating sphere are specified. By 
controlling the evaporation temperature one can control the rate of solute enrichment 
which, together with mass diffusion constraints, allows us to control the rate of solidi- 
fication. When the evaporation or isothermal solidification temperature is very close 
to the equi-evaporative temperature of 868. 7K, one gets a true picture of solute 
segregation due to solidification rather than the combined effects of solidification and 
unknown amount of evaporation, as is almost always the case in current studies. 

Further, in space, where surface tension is important but gravitational convection is 
negligible, one can start with a homogeneous sphere and freeze it isothermally from the 
spherical surface at a rate so slow that unique cellular and granular structures can be 
obtained. In addition, one can locally control the temperature of the molten sphere by 
means of a heated or cooled sting, or by means of concentrated electron or optical 
energy, so as to initiate single crystal growth in the manner of Walter and Snyder (1971). 
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Section 5. CONSTITUTIONAL MELTING AND SOLIDIFICATION DUE TO EVAPORATION 


f 
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s 



\ 


% 


\ 


I 

» 

.1 


Consider aluminum and tin. A1 melts at 933. 3K and Sn at 505. OK (see Table 1-1). 
This binary system has a equi-evaporative temperature of 868. 7K (Table 1-8), which lies 

between the melting points of the two components. At this unique temperature T (line 

s 

AB in Fig. 1-28), both elements are evaporating at exactly the same rate and no com- 
positional change takes place. At a temperature T > r g , A1 evaporates faster than Sn, 

and (surface) depletion of A1 and enrichment of Sn occur. At a T <T , the reverse is 

s 

true and A1 enrichment and Sn depletion occur. 

According to the Al-Sn phase diagram (Hansen 1958), the binary alloy is liquid 

if it cousins more than 0.26 mole fraction of Sn at T (point D). Let us select a liquid 

s 

composition below line AB and somewhat to the right of point D, e.g. , 0.265 mole 
fraction of Sn at 868. IK (point G). Since at 868. IK A1 is evaporating slower than Sn, the 
alloy composition will constantly increase in Al, and must move isothermally toward the 
left or Al-rich side. Upon reaching the liquidus on the Al-rich side of the phase dia- 
gram, i.e. , point H, the alloy solidifies, first on the surface but progressively toward 
the center of the sample because of mass diffusion. We thus achieve a new isothermal 
solidification of non-eutectic, binary alloy, because of surface constitutional or com- 
positional changes through evaporation. Note that isothermal solidification of alloys 
(not pure metals) can only occur at the eutectic temperatur3. Our constitutional solidi- 
fication experiment, is, however, not limited to the eutectic composition. 

If point G is sufficiently below line AB, we can heat up the molten sample G below 
the line AB to achieve Al enrichment and, hence, solidification on heating. 

Similarly, the solid composition M can be constitutionally melted at a constant 
temperature (along line MN) or at a decreasing temperature (along line MP). 

All four constitutional phase change phenomena, i.e., isothermal melting 
(path 1), melting on cooling (path 2), isothermal solidification (path 3), and solidification 
on heating (path 4), of non-eutectic, binary (ternary, or other multi-component) alloy 
systems appear to have never been reported. Yet these phenomena may well occur in 
space experimentation. In fact, these "anomalies" may greatly confuse the experimenter 
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TEMPERATURE, K 


PATH 1 CONSTITUTIONAL MELTING AT CONSTANT TEMPERATURE 
PATH 2 CONSTITUTIONAL MELTING ON DECREASING TEMPERATURE 
PATH 3 CONSTITUTIONAL SOLIDIFICATION AT CONSTANT TEMPERATURE 
PATH 4 CONSTITUTIONAL SOLIDIFICATION ON INCREASING TEMPERATURE 



Fig. 1-28 Constitutional Melting and Solidification of Al-Sn Alloys 
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if an evaporation analysis is not first performed for the particular alloy system within 
the temperature range of the experiment. 

For Au-Ce alloys. Table 1-8 gives equi-evaporative temperature as T = 131 OK 

s 

(line AB in Fig. 1-29), which is intermediate between the melting points of the two 

components. At a temperature T < T Ce evaporates slower than Au, while at a 

s 

temperature T > T g Ce evaporates faster than Au. 

An interesting point about the Au-Ce system is that, at T = T , the alloy is 

s 

liquid within the following several ranges of Ce atomic concentration: 0. 03-0.16, 
0.30-0.32, 0.37-0.38, and 0. 72 (Hansen 1958). 

An alloy starting at point F and isothermallv heated at 1323K will be initially 
liquid. The alloy will, however, lose more Ce than Au at this temperature, and its 
composition will therefore move along line EF toward the Au-rich side. Solid crystals 
of CeAu will separate out on the surface when the composition reaches point G. Con- 
tinued heating at the temperature of 1323K will cause more intensified surface depletion 
of Ce, i.e. , the alloy composition will further shift toward the left, on the surface, 
through evaporative segregation, but inside the sample through solid and liquid diffusions. 
Eventually, the alloy will have four solid phases representing, in order, the compositions 
CeAu, CeAu , CeAu , and nearly pure Au. The compositions of the phase boundaries can 

u O 

be read off Fig. 1-29 at the intersection points of line EF with the several liquidus lines. 
If the sample is in the form of a freely suspended drop or sphere in space, it will have 
seven, alternately liquid and solid, concentric spherical shells. Such a structure is 
novel and unique, but may be used to determine solid and liquid diffusivities and phase 
boundaries. Since Au has a density of 19.3, compared to 6. 75 of Ce, we may use this 
composite structure to study convection currents or other fluid movements due to gravity 
in spherical shell geometries. 
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Section 6. APPLICATIONS OF EVAPORATION TO SPACE MANUFACTURING 


In summary, surface evaporation of materials has the following important appli- 
cations or effects relative to space manufacturing of materials. 

6.1 MATERIAL LOSSES 

Data have been presented to show, for example, how the solute A1 in a uniformly 

_3 

heated, 6. 0 x 10 -m sphere containing 8 percent by weight of A1 in Ni evaporates. The 
losses of material, mostly A1 at, and the times to reach, various weight fractions and 
surface concentrations are given. The preferential evaporation of A1 makes it difficult to 
control sample weight, size, and composition. The A1 vapor coats the vacuum chamber 
walls and, more seriously, the nearby samples in the Skylab experiments. Worst of all, 
the electron beam source for heating these samples may be contaminated to the point of 
erratic electron emision and loss of control. As a result of these disadvantages, which 
were partly confirmed by actual simulation runs, the Al-Ni alloy was dropped as a 
sample material for the M512 experiment. 

6.2 EVAPORATIVE PURIFICATION 

Evaporative purification would benefit the following processes: floating zone 
melting; single crystal growth; thin film growth; and preparation of materials standards. 

As shown, we presented quantitative data on the purification of nickel alloys with aluminum 
as a solute. We further described, in some detail, the evaporative purification behaviors 
of nickel-iron and nickel-chromium alloys. The purification of beryllium containing nine 
impurities in the ppm range was separately treated. Three reports have been published 
on these data and the necessary predictive equations, and submitted for outside publi- 
cation. A paper on evaporation prepared under the contract was also presented at the 
AACG at Princeton University in July, 1972. 

Not fully appreciated, it appears, is the fact that evaporative purification is 
relative, in terms of speeds and results. It is also a two-edged sword. Thus, while 
with some elements evaporating faster than the solvent we have evaporative purification, 
with other elements evaporating slower than the solvent we have evaporative (surface) 
contamination or enrichment. A freely suspended molten drop in space may, for 
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example, have its surface solute concentration enriched a thousandfold or a millionfold 
by neglected and undetectable trace impurities within a fraction of a second of its 
deployment. From there on, evaporative purification may be completely negligible, 
for any practical length of time. Further, evaporative purification or (surface) con- 
tamination can be so fast as to be unavoidable, or so slow as to be practically non- 
existent. 

6.3 "ANOMALOUS" MELTING AND SOLIDIFICATION 

Because of evaporative segregation, anomalous "constitutional" or evaporative 
melting (on cooling) and solidification (on heating) cau occur, as shown above. We have, 
therefore, pointed out that any solidification study, particularly on nucleation, must be 
preceded by an evaporation compatibility study. We have also suggested unique and defini- 
tive solidification experiments that circumvent or minimize the evaporation effects. 

6.4 DIFFERENTIAL EVAPORATION '.ND JETTING ACTION 

Different rates of evaporation at various surface regions, such as those resulting 
from nozzle detachment, temperature gradient, cr surface co.iidmination, give rise to 
unbalanced forces and ’.r.umenta that may produce erratic or unwanted accelerations. 

Very large "equi-valent gravities" may be produced. Another effect of such jetting 
actions may be the formation of powerful convection currents with entirely new charac- 
teristics that can not be directly correlated with those of other convection currents. 

6.5 SURFACE CCX) LING 

Evaporation requires the absorption of the heat of evaporation. Surface cooling 
therefore results. For nickel at its melting point of 1726K, the evaporative surface 
cooling component is about four orders of magnitude smaller than the radiative cooling 
component and is, therefore, probably negligible. For low-melting materials, particu- 
larly those having high heat of evaporation, the evaporative cooling component may be 
very important. Work by Dr. The mas T. Kassal of the Grumman Research Department 
with water or brine streams jetting into vacuum demonstrates the impcr.ance of surface 
cooling in ice formation. 
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6.6 MATERIALS STANDARDS 

Properly controlled evaporation, particularly in space, allows us to obtain 
material purity or evaporation standards, thermal properties, and even basic thermo- 
dynamic constants. 


s 



I 
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PART II 


FREEZING SEGREGATION 


This section on freezing segregation represents, in part, a review paper the 
principal investigator has been asked to prepare for inclusion in "Advances in Materials 
Science. " While this work was not done under NASA support, it is included for the sake 
of completeness and because of its relevancy. 

Section 1. INTRODUCTION 

In this section, we will discuss the normal freezing process, and its application to 
steady-state solidification in the Ni-Sn system. 

In the first part of this section, the normal freezing process is described in some 
detail. The assumptions, its industrial application and importance, the normal freezing 
differential equation and its solution for the various cases are given. The meaning of the 
segregation coefficient k is then discussed. The best value of k and the errors involved 
in assuming constant k for various cases are studied. Finally, numerical computation 
of the normal freezing process is exemplified. 

The formation of spirals or screw dislocations and other defects has always received 
great attention because these defects make all materials very weak compared to their 
theoretical potentials. During this contract, new spiral growth observations on barium 
st itium niobates grown under controlled conditions were made. These observations 
cannot be explained by the classical core-first growth theory (or BCF model) A new 
mechanism was therefore developed to explain them. We believe that knowledge of the 
formation mechanism of these and other defects will enable us to control or eliminate 
them in space manufacturing and other operations. Details of this new mechanism is 
given in Appendix F. 
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Steady-state solidification is related to the important process of single crystal 
growth. In particular, the critical temperature gradient must be accurately determined 
in order to achieve interfacial planarity and single crystal growth. We used a non- 
constant k approach to study the steady-state growth of Ni-Sn alloys. Our results show 
that assuming constant k may introduce intolerable errors in predicting solidification 
behavior in general and single crystal growth in particular. 

Section 2. THE NORMAL FREEZING PROCESS 

In normal freezing, we assume that the liquid is completely mixed, that no atomic 
diffusion occurs in the solid, and that the compositions of the freezing solid and liquid, 
respectively, follow the solidus and liquidus curves of the constitutional diagram. 

Since the liquid is completely mixed, there is, effectively, infinite diffusion in the 
liquid. The liquid is, therefore, always homogeneous. Solute pile-up in front of the 
liquid-solid interface is not possible and solute concentration at this interface is, therefore, 
the minimum possible. The equilibrium solute concentration in the freezing solid is also 
minimal. Thus, the interface acts as if it were an efficient atomic filter, allowing pure 
solvent atoms to freeze out together with the minimum amount of solute atoms while 
keeping the maximum amount of solute atoms in the remaining liquid (for the case where 
the segregation coefficient is less than 1). Thus, normal freezing is the most efficient 
purifying or segregating process, as far as separating solute and solvent atoms in the 
liquid (according to the constitutional diagram) is concerned. 

In contrast to i ’finite liquid diffusion, the solid diffusion is zero during normal 
freezing. Thus, all the (solvent and solute) atoms stay exactly where they first freeze 
and remain completely and forever immobilized. Justifications of these assumptions, 
as well as additional comments and implications, will be more fully described. 

Normal freezing is a useful and important ideal limiting case, as is infinitely rapid 
quench or splat cooling. In the rapid freezing process, time is not allowed for any solid 
and liquid diffusion, or even equilibration, according to the constitutional diagram. The 
effective solid and liquid diffusion constants are thus zero. The solute and solvent atoms 
are instantly frozen, in situ, from the melt; and the constitutional diagram is completely 




irrelevant. There is also no solvent purification or solute segregation. In practice, when 

the characteristics distance d , defined as the liquid diffusion constant D, divided by the 

c 1 

growth rate v, falls appreciably under the resolution limit of the detecting instrument, the 

frozen material appears to be homogeneous and the freezing process is practically identical 

-4 2 , 

to that of an infinitely rapid quench. Since is generally about 10 cm , sec, a growth 

rate greater than 10 cm/sec from a severe quench yields a characteristic distance of 

d = D /Kv = 1 pm. Such a frozen material would thus appear homogeneous under the 
cl 

electron microprobe, which has a resolution of about 1 pm. 

Another useful and limiting freezing process is equilibrium freezing according to 
the level rule. Here, both the liquid and solid are homogeneous, and the liquid-solid 
interface is perfectly permeable to the solute and solvent atoms. At all times, the liquid 
and solid are in equilibrium according to the constitutional diagram, and unique, 
"equilibrium segregation" results obtained. 

Any actual freezing results must fall within these maximum segregation, no 
segregation, and equilibrium segregation limits. For some alloys systems, the limits 
between maximum and equilibrium segregation results may be only a few percent apart 
(Li, 1966). Further, the relative proximity of practical freezing results to either of these 
limits can often be at least qualitatively predicted bv a study of the freezing conditions. 

Among the three methods, normal freezing is generally the most important since 
it gives the maximum solvent purification or solute segregation, the purest possible 
single crystals, and the most defect-free devices or structures since many defects are 
often associated with segregated impurities. Also, it is particularly valuable in the pre- 
paration of many material specifications which call for precisely these limiting case con- 
ditions and require the maximum initial solute content, the maximum defect size of a 
given type, the maximum grain boundary thickness for grains of a given size, and the like. 

Normal freezing often does not represent the true freezing conditions. This is 
particularly true in dendritic growth when the growth velocity is high and solid diffusion is 
negligible (normal freezing assumption 2) but liquid diffusion is not < omplete (normal 
freezing assumption 1). With single crystal growth, the growth speed is often orders of 
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magnitude lower, solid diffusion is still negligible, and the liquid diffusion conditions are 
more nearly normal. In addition, if the melt is mechanically or induction stirred, the 
melt may be indeed nearly homogeneous and normal freezing conditions are realized. 

Even a lack of normality does not detract from the usefulness of the normal freezing 
methodology. This is partly because actual freezing conditions are more or less normal. 
The departure from normality can be predicted and actually evaluated, as will be shown. 
Further, normal freezing gives limiting results which are useful for writing materials 
specifier. ions. As will be shown, normal freezing conditions give, at each step during 
the freezing process, the maximum crystal purity possible (for the usual case where k < 1). 
Thus, normal freezing calculations allow us to prescribe the maximum initial impurity 
concentration to achieve a given minimum degree of purity in the final product, a given 
thickness of grain or cell boundary thickness, or even given mechanical properties such as 
creep and ductile-to-brittle transition temperature. In addition, abnormal conditions 
such as limited liquid diffusion and loss of materials through surface evaporation or 
chemical reaction can be accounted for, as will be shown. 

2. 1 NORMAL FREEZING AS USED IN THE SEMICONDUCTOR INDUSTRY 

Normal freezing techniques has been successfully used in the semiconductor 
industry and are indispensible in such areas as material purification, single crystal 
growth, dopant control, diode formation, and liquid epitaxy. Here, the conditions dis- 
cussed below are most favorable for the use of these techniques. 

The most common germanium and silicon materials are very pure. They generally 
contain only ppm or ppb of impurities and are, therefore, orders of magnitude purer than 
most other pure materials. We are thus dealing with extremely dilute, simple or ideal 
solutions. 

There is a critical need for control, at microscopic distances, of absolute depant 
IS 15 

levels to ppm (~ 10 /cc) or ppb (~ 10 /cc). Even the dopant concentration gradients at 
these extremely low concentration levels and small distances are critical. In fact, even 
ppb of wanted or unwanted impurities can radically change the characteristics of the 
resultant semiconductor devices. 
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Only a portion, often a small portion, of the solidified material is used; and the 
solidification range of interest is often exceedingly narrow. In the growth of single crystal 
of pure germanium containing 100 ppm of antimony in the melt, for example, the growth 
starts at 1209. 10K or 935.98°C (k to _. , At _ = 0.276), and ends at 1208. 12K or 935°C 
(k - 0.276), and over 98.0% of the melt by weight (representing more than the 

material of interest) is grown (Li, 1966). The liquidus and solidus curves of interest can 
be represented by two tiny straight lines extending over solute concentration ranges of only 
4,660 and 130 ppm, respectively, over a temperature interval of 1°K. The segregation is 
constant at 0. 276 to within the third decimal place. Again, in the growth of a germanium 
alloy diode with a 50-micron indium sphere, only the first 0. 1-micron linear distance 
(representing less than 0. 5% by weight solidified) is critical in the formation of the pn 
junction (Li, 1966). Here, the useful solidified portion is extremely small, as is the 
solidification temperature range of interest. The segregation constant k is, therefore, 
again constant. Pfann's constant-k, normal freezing equation (Pfann, 1952) thus finds its 
ideal application. 

Most germanium and silicon alloys have extremely low solute, solubility limits in 
the solid. According to the constitutional diagrams (Hansen, 1958), germanium alloys 
containing Ag, Al, Au, Cd, Co, Cu, Fe, In, Mg, Mn, Ni, Pb, Sn, Te, Tl, and Zn, and 
silicon alloys containing Ag, Al, As, Au, Ca, Ce, Co, Cr, Cu, Fe, In, Mg, Mn, Nb, 

Ni, Pd, Pt, Sb, Sn, Ti, and Zr, have invisible solidus lines. Only germanium alloys 
containing As, Bi, Ga, and Sb, and silicon alloys containing U, have some visible solidus 
lines that are very close to the vertical temperature axes. 

The segregation coefficients for the above germanium and silicon alloys are very 
small, being close to zero. Also, the solute material frozen in the solid is negligible, 
and does not appreciably affect the fraction solidified. It is therefore possible to con- 
sider multiple, non-interacting dopants in the same crystal for the determination of the 
pn junction depth. Extending this approach indiscriminately would be misleading or 
meaningless. Unfortunately, this has been done. 
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2. 2 EXAMPLES OF CONSTITUTIONAL DIAGRAMS 


Figure 2-1 shows several common types of constitutional diagrams that we deal 
with in solidification studies. Figure 2-l(a> represents a simple eutectic system typified 
by the Pb-Sn alloys. Here, solid or Pb crystals are shown to have different concen- 
trations of Sn, reaching their maximum Sn concentration of c gu at the eutectic temperature 
of 456. IK (183°C) as described by the left, or Pb-rich, solidus line. Similarly, solid 0 
or Sn crystals also have different concentrations of Pb, also reaching their maximum Sn 
concentration of c at the eutectic temperature, as represented by the right or Sn-rich 
solidus line. At each melt temperature, the « or 0 crystals are in equilibrium with only 
one single melt composition, as determined by the appropriate liquidus line. 

Figure 2-l(a) also shows that the four liquidus and solidus lines are far from being 
straight lines. The liquidus and solidus of the Sn-rich alloys may be approximately 
•presented by quadratic functions of the melt temperature. The liquidus and solidus of 
the Pb-rich alloys, on the other hand, would probably be more adequately fitted with 
cubic equations of the melt temperature. The exact degree of such fits should, however, 
be determined objectively by statistical tests (Anderson and Bancroft 1952). 

Figure 2-l(b) shows the complete constitutional diagram for the Au-Ag system. In 
this system, the Au and Ag are completely miscible, both in the solid and liquid phases. 
Hence, there is only a single liquidus and a single solidus for the entire system, instead 
of the two liquidus and solidus lines shown in Fig. 2-l(a). It is immaterial whether we 
consider Au or Ag as the solvent. It is even arbitrary or artificial to consider the major 
constituent as the solvent and the minor constituent as the solute, since such a designation 
depends on whether weight percentages or atomic percentages are used as the classifying 
criterion. Thus, an a Joy containing between 50 to 67 weight percent of Au would be 
considered as an alloy of Ag in Au if weight percentages are used and, yet, should be 
considered as an alloy of Au in Ag if atomic percentages are used as tne classifying cri- 
terion. In many other systems, such as W-C, Ta-C, and W-B, this disparity is even 
wider and the designation of one component as the solvent or solute is even less meaningful. 
It would, therefore, be better to identify an alloy of, e. g. , 10% of Ag in Au simply as 
such; or, interchangeably, as an i lloy of 90% of Au in Ag. 
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Fig. 2-1 Some Types of Constitutional Diagrams 
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A good system of solidification methodology, in particular the normal freezir g 
equation, should be a self-consistent one. It will be shown using this Au-Ag diagram 
and several other considerations that the constant-k method is not such a system. 

Figure 2-l(c) shows the Cu-Au constitutional diagram. This system is characterized 
by having very large solid solubilities at either terminal region. 

Figure 2-1 (d) shows the Au-Zn constitutional diagram, revealing the presence of a 

number of intermetallic compounds, i.e., AuZn, AuZn , and AuZn . Furt .er, at the 

2 * o 

compositions of the first two compounds, the liquidus and solidus reach their peak tem- 
peratures. In many such systems involving intermetallic compound formations, the liquidus 
lines, and even the solidus lines, have horizontal tangents (k = ®/® = indefinite) at their 
peak temperatures at the most important compositions. This condition has important im- 
plications in the determination of the segregation coefficients k, as will be described. 

Figure 2-l(e) shows a portion of the Fe-C diagram. An interes :ng feature is that 

a peritectic reaction occurs at a T f of about 1765K (or 1492°C). At this temperature, the 

solidification process is discontinuous. A small portion of the Fe-C system has two 

distinct liquidus and two distinct solidus lines, one set above and one set below T . The 

slopes of the two liquidus lines are discontinuous at T . Not only are the slopes of the 

two solidus lines discontinuous, but even the values of the solute concentrations in the 

solid, i.e. , c , are discontinuous. Specifically, c changes abruptly from 0. 10% at 
s s 

slightly above T^ to 0. 18% at slightly below T . The study of solidification of Fe-C 

alloys in this portion of the diagram thus requires dividing the solidification process into 

two regions, one above and one below T . In each region, a unique set of liquidus and 

r 

solidus lines should be used. 

2.3 THE NORMAL FREEZING DIFFERENTIAL EQUATION 

Let c, and c be functions of a modified melt temperature T = T - T , where T 
Is erne 

is the fixed melting point of the pure solvent and T is the variable melt temperature 

m 

(Fig. 2-2). Thus: 
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Fig. 2-2 The Normal Freezing Process 
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Cl = L(T) 


c = S(T) 
s 

where L(T) and S(T) are functions of the modified melt temperature, for the solidus and 
liquidus curves of the constitutional diagram. 

At the b< ginning of normal freezing, time t = t = 0, the fraction solidified p = P Q = 0, 

c, = c, = c , T = T , and k = k = S(T )/'L(T ), where T can be determined from 
1 lo o o i ' o ' o o 

c = L(T ). The frozen length isx = x =pH=0, where H is the total melt or crystal 

o o o o 

length. 

At t = t, the fraction solidified, frozen length, and modified melt temperature are, 
respectively, p, x = ph, and T. At the end of the next short time interval A t from t, these 

values will respectively increase to p + Ap, x + A x, and T + A T. 

Balancing the amount of solute Q in the remaining melt at time t according to a 
general procedure previously described (see, e. g. , Pfann 1952) gives 

AQ = - Cg Ap - A j(l-p)c 1 J -f x (c^ A t, T, . . . ) (24) 

where represents tne loss of solute from causes such as surface evaporation and 
chemical reaction. This loss may depend on c^, A t, T, melt volume, or free surface. 

In many cases, f^ = 0. If we then take the limit as A t, A p, A x, and A T all 
approach zjro, we obtain the following basic differential equation for normal freezing: 

dp/(l-p) = -d ln(l-p) = | L’ (T)/ [L(T)-S(T) ] JdT (25) 

where L'(T) is the temperature derivative of the liquidus curve L(T). 

In the above derivation, e and c may be in terms of weight fractions. Q then 

1 s 

represents the weight of solute in the liquid melt for each one gram of original melt. If, 
however, c and c are in terms of atomic fractions, Q then represents the number of 

L 8 

moles of solute atoms for each mole of combined solvent and solute atoms in the original 
melt. 
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2.3.1 Constant-k Normal Freezing Equation 

In a few constitutional diagrams, the liquidus and solidus curves can be approximated 
by straight lines. When we are dealing with very narrow ranges of freezing temperature, 
solute concentration, or fraction solidified, such as in the previous examples of growing 
pure germanium and silicon crystals and of forming germanium alloy diodes, the portions 
of the generally curved liquidus and solidus lines of interest may be so small that they 
can be represented by two straight lines. We then have: 

c = aT + e 1 

c = fT + e 
s s 

where a and f are the slopes of .he liquidus and solidus lines, and e and e are the 

1 s 

intercepts of these lines, when extended, on the vertical temperature axis (see Fig. 2-2). 
The normal freezing differential equation, Eq (25), then reduces to 


dp _ a 

(1-p) (a-f)Tf(e-e) 
1 s 


dT 


This equation is easily integrated. When proper initial conditions are included, we obtain 
the following general, constant-k, normal freezing equation: 

,ta/(f-a) 

(26) 


’(a-f)T + (e 1 * 
~ P= (a-"f)T o + (e x 


-v 

•Vj 


When pure germanium crystal or diode materials are dealt with, only the terminal 

regions of the constitutional diagrams are used. The two hort liquidu and solidus li-ies 

then usually both pass through the melting point of pure germanium, and both e, and e 

Is 

are zero. Equation (26) then reduces to the famllar, constant-k normal freezing 
equation: 


T/T = aT/aT = c / c = c /kc = (I-p) < f “ a)/a = (* .p^' 1 
o oloso ' v 


( 27 ) 


! 
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This equation is so important in materials processing, particularly in the semi- 
conductor industry, that it has been repeatedly derived (Gulliver 1922, Scheuer 1931, 
Hayes 1939, Luneau 1942, Mathier 1942, Schell 1942, McFee 1947, and Pfann 1952). 

An analogous equation, known as Rayleigh's equation, has also been derived and used 
for "simple" distillations in chemical engineering (Rayleigh 1902). 

2. 3. 2 Accuracy of the Constant-k Normal Freezing Equation 

According to Christian (1965), Pfann's equation applied to reasonably dilute 
solutions, is a useful approximation until some 80$ of the liquid has solidified, that is, 
until p = 0. 80. 

Dowd and Ro (e (1953) studied the distribution of radioactive indium in germanium, 

and concluded that the indium concentration obeyed Pfann's equation. On the other hand, 

Goruw and Jentoft (1967), in evaluating the effectiveness of normal freezing, column 

crystallization, and zone melting of n-C, H in n-C H . found that the segregation 

14 30 16 34 

coefficient is a linear function of the solute concentration in the melt. 

Kurgistsev (1967) analyzed .he ptoblem of determining k from experimental dat". 
and showed that the methods which reduce this problem to that of finding a solution of the 
equation on the assumption that k is constant lead to erroneous results. 

Brody and Flemings (1966) analyzed their data on the Al-Cu system using their 
semi-theoretical, constant-k normal freezing equation. They emphatically concluded 
that the constant-k equation gives "accurate" results, although they had to multiply their 
experimental data by an unfounded factor of three to achieve reasonable fit . Johnston 
and Tiller (1962) considered the problem of the segregation coefficient k, oj effective 
segregation coefficient k, varying as a function of the interface position (or p), and 
determined some of the functions required to produce particular solute distributions in 
the solid. Wilcox (1964) also considered k to vary with c , and set up the necessary 
normal freezing and zone melting equations for numerical integration. 
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Pfann (1958) stated that the linear normal freezing equation, Eq 27, does not 
describe the complete freezing behavior for the whole length of the ingot. The complete 
concentration curve, however, can be calculated by Hamming's method (Pfann 1958) 
by taking the dependence of k on the concentration into account. 

2. 3. 3 Non-constant-k, Normal Freezing Equations 

As pointed out above, Eq (27) has been very successfully used to determine gross 
concentration profiles in the useful crystal, and other solidification behavior, of pure 
germanium and silicon materials. This is true even when the liquidus and solidus lines 
are distinctly curved (Li, 1966) due to the unique conditions previously described. In 
other normal freezing studies related to general metallurgy, however, the entire 
solidification range must be considered. Such studies include the investigation of the 
segregation of trace impurities at grain, subgrain, cell, or subcell boundaries; the 
formation of freezing-in or induced point, line, surface, and massive defects in elec- 
tronic or structural materials (Li, 1967bl; the selection of homogenization, oxidation- 
precipitation, diffusion-reaction, or other heat-treatments from the calculated dis- 
tribution of alloying elements in polycrystalline welded or soldered joints; and*the magni- 
tude and character of voids or vacancies, and of » 3sidual stresses or strains, in castings, 
caused by variable solid and liquid material densities. 

2.3.4 Quadratic Liquidus and So l idus Lines 

When the solidus and liquidus can be represented by second-degree functions of the 

melt temperature T = T - T, i. e. , 
m e 

c « aT + bT 2 = a(T -T ) +■ b(T -T ) 2 
1 e m ' e m 

c = iT + gT 2 * f*T -T ) + g(T -T ) 2 
s - ' e m **' e n 

Eq (l) then becomes: 


-d ln(i-p) 


(a-2’oT) c’T 
(a-f)T + (b-g)T 2 
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Hence, 


l-p = (T/T 


o’ (a- 


-n + (b-£ 


f) + (b-g)T 


where Q = a/(f-a) and Q = (ab^ag-2bf)/(f-a) (b-g), as given by Li (1967a). 

Li further showed (1967a) that Eq (28) reduces to Eq (27) only under the restricted 
conditions of very dilute melts, nearly straight liquidus and solidus lines, and very 
small solidification ranges. Hence, when one or more of these conditions are not met, 
the use of Eq 27 in normal freezing studies requires caution. 

2.3.5 Cubic Liquidus and Solidus Lines 

An examination of binary constitutional diagrams (with k < 1) reveals that, of 171 
binary diagrams given in the ASM Metals Handbook (1948), 36 (or 217 C ) have distinctly 
cubic liquidus and/or solidus lines. It is, therefore, desirable to derive the normal 
freezing equation for cases with cubic liquidus and solidus lines, i. e. , when: 

c = L(T) = aT + bT 2 + dT 3 

c = S(T) = fT + gT 2 + hT° 
s 

For this case, Eq (25) reduces to (Li, 1968): 


-dln(l-p) 


where 


t aj2bT Od'fV dT 

(a-f)T + (b-g)T 2 + (d-h)T 3 T T+U T+V 


u = jb-g^b-g) 2 - 4(a-f) (d-h)j / 2(d-h) 
v = |b-g - ^ (b-g) 2 -4(a-f) (d-h)j /2(d-h) 


Q = -a/uv(d-h) 


Q 2 = - juv(3du-2b) + avj /us'(d-h) (u-v) 
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and 






Qg = - |uv(2b-3dv) + auj/uv(d-h) (u-v) 



( 29 ) 


2.3.6 Quartic and Higher Degree Liquidus and Solidus Lines 

Here, the liquidus and solidus lines, L(T) and S(T) respectively, can be repre- 
sented by two polynomials in T such that L(T)-S(T) is an nth degree function, and the 
degree of L’(T) is less than n. Then, Eq (25) yields (Li, 1968): 

-dln(l-p) = [l- (T)/ [l(T)-S(T)1 1 dT = 2 Q/(T-T ) dT = £ Q d ln(T-T ) 

{ L J ] i=l i=l 1 


from which 


n 

1 - p = n 

i=l 


j^(T-T.)/(T o -T i )] 



(30) 


the T^’s are the roots of the polynomial L(T)-S(T), and the Q.'s are the numerators of 
the various partial fractions for the function L'(T)/ [(L(T)-S(T)] . 

2.3.7 Ideal Binary Systems 

In an ideal solution having an ideal, dilute solid solution in equilibrium with a 
simple melt solution (Panish and Ilegems 1970), the liquidus curve (in mole fraction) is 
given, according to Prigogine and Defay (1954), by: 

c 1 = 1 - exp Ju/T o -l/T) A H/rJ = 1 - exp (a+b/T) = 1 - U 

where AH is the molar heat of fusion, T q is the absolute melting temperature of the pure 
solvent, T is the absolute melt temperature, b = - AH/R, a = -b/T, and U = exp (a+b/T). 

For an ideal solution in equilibrium with a dilute ideal solution, the segregation 
coefficient for the dilute component can also be shown to be an exponential function of 
the reciprocal absolute melt temperature, both theoretically (Hildebrand and Scott 1950; 


1 ! 


1 

i 
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Thurmond and Struthers 1953) and experimentally (Hall 1957; Orioni and Hall 1958). Thus 

k = exp (P + g'/T) 


and the equation of the ideal solidus curve is 

c = exp (a + b/T) exp (f' + g'/T) = exp (f + g/T) = V, 
s 

where f', g', f = a + f', and g = b + g' arc constants, and V = exp (f + g/T). 

The basic normal freezing, differential equation, Eq (25), then becomes: 


j, ,, . exp tato/iM-p/i )ai 

n P exp (a+b/T) - exp (f+g/T 


Hence (Li, 1971), 


exp (a+b/T) /exp f f+g/T 1 - 1~ 
exp (a+b/T^J/exp (f+g/T Q ) - 1 


b/(g-b) 


, D = [£XElg±/lL4 gx plftg/TI „-_l _/ U/V-l \ 

" p exp (a+b/T )/exp (f+g/T ) - 1 U /V -1 j ' ; 

_ 0 o _j \ o o / 

where T q = b/£ln(l-c o )-aJ, U q = exp (a+b/T Q ), = exp (f+g/T^), and = b/(g-b). 

2. 3. 8 Ideal Ternary Systems of the Pseudobinary Type 

Many ternary liquid solutions are ideal in that the A, B, and C atoms in the 
solution are solidified into ideal solid solutions of AB-AC with the A atoms on one 
sublattice and B or C atoms on the other. Such systems can be treated as pseudobinary 
systems. 

In such systems, the mole fraction for the AB compound in the solid is (Panish 
and Ilegems 1970): 

AB 

c = (1-V)/(U-V) 

S 

where U = exp (a+b/T), V = exp (f+g/T), b = -AH^/R, a = -b/T^, g = -AH^/R, 
f = “g/T^£» AH^g is the heat of formation of compound AB at its absolute melting 
point T^g, and AH^, is the heat of formation of compound AC at its absolute melting 
point T . The mole fraction of B atoms in the liquid solution is (Panish and Ilegems 1970). 
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B _ _ AB, 
c = U c /2 
1 s 

A B C 

By introducing chemical molar concentrations m , m , and m , we can balance 

the amount of B or C in the remaining liquid. Thus, by accounting for the amount of B 

in the liquid, Q , Eq (25) becomes (Li, 1972): 

B 


-dln(l-p) = 


;[b(i-v) - K (i-al d(l/T) 113, -^_] d(1 / T , 

L-U) (1-VO (1-V/U) d(1/ ' Ll-U 1-V l-V/UJ dU/1) 


Hence, 


Q 1 Q 2 
-i -ii 1 

, l-U 1-V 1-U/V 

^ “ p ~ -1 _1 

l-U 1-V l-U /V 

o o o o 

•4 p 


where 


Q 1 = q/(q-l). 


Q 2 = l/(l-s). 


Q 3 = r/(l-r). 


q = exp (f-ga/b), 


= exp (ag-bf)/ (g-b)j 
s = exp (bf-ag)/ gj 


T o , the initial freezing temperature, can be determined from the pseudoliquidus 

AB AB 

line by substituting the initial melt composition c for c 

l,o 1 

Complete details of the derivation of Eq (32) and its application to the theoretical 
computation of the normal freezing behavior of GaAlAs system is given in Appendix D. 
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2.3.9 Summary of Normal Freezing Equations 

The closed-form, normal freezing solutions described above, but expressed in 
slightly different torms, are given as follows: 


• Linear 


_1_ 

k-1 


8 „ >1 i 

— = (1 - p) or 1 - p = 1 
KC o 


-J?l -(f) 


k-1 


• Quadratic 


Q, Q, 


l - p = 


/_T\ / a + bT \ 
l T oj l a + bT oj 


• Cubic 


1 - p 

• General Polynomial 


[ T - 'if 1 r T - t 2 tV - t 3 f 

l T o- T J L T o- T 2 J L T o - T 3 J 


1 -p = 


» r T ~ T i 

i-1 T “ T - 
i-l o i 


Q. 

- i 


Simple Binary Solutions 


1 - P = 


r U/V-l I 

lVV‘J 


• Pseudobinary 


1 - p = 


l l-U 

1-1 


-1 




-11 


p7(M 


( 27 ) 


(28) 


(29) 


(30) 


(31) 


(32) 
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2.4 APPLICATIONS OF NORMAL FREEZING EQUATIONS TO SOLIDIFICATION 

Application of the various normal freezing equations to solidification problems 
will now be described, with emphasis placed on the constant-k normal freezing equation 
(or Pfann's Equation) which is so important in the semiconductor industry. Specifically 
considered are several methods of fitting straight lines to liquidus and solidus curves; 
the meaning of the constant segregation coefficient; the best value of a constant segregation 
coefficient for a given solidification problem; and errors in the use of a constant k when 
the liquidus and solidus are curved. Next, we will show that the last portion of a freezing 
melt must be a eutectic having a fixed composition, i. e. , the eutectic composition. 

Finally, we give two useful numerical computation methods, i.e. , the Runge Kutta Method 
and The Finite Difference Method, to treat cases not covered by the various normal 
freezing equations given above. 

2. 4. 1 Fitting Straight Lines to Liquidus and Solidus Curves 

Most liquidus and solidus lines are more or less curved. There are many ways a 
straight line can be fitted to a curve. The most obvious method is to fit a straight line 
through the discrete data points (obtained experimentally or from the constitutional 
diagram) to minimize the sum of square of the residual errors (line 1, Fig. 2-3a). This 
is the standard, least-square fitting technique (Anderson and Bancroft, 1952). For most 
solidification work, however, this technique is improper. In the case of pure Ge, Si, or 
other materials, particularly during single crystal growth, we are mainly interested in 
the terminsl region (around the point E) of the constitutional diagram. Line 1, however, 
often gives the highest errors in these regions. Further, these fitted lines often cannot 
be used in the various freezing equations, failing even to meet at point E for zero solute 
concentrations. These lines are thus physically unsatisfactory and often meaningless. 
Moreover, line 1 does not give the correct eutectic compositions. 

A second fitting method is to pass both straight lines through the melting point of 
the pure solvent (i.e. , point E), while still minimizing the sum of squares of the residual 
errors (line 2, Fig. 2-3b). Unfortunately, most liquidus and solidus lines have curvatures 
around point E. If these lines are convex upward, line 2 will consistently underestimate 
the melt temperature for a given solute concentration, and vice versa. Further, the 

2-19 



i 

& 

i 

*>>• 


r- 

i* 

r r 

*• 



Fig. 2-3 Methods of Fitl 
and Their Assi 


E 



(b) 




tight Lines to Curved Liquidus and Solidus 
Segregation Coefficients 


T 



all-important segregation coefficient k is markedly off, particularly when one of the 
liquidus and solidus lines is convex upward while the other is convex downward. The 
estimated eutectic compositions may be erroneous and useless in eutectic solidification 
work or in estimating the eutectic thickness at the grain boundaries. 

Another method, particularly useful for dilute solutions or pure materials, is to 

pass the straight lines through the point E and tangent to the liquidus or solidus curve 

(Fig. 2-3c). These fitted straight lines, however, generally give the greatest errors 

in the estimated eutectic compositions. Further, for many liquidus curves, particularly 

those for intermetallic compounds, the liquidus is often nearly horizontal at the terminal 

region. The segregation coefficient k=k must then be zero or nearly so. The 

constant-k, normal freezing equation, Eq (27), then reduces to the meaningless form 

c = 0. 
s 

In Fig. 2-3d, the two straight lines pass respectively through the c gu and c 
points while minimizing the sum of squares of the residual errors. While this set of 
lines is sometimes sufficient for eutectic solidification work, it is evidently inadequate 
for work on pure materials or dilute solutions. 

Figure 2-3e shows a set of fitted, straight liquidus and solidus lines that intersect 
the same c gu and c points, as in Fig. 2-3d, but are tangent to the respective liquidus 
or solidus curves. For most eutectic solidification with initial melt compositions near 
the eutectics, this set of straight lines is superior to that of Fig. 2-3d. 

Figure 2— 3f shows a set of fitted straight lines that pass through the points E and 

c ore . The segregation coefficient k = k is seen to be greatly overestimated 
su lu tu 

here. 

There are other methods of fitting straight lines. We could, for example, assign 
different weights to different temperature or concentration ranges. These methods are, 
however, too specialized to be discussed here. 
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2.4.2 The Meaning of Constant Segregation Coefficient 


fr; 

£ 


% 
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With a few exceptions, such as Tiller (1962), Wilcox (1964), and Rubenstein (1971), 

practically all researchers in solidification assume a constant segregation coefficient k 

and, hence, perfectly straight liquidus and solidus lines on the constitutional diagrams. 

This coefficient k is generally defined as the ratio of solute concentration in the just 

frozen solid to that in the remaining melt, i.e. , k = c /c . The values of c and c may 

si si' 

be obtained experimentally or from the constitutional diagram. When dealing with dilute 

solutions or pure materials such as Ge, Si, W, or Ti, we find that both c and c 

s 1 

approach zero and, hence, k = 0/0 which is undefined. In this case, we define k as the 
ratio of the slope of the solidus at point E (Fig. 2-3) to the slope of the liquidus at the 
same point. Thus, k becomes k Q = k £T of Fig. 2-3c. If c g and c are fitted by two 

” i n i 

n-degree polynomials of the melt temperature, i.e. , c - Z f T and c = Z a T , 

8 i=l 1 i=l 1 

then k = k is the ratio of the first fitted coefficient in the solidus equation to that 
o El 

in the iiquidus equation. Thus, k = f /a, . 

oil 

It has been shown that none of the several methods of fitting straight lines to 
curved liquidus and solidus generally give accurate results in solidification work. 
Basically, it is impossible to accurately fit straight lines to these curves. The six 
segregation coefficients, each given for a particular type of fitting in Fig. 2-3, usually 
differ widely in value. Each is useful only in its limited range of application. Given a 
specified solidification problem, it is often difficult to decide which constant k should 
be chosen. Even with a properly chosen k, when a specific initial melt composition is 
given, the beginning of solidification and, hence, the solidification range, cannot 
generally be correctly ascertained. The upper ends of the fitted straight lines in Fig. 
2-3a, 2-3d, and 2-3e do not even intersect at the point E, as they should. Ridiculous, 
negative solute concentrations and, hence, negative segregation coefficients may therefore 
appear. 


\ 
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Each of the several types of fitted straight lines in Fig. 2-3 gives a unique type 

segregation constant valid only in certain solidification applications. None is generally 

useful, each being often greatly different from the others, This is easily seen, e.g. 

in the Ge-Sb system. The computed segregation constant from the constitutional diagram, 

i.e., k = c /c , varies from 0.0276 at the melting point of germanium (i.e., 1209K 
CD si q 

or 936°C) to 0.0390 at the eutectic temperature of 863K or 590 C. This represents an 

increase of 41.3% (see Li 1966a). In a very dilute solution, e.g., in most semiconductor 

work including the growth of nondegenerate germanium crystals, a segregation constant 

k = k =■ 0.0276 should be used for constant-k freezing equations, such as Eq (27). In 
o 

a Sb-rich melt near the eutectic composition, however, a k = k = 0. 0390 should be 

L 

chosen. 


In the case of dilute Ge-Sb solution, if we use k^ instead of the correct k Q , then 

Ak ■= k - k = + 0. 0114, and the percentage error in k is Ak/k = + 0.413. The percen- 
U o ° A k Ak 

tage errors in c computed according to Eq (27), i.e., Ac/c = (1 + ) (1-p) -1, 

are 41.1%, 40.9%, 40.7%, 40.5%, 40.2%, 39.8%, 39.4%, 38.7%, 37.6%, and 36.6% 

for p = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 0. 95, respectively. The 

percentage errors in p computed from Eq (27), i.e. , 


, Ak/k-Ak , 
Ap/p - P "I 


are -2.69%, -1.89%, -1.42%, -1.08%, -0.8i9%, -0.604%, -0.422%, -0.264%, -0.125%, 
and -0.00608%, respe. lively, for the same values of p. 


If, in the case of melts of near-eutectic compositions, we use k Q instead of the 

correct k ,, then Ak = k - k T = -0. 0114 and the percentage error in k is now Ak/lc. 

U o U u 

= -0.292. The percentage errorsin the computed c's are then -29.1%, -29.1%, -28.9%, 
-28.8%, -28.7%, -28.5%, -28.3%, -27.9%, -27.3%, and 26. 8%, respectively, for the 
same ten values of p. The corresponding percentage errors in the computed p's are 
now 2.74%, 1.90%, 1.42%, 1.08%, 0.816%, 0.601%, 0.419%, 0.262%, -.124%, and 
0. 00602%, respectively. 


Yet, the Ge-Sb system should be considered as a fairly well-behaviored system. 
The possible errors in k and, hence, computed c’s and p’s using any constant-k freezing 
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equations, are even much worse with many other systems. Consider, for example, the 

Au-Ag system having a constitutional diagram of the type shown in Fig. 2-lb. Any 

Au-Ag alloy can, therefore, be treated either as an Au alloy containing Ag as the 

solute, or as an Ag alloy containing Au. Considering *. .~i as an Ag alloy containing 

Au, we see that the c . is always greater than c, . except •«?., the end points, where 

s , Au 1 , Au 

they are equal. The segregation coefficient k^ u = ^ is thus always > 1. In 

fact, k varies from 1.000 to over 2.080, for the temperature range of 1336K or 

AU 

1063 C to 1234K or 960.5 C. Considering the same alloy as an Au -Jloy containing 

Ag, we see c . never exceeds c, . . Hence, k. - k. . . never exceeds unity. 

° s,Ag l,Ag Ag Ag in Au J 

In fact, k^ varies from less than 0.1986 to 1.000 over the same temperature range, 
representing an increase, of over 600 f /( (Li 1968). 


Li (1968) has fitted, by least-square techniques, the Au-Ag liquidus and solidus 
curves using straight lines passing through the Au terminal point, and also by quadratic 
or cubic polynomials through both the Au and Ag terminal points. According to the 
best-fitting equations presented, the constant segregation coefficient at the Au terminal 
point varies from 0.8781 through 0.7015 to 0.7287 for the linear, quadratic, and cubic 
polynomials, respectively. Thus, even at the same Au terminal point, the k = f^/a 
can vary by as much as 25.2% depending pn the degree of the polynomials chosen. The 
order of these k values is also inconsistent with the degree of the polynomials. 

Thus, we are often faced with a bewildering array of unique segregation constants 
including: k^, Fig. 2-3a; k £ , Fig. 2-3b K £T , Fig. 2-3c; k^ Fig. 2-3d; k , Fig. 
2-3e; k £u . Fig. 2-3f; and k^'s such as k oT k QT ^ and k^ rj . Any indiscriminate 
use of unspecified k is likely to confuse and mislead, particularly with melts of medium 
high solute concentrations such as 4. 5 % by Cu in A1 (a favorite alloy system for 
solidification study), for which the choice of k is the most ab undent. 


For any binary alloy system A-B in general, and the Au-Ag system in particular, 
the segregation coefficient of one component in the other, e.g. , of Ag in Au, is: 

k. = c . /c, . 

Ag s, Ag 1, Ag 
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The other segregation coefficient for Au in Ag in the same alloy is: 


k Au" C s,Au °l,Au " a_C s,Ag ) (1 " c l,Ag ) = (1 ' k Ag C l,Ag ) (1 ' C l,Ag ) ‘ 


1-K. L. (T) 
^Z.Ag 


Ti-l.JT) 


Ag' 


1 


If c . S c, . , then K. <1 and 1 - k. c, . M-c . Hence, k must be 
s, Ag 1, Ag Ag Ag 1, Ag 1* Ag Au 

> 1, as observed. Thus, if ' . _ < 1, then k . . must be >1. 

A in B o m a 


Further, since c. . = L 4 (T), k. cannot, from the last equation, be constant 

1. Ag Ag Au 

unless k. slsk. . This condition is tiue only when the liquidus and solidus are one 
Ag Au 

and the same curve, in which case Eq (27) reduces to the meaningless form c = c = c . 

s o 


Notice that all the above k's refer to a system of c 's and c 's on either the 

S X * 

weight or the atomic scale. That is, the k's are either k wt 's or k^'s. If can be shown 

that k cannot be the same as k 4 for any two alloy components unless c = c on both 
wt at si 

the weight and atomic basis, i.e. , only when both the liquidus and solidus are, again, 

one and the same. Further, if k is constant, k . cannot be constant, and vice versa. 

wt at 

Hence, both k. and k. in the Au-Ag system or k, . _ and k . in any binary 
Au Ag 6 A in B B in A 

system, whether on the weight basis or on the atomic basi r cannot be constant but must 
vary with the melt temperature. 

2.4.3 The Best Value of the Segregation Coefficient 

The best value of constant segregation coefficient k depends, of course, on the 

c 

alloy system. For a fixed alloy system, the best k £ varies, as shown above, with the 
selection of atomic or weight fractions for expressing the solute concentrations, the 
initial melt composition, and the degree and number of terminal restrictions in the 
polynomial fitting of the liquidus and solidus ci rves. A change in the least-square 
fitting method (to account, e.g. , for the fact that errors occur in both the solute con- 
centrations and melt temperatures rather than only in the concentrations or temperatures) 

also changes the value of k . The best k also depends on the use of the resultant 

c c 

solidified material. In the manufacture of critical electronic, optical, or mechanical 

devices, the objective may be to minimize either absolute or percentage errors in 

solute concentrations c, concentration gradients v = dc/dp, or other properties related 

P 
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to c or c . The minimization is to be achieved either over a specified region (i.e. , 

P 

between two fractions solidified, p or p ) or at a critical location (i.e., at the fraction 

Jl 4 

solidified, p). The best depends not only on the magnitude of allowable error, but 
on whether the error is positive or negative for the same magnitude. Specifically, 

I.i (1967) has shown that if overestimated c values cause more serious damage than 
underestimated c's, then a slightly inflated should be used for safety, especially 
when k and p are large. 

V 


However, it is always possible to find the optimum k^, for use in Eq (27) for 

example, to best St any given normal or abnormal freezing data. This is true even 

with alloys having distinctly curved liquidus and solidus lines. Further, the errors due 

to the assumed constancy of k c can usually be estimated either analytically or by 

numerical methods. In the very narrow melt temperature range around T^ T - T, 

for example, the optimal k can be determined from Eq (27) and (28). 

c 



where cr = l/(l-k ), and 8. and r are defined in connection with Eq (28). Hence 

c c 'op 

the optimized constant k c , ’.e., k^ , is such that: 




8 In 


1 + rT 
1 + rT 



(34) 


or k 


op 


°P , 
- 1 


l n 


T 


op 


- = 1 -■ 


ty tn 


I* tn 


1 + rT 
c 

rT 


for T - T 


(35) 
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In the above derivation, the liquidus and solidus curves are assumed to be 

adequately representable by quadratic functions of the melt temperature. In the 

general case when the liquidus and solidus are representable only by nth degree 

polynomials, the optimized k is 

c 


op - 1 
k °p = fc 

c op 


= 1 + 


for T ~ T 


V“~ + Ev n 

° i-2 


For a wide solidification range, the optimized may be determined by least- 
square techniques for the general polynomial case, such that: 




k° P -l 


[•■©]© 


k° P -l 


n / T-T. v 

■ V* 


In practice, the k v may often be more easily obtained by numerical methods 

C op 

using a computer. Table 2-1 lists the numerically obtained values of k^ for the system 

of Sb in Ge. The k° P is seen to increase with initial melt concentration c and fraction 
c o 

solidified p, the increases being greater with larger values of c q and p. The standard 

errors in the computed p, computed with these k° P according to Eq (27) and compared 

c -12 

with corresponding p values from Eq (28), can be as small as 3 x 10 for c/o = 

0. 001 and p = 0. 01 to 0. 20. 
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TABLE 2-1 OPTIMIZED CONSTANT k, k° P FOR Ge ALLOYS 


CONTAINING Sb AS SOLUTE 





2.4.4 Errors Due to Curvature of the Liquidus and Solidus Lines 


We consider here the errors involved in the use of the constant-k, normal freezing 

equation when, in fact, the liquidus and solidus curves are not straight lines. Three 

cases are considered: quadratic, high-degree polynomial, or exponential functions of 

the melt temperature T . 

m 

In the first case, the liquidus and solidus curves can be represented by quadratic 

function of T , i.e. , 
m 

2 2 
L(T) = c„ (T) = aT + bT = a(T -T ) + b(T -Tm) 

' 1 erne 

2 2 
S(T) = c (T) = fT + gT = f(T -T ) + g(T -T ) 

^ s ' e m ° e m 

The segregation coefficient k used in Eq (27) should, particularly for such dilute 
alloys as semiconductor germanium or silicon and pure Be or Ti, be: 

k o = S(T)/L(T) Ts0 = f/a 

Equation (27) then shows that the proportion of liquid remaining is: 


1 , 1 , „ .l/(k-l) , . l/(k-l) a/(f-a) 

P, = 1-P = (cAc ) ' = (c,/c ) ' = (T/T ) 

*1 s o ' 1 o o 


According to Eq (28), the proportion of liquid remaining is, however, 

B 


a_ 

f-a 


a 

J*-r 


4 ■ ■ ( T /T/ /(t - a> p^:|>^ t ) ] 1 jl + I tS (T -T o) J,^ '[l-r,T-T o ,j 

where B = (ab+ag-2bf)/(f-a)/(b-g) and r = (ab+ag-2bf)/(a-f) 2 
2 1 

To estimate by means of p^ of Eq (25) thus involves percentage errors 
amounting to: 
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(38) 


12 1 2,2 

e = (P 1 ~P 1 )/P 1 = 


oKT-T o ) 
a + o, (T-T q ) 


= o(T-T ) = ^ - + - 2bf A T -- IT 


(a-ff 


This error is thus roughly proportional to the solidification range AT = T-T q . 

Estimating the amount of eutectics thus represents the worst-case condition for the use 

of the constant-k normal freezing equation, Eq (27), particularly when initial eoneen- 

12 

tration c q is small. Error f is also inversely proportional to the square of the slopes 

of the liquidus and solidus lines, i.e., (a-f) 2 . When both a and f are small, i.e. , when 

the constitutional diagram has very steep liquidus and solidus curves and, in particular, 

2 

when a-f is small, l/(a-f) is very large and, hence, very large errors may result. 

Examination of Hansen's (1958) binary constitutional diagrams discloses many 
systems whose liquidus and solidus are very close at the melting point of the pure 
solvent (i.e. , k Q = 1. 0), or are very close at the eutectic temperature when eutectic 
growth is of prime interest. For example, binary alloys having k Q ^ 1 include silver 
alloys containing Cd, Li, and Zn; gold alloys containing Ag, Al, Cd, Co, Cu, Ga, Fe, 

Pd, and Zn; iron alloys containing Be, Co, Mn, Ni, Si, and V; and nickel alloys of Cr, 

Mn, Ta, Ti, and V. In addition to k Q «s 1, many systems have liquidus and solidr 
adequately representable only by cubic or higher degree functions of the melt temperature. 
Such systems include Ga in Au, Cd in Mg, Mg in Cd, Co in Pd, Mo Cr, Cr in Pd, Cr in 
Ti, Cu in Pd or Pt, Mn or Ni in Fe, Cr, Cr, Fe, and Ni in Pt, and Ta, Ti, and V in 
Ni. 

-3 

As an example, consider the case of Sb in Ge. Here, a = 4.668 x 10 , b = 

—fi A —ft 

-6.047 x10 , f = 1.288 x10 , and g = -8.222 x 10 (Li, 1967). Hence, 

e 12 = -0. 00133 x (T-T ) = -0. 133 A T % 

o 

The maximum solidification range to be considered occurs when a semiconductor- 
grade germanium starts to solidify at about the melting point of pure germanium (i.e. , 
1209K or 936°C) but finishes the solidification at the end of a rod or at the boundaries 
between subcells or subgrains at the eutectic temperature of 863. IK or 590°C. In this 
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case, T q = 0 and AT = (1209-863. 1)K = 346. IK. The maximum error in the growth 
of Sb-doped germanium crystals is thus 

e 12 = -0. 00133 T % = -44. 7% 


That is, the constant-k equation consistently underestimates the actual p 1 in the Ge-Sb 
system. This underestimation may reach a maximum of about 44. 7%. 

As another example, consider Co in Ti. According to Hansen's (1958) diagram, 

-2 -5 -2 

coefficients a, b, f, and g are 2.3552 x 10 9 4.3175 x 10 9 1.4624 x 10 9 and 

6.7989 x 10 , respectively. The value of a may then be computed to be -0.10764%. 
Since the maximum solidification range is T = 908. IK, the maximum error is 

e 12 = oAT =-68.35% 

max max 


Thus, the amount of the critical eutectics at the grain or cell boundaries in these 
alloys, estimated using the constant-k equation, is always lower than the true value by 
up to 68.35%. 


Often, the liquidus and solidus can be represented only by nth-degree polynomials 
of the melt temperature, i.e. , 


n i n 
s, =£. aT 1 =T a (T 
l i=l i i=l 1 


T ) 
e m 


n i n i 

c f, T = £ f . <T -T > 
8 i=l i j=i » e m 


The proportion of the remaining liquid is then given by Eq (29): 

/ Tv 

1' * >1 \ T „-V '( T „) ( T o) 


n q.- 


1 + (T-T.)E 


i i=2 T 


where T^ and q^ have been previously defined in connection with Eq (29). 
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Hence, the percentage error involved in the estimate of by means of the 
constant-k equation (Eq 27) when the liquidus and solidus are high-degree polynomials 


is: 

n i 

-r t 

In . 1 n. , n i=2 
e = (Pj-P^/Pj = 


n Q, 


n q i 


n 

1 + (T-T ) £ T 

o i=2 i 


— (T -V“ i n i 2 n tT i39 » 


2 i 


which is again proportional to solidification range A T = T - T q . 

For ideal solutions, the molar concentration for the solvent A atoms is 


In x. = 
A 


tH A /1_ J_\ 
R (t -T A j- 


4H » / T V* ‘‘"A 

< T ' =7 


H 


T_ _T 

T . + 2_ 2 
A T, 


+ . , 


The molar concentration for the solute species is 


X B = 1 ~ X A 


AH A ^H 9 n i n , 

— T + — t + . . . * E a.T or T f.T , 
RT, RT. n=i ‘ 1-1 1 


for the liquidus or solidus lines, respectively. 

Hence, an exponential function of the melt temperature in the liquidus and solidus 
equations can be approximated by appropriate polynomial functions of the melt temperature. 
The percentage errors involved in the use of the constant-k equation can, therefore, be 
similarly estimated or computed. 

2.4.5 The Last Portion of Freezing Melt 

Chalmers (1969) discussed the fqct that normal freezing of an alloy melt results in 
the last material solidifying at concentrations higher, sometimes considerably so, than 
that of the original liquid. In 1939, Hayes and Chipman stated without proof that the last 
portion to freeze may be a eutectic. 
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To understand why the last portion of a melt to solidify under nonequilibrium 
conditions always takes place at eutectic temperature T y (Fig. 1-13) yielding, in- 
variably, the eutectic, we can intuitively reason as follows. We consider only the 
extreme case of a very dilute melt whose initial solute concentration is =*0. During 
one step of the freezing process the melt temperature decreases by perhaps AT m = 

- AT = 10K or IK, with an always incomplete or fractional loss by solidification, i.e., 
Ap or Ap/(l-p) respectively, of the melt. Simultaneously, both the freezing solid and 
remaining melt are solute-enriched according to the solidus and liquidus curves (for 
the usual case when k < 1, as shown in Fig. 2-2). This stepwise nonequilibrium 
freezing process may be repeated many times, but eventually the melt temperature 

T must come down to T . At T , there must still be a finite but perhaps very small 
m u u 

amount of the liquid melt left (since a fraction of a small fraction is still a finite 
fraction). This liquid, being at T^, has the eutectic composition = c^ u « 

One can also obtain a proof of this conclusion from the normal freezing equations. 

Since at T = T , c = c , the proportion of the normally-frozen eutectic P is from 
mussu u 

Eq (27), for the constant-k case with linear liquidus and solidus lines. 


1 

k-1 



When the liquidus and solidus curves are quadratic functions of the melt temperature 

i.e. , c = a (T -T ) + (T 
t e m' ' t 

normal freezing equation gives 


T , i.e. , c = a (T -T ) + (T -T ) 2 and c = f(T -T ) fg, (T -T ) 2 , the appropriate 
m t e m e m s e m 1 e m 



"(a-f) + (b-g) (T e -T u )] q 
_(a-f) + (b-g) (T e -T.) _ 


(41) 


where r and q are constants for a given alloy system, as defined by Li (1067a); and a, 

b, f, and g correspond, respectively, to A , b . a , and b . 

f f s s 

Since c , c , k, T , T , and T. all have real, positive values, and since a, 
su o e u i 


2-33 



f, b, g, r, and q are all real constants for a given alloy system (with a > f | b | 
s* | g | ), P ; must therefore be real, positive, nonzero, and calculable. 

Similar mathematical proofs can be obtained for cases where the liquidus and 
solidus curves are cubic and higher-degree functions of the melt temperature (Li, 1968); 
are exponential functions of the reciprocal absolute melt temperature, i.e. , ideal 
solutions (Li, 1971); or are portions of a ternary constitutional diagram of the pseudo- 
binary type (Li, 1972). 

It should be emphasized that determining p y from Eq (27), particularly when c q 
is small, represents its worst-case usage. As will be pointed out, such a usage is not 
only pedagogically incorrect but may yield erroneous and completely misleading results. 
Such not uncommon usages (see, e. g. , Flemings and Mehrabian, 1971) are, therefore, 
to be discouraged. 

Mathematical proofs with the aid of the non-constant-k, normal freezing equations 
are, therefore, necessary for most cases where the liquidus and solidus lines are 
curves rather than straight lines. The intuitive proof given above is particularly desir- 
able because it is universally true and applies to all normal freezing conditions, regardless 
of the types and shapes of the constitutional diagrams. All these proofs are, however, 
based on the assumptions of zero solid diffusion and infinite liquid diffusion. While 
solid diffusion is negligible in most freezing studies, as will be more full discussed, 
the liquid diffusion in actual freezing processes is never infinite. This condition results 
in a solute-enriched layer next to the solid-liquid interface (Tiller et al, 1953). The 
higher solute concentration in this layer, rather than the lower concentration in the main 
body of melt, now determines the concentration oi the freezing solid. That is, the 

effective (variable) segregation coefficient k is now increased toward unity. However, 

6 

k never reaches unity at finite growth velocity (that is, solute rejection by the freezing 
© 

solid into the remaining liquid is unavoidable) and, above all, there is always a final 
solidifying portion in which the freezing conditions are nearly normal (Tiller et al, 1953). 
Hence the last solidifying eutectic weight or volume is never completely eliminated. 
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The above discussion thus shows that nonequilibrium freezing always results in 
the formation of a finite amount of eutectic at the eutectic temperature no matter how 
small the initial solute concentration. 

2.4.6 Numerical Computation of Normal Freezing 

The array of normal freezing equations presented above does not cover all con- 
ceivable types of constitutional diagrams. Often, the liquidus and solidus lines are not 
simple polynomial functions of the melt temperature. Computations, using closed-fom 
equations, of the normal freezing behavior directly from basic thermodynamic constants 
are impossible for cases where the liquidus and solidus lines are not simple, logarithmic, 
or exponential functions of the melt temperature. For example, Thurmond and Kowalchik 
(1960) has shown that, when the heat capacity of the pure supercooled liquid solvent is 
constant but that of the pure solid solvent is linearly related to the melt temperature, 
the logarithm of the solute concentration in the liquid has three additive terms in addition 
to the one previously given for the ideal binary solution case. Here, the normal freezing 
differential equation cannot be solved exactly. Reisman (1970) has presented a large 
number of complicated equations for the liquidus lines, the normal freezing differential 
equation of which cannot be solved exactly and numerical methods must be used. 

In most cases, the liquidus and solidus lines can be adequately represented by 
polynomial functions of the melt temperature. Yet, determination of the roots of the 
polynomial functions for insertion into the closed-form, general normal freezing equation, 
Eq (30), must be accomplished numerically if the polynomial functions are of quatic or 
higher degrees. Hence, there is a definite need for a reliable numerical methods of 
known accuracy that can be applied to any type of constitutional diagram. 

A literature search failed to reveal such a numerical method. Gulliver in 1922 
and Brody and Flemings in 1966 presented the same numerical method for variable-k, 
normal freezing computations. This is a combined graphic and finite-difference 
(Butler and Kerr 1962) method that approximates the liquidus and solidus curves by a 
series of straight-line segments. Pfann's equation, Eq (25), is repeatedly used, once 
for each pair of the segments between two consecutive isotherms. This method is 
simple and convenient, but its accuracy is limited by the large (at least 5 or 10%) 
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errors introduced when graphically estimating the local liquidus and solidus slopes 
and by the repeated error accumulation or buildup through the many computation steps. 
The accuracy, therefore, cannot be high or even satisfactory for critical solidification 
work, no matter how small the size of each computation step. Further, the accuracy 
has never been compared with closed-form, normal freezing solutions, other refined 
numerical methods, or even the same method with reduced computation step sizes. 

Note that Gulliver (1922) did disclose a useful method for metallurgists at a time when 
numerical methodology was not advanced and when modern computers or even calculators 
were not available. 

Two numerical computation methods are briefly described below. These are 
the Runge-Kutta method and a finite-difference method. 

2 . 4 . 6 . 1 The Runge Kutta Method 

The Runge-Kutta method (Kunz, 1957) can be used to numerically integrate the 
normal freezing differential equation and to produce the final, computed value of fraction 
solidified p at any melt temperature T. This method is stable and self-starting; that 
is, only the functional values at the initial, or a single previous, point are required to 
obtain subsequent functional values (i.e. , p and T). For this reason, we can easily 
change the step size h at any step in the computations. This is particularly desirable 
when the computed results must be more accurate in some solidificat on ranges than in 
others. 


The ordinary differential equation for normal freezing, given as Eq (25), can be 
written; 


dE _ (1_-P) L’(T) 
dT L(T) - S(T) 


= f(T,p) 


where L'(T) is the temperature derivative of L(’i ), L’(T) is the temperature derivative 
of the fraction solidified p(T), I.e. , dp/dT, and f(T,p) is a new function of both T and p. 

For the case with quadratic liquid’ s and solidus lines, for example, the above 
equation becomes: 
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dp/dT = f(T,p) = 

(a-f)T + <b-g)T 


(42) 


The initial conditions are p * 0, T = T , L(T ) = c , and S(T ) = c (Fig. 2-2). 

The above normal freezing equation is solved numerically by using a single-step, 

fourth-order Runge-Kutta integration process (Kunz, 1957). In this process, the value 

of f(T,p) at T = T is used to compute f ,(T,p); earlier values f , , f ., etc., are 
n n+i n-i n-4 

not used. 

The relevant equation is: 

f n + l =f n + <k l +2K 2 +2 V k 4 ,/6 < 43 > 

where, for step size h; 

k =h f (T ,p ) 

1 n n 

k =hf (T +h/2, p +k /2) 

4 n n 1 

k = h f (T +h/2, p 4k /2) 

3 n n 4 

k 4-hf( V h.P n ^» 

2. 4. 6. 2 The Finite-Difference Method 

The basic differential equation for normal freezing, Eq (25), is 


J*E = 

1-p 


-din (1-p) 


L’. 1T) 

L(T) - S(T) 


dT 


where C = L(T) and c = S(T) are function of the modified melt temperature for the 

F S 

liquidus and solidus, respectively. L’(T) is the first derivative of L(T) with respect to 
the modified temperature T and can be found from L(T) = c^. 

During a short time interval At, T, c = L(T), and c = S(T) are all fairly con- 

v 8 

slant, and usually both c and c are increasing. Their difference, Q(T) = L(T) - S(T), 

f 8 

can also be considered as constant, which is the essence of the finite-difference method. 
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Hence, 


In (1-p) 






L'(T) 

Q(T) 


dT = - 


1 

Q(Tj) 


MT) 



1 

Q(Tj) 


[ L(Tj)- L(T i+1 > ] 


L(T W 

Q(Tj) 


(44) 


where 

MT t ) 
Z i “ Q0\) 


(44a) 


Knowing at T ■ T^, we can estimate p J+1 at T = T (so that AT = T J+1 - T.^), 


as follows: 


P, +1 - 1 - (1-Pj) exp 



(45) 


Equations (44) and (45) allow us to compute the complete normal freezing behavior 

of an alloy having an initial solute concentration in the melt c = c . The starting 

o to 

point for this computation is the point at which the liquid begins to greeze, i.e. , 

p = 0, T = T , c , = c, = c , c = c , and Q(T ) = L (T ) - S(T ). From these known 
o o l lo os so o o o 

values of parameters, Z q can be counted by Eq (44a), from which 

Z - MT o ,/Q(T ) V, 
lo BO 

We can then specify a constant (or variable) value of incremental temperature 
AT * - Tj, and And the successive values of T. From these T's, we can compute 

the values of L(T), S(T), and Q(T) » L(T) - S(T) for substitution and use in Eq (45). 
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Note that, in deriving Eq (45), the value of L(T), S(T), and Q(T) = L(T) - S(T) 
are assumed constant during the very short time interval At or temperature interval 
AT ■ T - T . This represents the zero-th order approximation of polynomials L(T) 
and S(T). This approximation is selected to integrate the basic normal freezing 
equation, Eq (25), for use in the finite-difference method. As will be shown, this 
method gives useful results when compared with exact solutions from Eq (28) - (30). 

Another way to integrate Eq (25) for use in the finite-difference method is to con- 
sider S(T) constant during the very short time or temperature interval. We then have 

- -dJn (1-p) dT 55 din [L(T)-5(T)3 (46) 

The results, computed according to Eq (46) are, however, no better than those 
obtained from Eq (44) - (45); they are therefore omitted from further discussion. 

Added accuracy in the computed results may be achieved if we replace Q(T { ) 

in Eq (44) and (45) by Q(T ), the average of the two values of Q(T) at T = T and 

1 | 1^1 1 

T = Tj + ^, respectively. We therefore have a choice of the following two computation 
schemes, which are arranged in the general order of decreasing accuracy in results but 
increasing ease in usage: 

I. Equations 44 and 45, together with Q(T { j+ ^) and Q(T j+ ^ i+ g) 

II. Equations 44 and 45, together with Q(T { ) and Q(T^). 

2.4.7 Computation Example 

The following computations illustrate the above procedures. The alloy melt con- 
sidered here consists of 10% Sb in Ge, i.e. , c = 0.10. 

o 

The liquidus and solidus curves, L(T) and S(T) respectively, for the Ge-rich 
portion of the Ge-Sb system have been determined and represented by second-order 
equations, from which we find: 
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22. 05K 


T = 
o 

c = G(T ) = 0. 00260 
s o 

o 

c 

Z = 1.03 

o c - c 
0 s 

0 

Table 2-2 gives the computed values of the fraction, p, of solidified melt at 17 
different melt or interface temperatures, for 10% Sb in Ge. Both computation uchemes 
were used, together with a AT of from 10K down to 0. 0001K in five decades. The number 
of sets of computed results thus varied from 32 (for AT = 10K) to 324,104 (for AT = 0. 001K). 
Only a small fraction of these recults are tabulated. Table 2-2 also includes the corres- 
ponding values of p computed with _ closed-form equatior derived for cases where both 
the liquidus and solidus curves can be expressed as quadratic equations o' T, Eq (27), 
such as those given above for the Ge-Sb svstem. 

In comparing the p values at a particular T, wo notice that: 

• The two computation schemes gave results close to the corresponding values 
obtained from the closed-form solution. This indicates that the numerical 
method is useful, at least for 10% Su in Ge 

• Scheme I gave much better results than scheme II. Thus, the . Jue of 
averaging L(T), S(T), and Q(T) during the computations is substantial, 
and averaging procedures should generally be used. 

• Decreasing AT generally improved accuracy or decreased deviation from 
the closed-form solution values, as would be expected 

• For a given scheme with a specified AT, accuracy consistently improved 
with the decreasing melt temperature, Tm. 

These conclusions are quantitatively disclosed in Table 2-3, which gives the 
average and standard deviations from the p's by the closed-form solution, for the 
different schemes. This table shows that for 10% Sb in Ge: 
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• Minimum discrepancy from the closed-form results for each scheme was 
reached with a medium AT. Scheme I reached minimum standard deviations 

_7 

of about 3x10 with a AT of 0. IK, while scheme II had minimum standard 

-5 

deviations of about 1x10 with a AT of 0. OIK 

• Except with a AT of 0. 001K, scheme I generally gave results one to over 
two orders of magnitude more accurate than corresponding values obtained 
from scheme n 

• The average deviations of the results from p (from Eq (27) computed with 
scheme II were consistently positive; that is, this scheme usually over- 
estimated p. Scheme I, on the other hand, changed from underestimation 
to overestimation as AT decreased from 10K to 0. 001K 

-7 

• Average deviations were as low as 1.6 x 10 for scheme I with AT = 0.1K. 
Even with AT = 10K, the same scheme had an average deviation of only 
-1.3 xlO -3 . 

Knowing the standard deviations and number, n, of results per set, one can 

easily compute the maximum probable deviations for any single result. Thus since 

-5 

scheme I had an average deviation from p of -1. 34 x 10 and a standard deviation of 
—6 

6. 50 x 10 with a A T of 1. OK and n of 324 , 95% of the time one would expect any single 

-5 

numerically computed result to be accurate within the range -1. 34 x 10 + 2. 0 x 6. 50 

_g —7 —5 

x 10 or -4. 0 x 10 and -2. 64 x 10 ; that is, under these conditions the numerically 
computed result is almost always underestimated, by from 0.40 to 26.4 ppm in 95 
times out of 100. 

The above finite-difference procedures are well- adapted to rapid, automatic 
machine computations. Further, all these procedures apply for any liquidus or solidus 
curves that can be represented by reasonable mathematical functions, F(T) and G(T), 
respectively. These functions may be power series of T and include, as special cases, 
the simple straight-line (constant k) or quadratic equations shown above. These 
functions may also be exponential, logarithmic, or otner types. Even several different 
functions for the same liquidus or solidus curve, but in different freezing ranges, can 
be handled with nearly equal facility and computational speed. 
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PART III 


NUMERICAL PROGRAM 
Section 1. INTRODUCTION 


To understand solute segregation and defect formation in Single crystal growth 
thoroughly, one requires a complete characterization of the (mass) diffusion and 
temperature fields in the solid and remaining melt. The zero-gravity effect on solidi- 
fication appears to be subtle and may be overshadowed by other effects invariably in 
any such growth process, a condition necessitating that such characterization be 
accurately defined. Unfortunately, the coupled partial differential equations for the 
diffusion and temperature fields are generally not solvable. Although special-case 
solutions have been given for some types of physically nonsatisfying two-phase Stefan 
problems, we must resort to numerical computations for the general case solution. 

E xisting numerical methods are always subject to such unrealistic assumptions as 
constancy of interfacial temperature, segregation coefficients, diffusion constants, and 
other material thermophysical properties. 

A program for the numerical solution of unidirectional solidification (or crystal 
growth) from a semi-infinite binary alloy (with nonconstant segregation coefficient) has 
therefore been initiated. We numerically determine the interphase temperature, and 
temperatures and concentrations of the mesh points, for the interphase point after each 
time interval. The interphase temperatures and the temperatures and concentrations 
for the two adjacent mesh points are solved, iteratively, only after convergence has been 
obtained. The temperatures and concentrations of the remaining mesh points are deter- 
mined by solving two tridiagonal systems of linear equations. The semi-infinite nature of 
the problem is handled by admitting as mesh points enough points, so that conditions at 
the penultimate (next to last) point have not changed significantly from the original ones. 

A scheme with variable time intervals is adopted to achieve convergence within reasonable 
computing time. 
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We have checked out the program for a binary alloy. The results compare favor- 
ably with those of the exact analytic solution, differing mainly in that the interphase 
temperature is not constant as assumed by the analytic solution. 

Section 2. MATHEMATICAL DEFINITION OF THE SOLIDIFICATION PROBLEM 


We consider here a liquid binary alloy directionally solidified into two phases, 

liquid and solid. We consider the liquid alloy to be semi-infinite with original (at t = 0) 

temperature T and concentration C . Solidification occurs when the temperature at 
o o 

x = 0 is changed from T q to a lower value T^ either instantaneously or gradually, 
such that T^ is below the temperature T^ at which the liquid mixture at concentration 
C q can be in equilibrium with a solid phas". As solidification occurs, the solid phase 
grows and its boundary is given by x = y(t), the interface point, and the interface tem- 
perature at this point is T. (t). The partial differential equations describing the solidi- 
fication process are: 


a 2 T 

s 

aT 

s 

a 2 c 

n s 

ax 2 

at 1 

8 ax 2 

aV 

3T. 

a 2 e 

i 

1 

i 

• n ^ 

>x 2 

at 

’ °t . 2 
a x 


a <.TT’Tr iD t7~2 = — for y(t) < x < - ( 

a x ax 

The variables T and C represent the temperature and concentration (of solute 
in solvent) and the subscripts t and s denote the liquid and solid phases, respectively. 
The quantities A s » A^, D^, and D^ are assumed constant. The following conditions 
are assumed throughout: 


T (x, 0) = T and C. (x, 0) = C 
x o i o 


T.(», t) = T and C (®, t) = C 
X ox o 


I 


I 



T g [y(t), t] = [y(t), t] = Tj(t) (49c) 

C s [y(t), t] = f s [T 1 (t)] (49d) 

C x [y(t), t] = [T t (t)] (49e) 

ST ST 

PYy (t) = k 8 57 - k ^ for x = y(t) (49f) 

ac ac 

| f s 1 - r 1 ( T , <t >]} = D iir- D s TT torx ' « 49 «» 


Equation (49a) specifies that the original mixture is all liquid at temperature T and 

o 

concentration C q . Equation (49b) is a consequence of the semi-infinite nature of the 
mixture so that, at any time t, the portion near infinity is unchanged. Equation 
(49c) assumes that there is an interface temperature T { (t) at the solid-liquid interface 
plane and that both the solid and liquid phases at x = y(t) have this temperature; there 
is no jump in temperature. Equations (49d) and (49e) state that the concentrations of 
solid and liquid at the interface are given by the solidus and liquidus curves, respectively, 
of the constitutional diagram for the alloy. Equation (49f) connects the derivative of the 
moving boundary with the redistribution of temperature, and Eq (49g) connects the same 
boundary with that of the concentration. 

The conditions on T (o, t) and C (o, t) are not fixed in our discussion and a number 
s s 

of alternatives are considered: 

• T (o, t) = T (t), with T (t) equal to a constant for all t, 

SIX 

• Unear, T^t) = T q + M^-T^/b for t<s and T^t) = T for t * s, 

• Exponential, T (t) = T + (T - T) e ,soT,(0) = T and T, (») = t. 

xiox 101 1 
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For C (o, t) the conditions considered are C (o, t) = c , usually taken C (T ), 
s Cl 8 2 

or at times a condition conserving mass between 0 and 00 . 

Our two approaches may be designated as analytic and numerical. The numerical 
approach can be applied to all three conditions on temperature, whereas the analytic 
approach holds only the case of constant temperature instantaneously applied. A variant 
of the analytic method, i.e. , applying linearly varying temperature, should be investi- 
gated. 


Section 3. THE ANALYTIC SOLUTION 
A similiarity solution can be written m the form 


= 2B^r 

(50a) 

= T 1 + (Tj - T^ erf (x/2a^)/erf Jy(t)/2a^f| 

(50b) 

= T q + (T t - T o ) erfc (x/ 2&jjt )/ erfc jy(t)/2a^J 

(50c) 

= c i + ^(Tj) - Cj erf (x/^/Dl)/erf jy(t)/^jD^Tj 

(50d) 

= c 0 + - C 0 ]erfc (x/?JC^t)/erfc [y<t)/3jDjtj 

(50e) 


The quantities T^ and fi are determined from the boundary conditions. It is clear that, 

if T, sT = f ,* (C ), then there is no solidification, C (x, t) = C , and T (x, t) = 

1 2 c o i o 1 

T + (T - T ) erfc (x/2a.+Jt). The solution given above satisfies boundary conditions 

O 1 O v 

(49a, b, c), and equations (49f) and (49g) are employed to determine B and T^ as 
follows (Rubenstein, 1971): 

y = B/y/t (51a) 

2 

x 

dT 2 

sr- ' ,T i - T i>^ 488 ' • igr /erl K 2 vf] ' 5lb > 


3-4 



(51c) 




V 


2 

X 




>c 

sr ■ [w - C J ^= 6 4Dst K 2 Vv] 


(51d) 


jC - ~ 

5T“[W- c o]- -ft e ^ ^ /erfc Ksp?] 


(51e) 


Substituting in Eq (49f) and (49g), eliminating the factor , and replacing x by y(t) and 

f 


y(t) by 2 a JT7 we get 


B 


B 


PvB = 


2 2 

k (T. -T )e a ® k. (T -T )e ^ 
si 1 -c i o 

a^n erf (B/a g ) erfc (B/a 2 ) 


(51f) 


[f 8 <Tj) - fj(Ti)] 6 = JTTerfe 


By rewriting Eq (51f) solving for T { , we get 


a 2 2 
"d" - 1 - 

sF,f (f t (T i>- c J e ' >FIf f s ' T i> - c i] e Ds 


Jnerf (B/^IT) 


(51g) 
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3.1 NUMERICAL METHOD FOR THE ANALYTIC SOLUTION 


The successful method for solving the pair of equations (51g) and (51h) involved 
starting with an arbitrary (usual small) value for 0. Substituting in Eq (51h) gave a 
value for T., and using the values for 0 and T. in the right-hand side of Eq (51h) and 
in the multiplier of 0 gave a new value for p . This value was checked against the 
previous value, and if they agreed to a desired accuracy the solution was accepted. 
Otherwise, the two values were averaged and the procedure repeated. Failure of the 
procedure to converge after a reasonable number of iterations suggested that a dif- 
ferent starting value be tried. Sometimes, iterates which failed to converge suggest a 
better starting value: if they diverge immediately try a smaller value, and if they 
oscillate try a value close to the point about which they seem to oscillate. 

A modification of the analytic method, based upon the assumption that 0 need not 

be constant and that T also varies, has been included in the program. The assumption 

* 2 i 

is that y(t) = 2 0(t)^7, where p(t) = + Pj* + 0£* + • • • 0j ^ and 


T.(t) = T. + tT., + 
i 10 if 


+ t J T 


(52) 


This procedure was originally used to determine whether solutions of this more complex 
form existed, but was inconclusive because the values of the higher-term coefficients 
were sufficiently small to be explained away by computation error. However, this pro- 
cedure may be useful in connection with a solution to the more realistic problem where 
T^(t) is variable. This is yet to be investigated. 

The method of solution follows lines similar to those given above. 


y (t) = ^ + 33l^ + ... + (2 H)0jt = 7= 


where 


0. = P +30,. t + +(2J+1) 0. t J 

d o J 


(53) 
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First, assume initial estimates of 0 Q , 0^ ...... 0^ and values of t = 0, &t, , 

j At. Next, determine values of T { (t) and 0(t) by substituting into (51g) and (51h), taking 
care to replace 0 by 0^ in those terms obtained from the derivative. Using *he values 
of 0(t), we can find the coefficients by solving a linear system; these coefficients can 
then be compared with the previous estimates. If they agree, we can find the values of 
Tyt) for the coefficients T. q , . . . . , T . If they do not agree to the required accuracy, 
the coefficients of 0(t) are averaged and the process repeated. To find estimates for the 
starting values 0 q , 0^, . ..0^, we first solve for j-1 and use the coefficients so obtained 
as a first estimate for j, with g^ set equal to -0^ . 001. Thus, we solve a sequence 

of problems j = 0, 1, . . . , up to the desired value of j. The numerical results using this 
technique has been good, although its application to the case of constant T haB not been 
fruitful. 

The only boundary condition for C (o,t), discussed above, has been C (o, t) = C. » 

s S A 

a constant for all t. A variant condition is one based upon conservation of solute: 


q < C s- C o )dx+ J 


(C -C )dx = 0 
..to 


Using the expressions for C and C in the analytic solution (Eq 50 c, d), we have 


f(T)-C_ ( J "' f.(T.) - C 

(c r c 0 ,y<,) * a fr&jn, J + eTferkiffi 


1 o' 


erf (y/2. 


S x 


erfc (y/2JDl) 


QO 

f erfc (x/2jo” )dx = 0 

Jy(t) 

U 2 

f 1 -u 

j erf v d V = u erf u + (e -1) 


and 1 erfc v d v = - u erfc u + — ■==- e U 

\ -F 
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we obtain the relation 



By comparing this condition with Eq (49g), we have 


(57) 


[ f s <T i> 


V T|I ][ y(t) - 2ty(t) 


. n 1 

y(t)] -jf 


f <T.) - C 

8 ' .2J5T 

erf (y/^jDt) N s 


(58) 


For the case where y(t) = 2P,jr, y(t) = 2ty(t) and the condition is f (T ) - C = 0. This 
means that the concentration of the solid solute at x = 0 is the same as the concentration 
at the interface, and hence the concentration of the solute in the solid is constant every- 

ac 

g 

where. This is equivalent to the condition given by Rubinstein (1971): — — = 0 at x - 0, 

OX 

» 

which implies that C, = f (T.). However, if (y(t) - 2ty(t) is not zero, f (T.) - C can be 
1 s i s i 1 

computed to determine at every time interval. 


The program for computing temperature and concentrations for the solidification 
of a binary alloy system is named BISCL. BISCL employs a subroutine BETTI to deter- 
mine the values of $ and T^. The main program computes, for each instant of time t, 
y(t), y(t), T , f (T.), and f (T ) and then, for each (x, t), T (x, t) and c(x,t). The spacing 

1 8 1 1 1 j_i 

for x may be equal; j A x for j = 0,1,.... Ix-1 given by IAM = 0, or 0, 2 ax for j =1, 

. . . , Ix-1 where a x given by DELX and IAM = 1. To compute f (T ) and f (T ), the point3 

Si li 

on the solidus and liquidus curve are input as a table of triplets (T, f , f ), and 

S 1 

Lagrangian interpolation is used to interpolate between tabular points. The number of 
interpolating points L is less than or equal to the number of triplets IF. The time inter- 
vals are equal, given by DELP, and are terminated when time exceeds TF. The number 
of iterations permitted before the program was stopped for further iteration is given by 
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NIT. The number of coefficients of 8 (t) = + Bj t + . . . . t^ is given by KT = j +1. 

In the ordinary case, KT = 1. To determine the first approximation of &, the quantity 
FRX is multiplied by DELX. In the case of a failure to converge, an error indicator is 
set in BETTI and tested in the main program permitting a change in FRX. If FRX is set 
to zero, the program will stop. As a boundary condition at C (o,t) if C. >0 is input, this 
is held constant. If = 0, then it will assume a conservation of solute concentration 
and compute C. in the program as required (see above), a is designated by AS, a by 

1 8 Jr 

AL, D by DS, D by GL, T by TO, T, by Tl, C by CO, K by KXS, k by XKL, p by 

S l o 1 O S l 

RHO, and V by GAMMA. Quantities needed by BETTI and BISCL are transmitted by a 
common statement. 

3.2 GLOSSARY 

The following is a glossary of additional terms used by the program. 


SQDS 

- 


FAST 

- 

erf (YAST) 

SQDL 

- 


YDST 

- 

YT/DST 

SQT 

- 


FDST 

- 

erf (YOST) 

BETA 

- 

s 

YALT 

- 

YT/ALT 

BETH 

- 

9 d 

FALT 

- 

erfc (YALT) 

TMIN 

- 

T 1 

YDLT 

- 

YT/DLT 

YT 

- 

y(t) 

FDLT 

- 

erfc (YDLT) 

DYDT 

- 

y(t) 

FLT 

- 

f ^ (TMIN) 

AST 


2 D „>F 

FST 

- 

f (TMIN) 
s 

ALT 

— * 

2 a «>F 

FLTC 

- 

FLT -CO 

DST 


M ~ 

FSTC 

- 

FST-C1 

DLT 

“ 

2 ^T 7 

TEMP 

- 

T(x,t) 

YAST 

- 

Y T /AST 

CN 

- 

C(x,t) 
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Section 4. THE NUMERICAL APPROACH TO SOLIDIFICATION 

In order to solve the solidification problem numerically, we devised a technique 
for separating the two major aspects of the computation. One aspect involves the deter- 
mination of the moving interface boundary, x = y(t), and also the interface temperature, 

Tj; the other aspect concerns the determination of the temperatures, T, and concen- 
trations, Cg and C on both sides of the interface. Although it is clear that the cwo are 
interconnected for computation purposes, we first extrapolated for a new set of tem- 
peratures and concentrations at the mesh points, then determined the interface point 
and temperature, and finally recomputed the temperature and concentrations based upon 
the interface point and temperature. In case of a failure to converge quickly to the inter- 
face point and temperature, we interposed new values of the temperatures and concen- 
trations at the mesh points to help expedite the process. With this general strategy in 
mind, we shall explain in detail the various aspects of the program. 

4.1 THE INTERFACE PROCEDURE 

At time t, we computed the values of T(x, t) and c(x,t) at the mesh points and the 
values of y(t), y(t), and T^(t). In addition, we obtained the values of T(x,t + it) and 
c(x,t + A t) for all mesh points by extrapolation, and approximated y(t + A t), y(t + A t), 
and T (t + At). We next computed f [T.(t + A t) ] and f (t + At)] from a Lagrangian inter- 

1 8 1 A 

polation on the solidus and liquidus curves. Recognizing that the interface point and tem- 
perature are crucially linked with the values of temperatures and concentrations of the 
adjoining mesh points, we recomputed the values of C and T at these points using a 
difference form of the partial differential equations. 

Based upon the values of the concentrations at the interface and the adjoining mesh 

3C 3C 

S A 

points, we computed the partial derivatives - — and r and determined y(t At) from 

o x o x 

the boundary condition Eq (49f). Substituting in Eq (49g), we determined a new value for 
the interface temperature. When the values of the interface derivative and temperature 
agreed sufficiently with the previous estimates, these values were accepted and y(t + At) 
is determined from y(t + At) = y(t) + | [ y(t) + y(t + At)] At. With y(t + At), Tj(t + At), 



} 

t 


I 

i 


i 
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and the values of concentrations and temperatures at the adjoining points determined, we 
computed f g jiyt + At)J and ^T^t + At)j and solved for the temperature and concen- 
trations at the other mesh points. If y(t + A t) and T^t + At) did not agree sufficiently, 
we averaged them and the iteration repeated. If conversion was not achieved after a 
designated number of iterations, new values for the mesh points are determined. If 
convergence was not achieved after a similar number of iterations, we halved interval 

At and determined values for T(x, t + ^~)* C(x, t + — ) t y(t + ^ ), yit + ~), and 

T.(t + ^~) by interpolation. The process was repeated until convergence was attained 

or the time interval cut down below an acceptable level. In the latter case, the program 

either accepts the values attained or is stopped. An additional device which aided in 

inducing convergence was based on a check of montonicity. The temperature must 

monotonically increase from x - 0 to x = », as must the interface temperature. Thus, 

if values that violated this monotonicity appeared, they were replaced by those found by 

linear interpolation. Similarly, th? concentrations in the liquid will decrease from 

x = y(t) to x = ® . A monotonic behavior, either nonincreasing or nondecreasing in the 

solid phase depending upon the condition upon C (o,t), can alBO be expected. 

s 

The values of the concentrations and temperatures for the mesh points adjoining 
the interface point were also determined using this iteration process. Although con- 
vergence checks were not made on concentration and temperature values it is clear that, 
since these values ultimately involved in the expressions for computing the interface 
derivative and temperature, their failure to converge would have prevented the con- 
vergence of the derivative and interface temperature. 

4 . 2 NUMERICA L APPROXIMATIONS 

In order to determine the values of y/t + At) and T { (t +■ At) from the boundary 

conditions at the interface, one must have a means for approximating partial derivatives 
> u 

of the form r — . where u is a function of x and t. In addition, computation of the con- 
o x 

centrations and temperatures at the adjoining mesh points require approximations for 

5 2 

expressions of the form ~ and — 5 . We have already indicated that y(t + At) is found 
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from y(t + At) = y(t) + \ y(t) + y(t + At) At, sometimes called the modified Euler 
equation. Although this is a lower-order method than the usual Runge-Kutta employed 
by other researchers in similar situations, it is preferable here because y(t) is computed 
by an elaborate iterative scheme. Runge-Kutta involves the use of intermediate values 
of y(t), which would have to be computed in the same complicated manner. Only when 
explicit expressions are available for y(t) can Runge-Kutta be recommended. 

The expressions for the partial derivatives are base upon difference expressions, 
assuming unequal spacing, because the presence of a vary,ng interface point requires 
approximations involving the values at this point. Although economies in computing are 
available if the subdivisions of the mesh are equispaced, they must be sacrificed because 
of the varying interface point. To require that the interface points be at mesh points 
would necessitate a variable time step and a very small space interval in order that 
convergence on y and T. be reached after reasonable number of iterations. In addition, 
the semi-infinite character of the problem, with little change in variables away from the 
origin, suggests that the space intervals increase in size. To accommodate such a 
situation, the key expression for the partial derivative at x = Xj is given by 

uiXj + Ax, t) - ufx^ - Ax,t) 

— , i.e., by a second order approximation. Thus, 


u(x 2 , t) - UfXj, t) 

X - X 

2 1 


is an approximation for at x = 


X + X 
1 2 


Similarly, the 


expressions for — r based upon the values of u(x , t), u(x , t), and u(x_, t) is given by 
K ^ 12 3 

dx 

r u(x 3 , t) - u(x 2 , t) u<x 2 t) - u(x lf tn 
L x 3 " *2 X 2 - x i J ' '* 3 ' * 1>/2 

As a consequence, if the partial at x = is required for given values at x^, x 2> 

/hu\ U(X 2’ l) ‘ U(X 1’ l) X 1 " X 2 a 2 u 

and x„, the formula [ r - ) , = + — - is used, where 

3 \Sx/ x 0 - x, 2.2 

\ / x 2 1 


u(x 2 , t) - UtXj, t) Xj - x 2 
*2 ~1 + 2 


is used, where 
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\«*t" > arwma 


z is computed as above and gives second-order accurate results. An arithmetic 
ax 

function statement of the form D2(X, F, Y,G,Z,H) = «H-G)/(Z-Y)-(G-F)/(Y-X))/(Z-X)/2) 

a 2 u 

is used to compute — using the six quantities x , u(x , t), x , u(x , t), x. and 
dx 112 2 3 

u(x 3 , t). 

a u a u 

We use a similar approach to interpret — — . — - is approximated by 

o * o t 2 

u(s,t+At) - u(x,t) At .. oU a u 

— — *— 1 , evaluated at t + — . Thus, the equation — = q — — is approx- 

Z ot Z 

o x 

mated by 

u t «, 1 - U|x, 1) . j q /V|\ J£l\ 

\ 3X ' t ' Sx f t + At 

2 

Here at t + has been replaced by the average of the values at t and t + At. This 

dX a 2 

is called the Crank-Nicholson method (Ref. 81). Since the evaluation of — ~ at t + At 

ax 

involves the values of u(x, t + At) at three values of x, this equation has three unknowns. 


However, as an iterative method, one can assume the values needed for 


are known 


and use the equation to find an improved value for u(x, t At), in this way the values v*f 
the concentrations and temperatures at the interface-adjoining mesh points are determined. 
This is a reasonable approach since many iterations may be required to compute 
y(t + At) and T { (t + At). 

Evaluation of y(t + A t) is obtained from 


. / bC SC \ 

(w - w] 


The partials are evaluated from the formula for . To find an iterate of T.(t ■* At) 

dx i 
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the equation py y (t + A t) = K — - K 

8 


with partials at t + At, we replace 


3T T -T x -y(t + At) 2 3T. T - T. 

? b ! + J i_T . — 1 b hi 1 

ax Dy y(t+At)-x 2 . 2 5x y x, - y(t + At) 

J «X J-l 


y(t + At) - x. 




, where the subscripts j and j-l indicate the mesh 


point after anH before the interface point, respectively. We can solve for T., excepting 


the T. terms involved 
i 




This gives an iterative formula for 


T (t + At). We see that the values of concentrations and temperatures at the j and j-l 
points play an important role in determining y(t + At) and T.(t + At), and their convergence 
presumes the convergence of the four other values. 

4.3 THE SOLUTION FOR REMAINING MESH POINTS 

When the values of interface derivative and temperature converged, we determined 
the interface point and concentrations of liquid and solid at the interface. The values for 
the concentrations and temperatures must then be computed for the remaining mesh points. 
To determine these values, four systems of equations, one each for concentrations and 
temperatures for each side of the interface point, must be solved. Since the values on the 
infinite side of the interface point change gradually, one need not consider too many mesh 
points. By starting with a small number of points, chosen so that the last few show no 
change from the initial values of C q and T q , the number of points is increased whenever 
the values at the pentultimate point shows some nontrivial change. The new points again 
have values C and T . The values C (o,t) and T (o,t) are determined from the initial 

O O 8 S 

conditions so that they, too, .re already known. The values of C (x ,t) and T (x. , t), 

s j-l s j-l 

where j-l indicates the mesh point before the interface point, and the values of C.(x , t) 

^ J 

and T.(x. t), where j indicates the mesh point after the interface, are also known from 

j 

the iterative process. At each of the interior points, the Crank-Nicholson equation 
involves the values at three mesh points: before, at, and after the point for which the 
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equation is given. Each system of equations is not only linear (since the coefficient of 
2 

^ is assumed constant), but also tridiagonal. A subroutine (TRIST) is used to set up 

ax 

the coefficients of each of the systems and calls upon a standard Grumman subroutine 
(TRIDI) to solve them. 

The same procedure is used when convergence has not been achieved and an 

updated set df values for concentrations and temperatures are desired for continuing the 

2 

i u 

iterations. This need results because the expressions for — - at j -1 or j depend not 

5x 

u:ily on the points subscripted by j-i anu j but also on those whose subscripts are j-2 
and j+1. 

4.4 STARTING PROCEDURES 

The previous discussion assumes that y(t), y(t), T(t;, and all values of Mx,t) 
and T(x,t) at the mesh points are known, and that first approximations are available 
for time t+At. However, all we know at time t = 0 is that y(0) = 0, but we have no 
knowledge of y(0) and T.(0). Indeed, it may even be ture that no solidification begins at 
time t = 0, because the temperature at x = 0 may be above the temperature T at which 
solidification can begin for the given initial concentration C^. In the latter case, until 
the surface temperature reaches T the problem is that of pure heat conduction, and for 
each time interval one solves a single tridiagonal system. After each interval, we 
ascertain whether the temperature at the pentultiinate point has changed nontrivally so 
that points should be added to the mesh. 

However, when the temperature at the surface is below T solidification occurs 

and, hence, a need arises for a starting procedure to determine y(t) and T.(t). An 

appreciation of the problems involved can be gained from the analytic solution for the 

case of constant surface temperature. In this case, as was discussed earlier, y=20Vt^ 

y= -f=, and V interface temperature is constant at T.. Then y(0) is infinite and, hence, 
vt 1 

no value that will give useful results can be assigned to y(0). One way to begin the 

numerical problem is to assume that the analytic solution holds at least for the first time 

interval. This would only be useful for the case of constant surface temperature. In 

addition, the purpose of the analytic solution, i.e. , to check the accuracy of the numerical 
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solution, would lose some of its validity. As an alternative, w~- may use some aspects 

of the analytic solution to assist in starting the solution. One such aspect is that 

• dy 

y(t) = 2y(t) . t. This enables us to estimate y(t) if an approximation of “ is given. 

However, having no value to T. and y at t = &t to begin with, we have no means of 
dy 

computing using boundary equation (49f). As a variant, we may begin with an initial 
estimate of Ay, determine T. by assuming a linear law for temperature from x = 0 to 

1 . . dv 5C ac » 

x - 2Ax, and determine At from 7 ^ = 7^ = (K ' — - K ) ~ ff (T.) - f (T.)l 

At dt s 3x l ?>. Is i £ i J 

From this At, one finds another which is of the form DELT/2 n and less than previous 

at, where DELT is an input to the program. A corresponding value for y(A t), T(Ax, At) 

is estimated linearly. From then on, the iteratives continue in the usual manner, except 

that y is computed as 2y At. Convergence occurs, or else At is cut back. The next 

time step no longer requires a special starting procedure. It is useful to choose the 

original Ay somewhere between | and A of the length of the first space interval A X. 

This removes y far enough from the mesh points to make the differences meaningful. 

4.5 PROBLEMS CAUSED BY INTERFACE APPROACHING MESH POINT 

One of the greatest obstacles in the development of the program occurred when 
the interface point approached a mesh point. The time interval was halved, down to the 
limit, until it failed to converge for the smallest interval allowable. Relaxing the 
minimum time interval resulted in the use of an inordinate amount of computer time, 
nor did it always succeed in breaking through this mesh point barrier. In addition, dis- 
continuities in interface derivative and interface temperature appeared after the mesh 
point was passed. There was a tendency to correct itself, however. Hence the dis- 
continuity due to the definition of the partial derivatives with respect to different points 
was an unavoidable but self-disappearing evil. 

The first solution to the mesh point approach problem was the development of an 
extrapolation procedure whenever the interface point neared the mesh point. For 
example, when the interface traversed 9/10 of the mesh width, the results of the last 
two points were extrapolated linearly to the mesh point. This carried the program 
through the difficulty with only a small error due to linear extrapolation. 
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Another solution to the approach problem was to change those points adjoining 
the interface point. Whenever the interface point came within an input fraction (1-S) 
of the mesh width, the adjoining point on the infinite side was shifted to the next point for 
all iteration considerations and the values of concentration and temperature at the skipped 
point were determined by a linear interpolation from the interface point to the assumed 
adjoining point. On the other hand, the adjoining point on the zero-side of the interface 
was not shifted until the interface point had passed a distance a ratio S of the mesh width. 
For example, if S = .15 and the interface point was .85 Ax, the adjoining point on the 
infinite side was shifted; and when it was 1.15 Ax. the adjoining point on the zero -idc 
was shifted. With this strategy, extrapolation was replaced by linear interpolation for 
the skipped point. Of course, a change in reference points made for discontinuities in 
derivative and interface temperature, but they seemed less severe and were quickly 
corrected as the interface moved away from the critical region around the mesh point. 

4.6 TIME INTERVAL CONTROL 

As indicated in the program, the time interval is subject to variation. The program 
allows for large time intervals when convergence occurs quickly, and smaller intervals 
when convergence is slow. The two factors are the initial DELT and the number of 
iterations NIT. DELT represents the largest allowable time interval, and DELT/128 is 
the smallest. If the number of iterations to convergence is less than or equal to 
NTT/4 then, for the next time interval, a time step twice as large is permitted. If, 
on the other hand, after NIT/2 iterations no convergence has been achieved, then the 
time interval for the current step is halved, if allowable, and maintained for the next 
step. If after NIT iterations no convergence has occurred, then the step size is halved 
again, if allowable, and iteration count restored to NIT/2 + 1. If it is not allowable, the 
program may stop or accept nonconvergent results as set by input. Another situation 
leading to step halving is the occurrance of a negative . Since this is not physically 
realizable, it is assumed that the time step being used is too large, and a smaller one 
substituted for accurate results. 

Another option that may affect time interval control is the desire to print results 
at fixed intervals. Thus, if results are required at an instant of time of t + At exceeds 
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that time, then At is adjusted so that t + At will be exactly as required. Also, for the 
programs where the surface temperature is variable and no solidification takes place 
until temperature T is reached, then the solidification time is used by the program for 
computation. If this is not the desired time for printout, the user will not be aware of 
these results. 

Section 5. COMPUTER PROGRAMS FOR NUMERICAL SOLUTIONS 

Four main programs have been developed in accordance with the ideas described 

above. Two programs solve the constant surface temperature problem, one (B INCR 3) 

using extrapolation to the mesh point and the other (B INCR 2) using the shifting of 

adjoining point method. The other two programs both use the shifting of adjoining points 

to solve for the linear surface temperature (B INCR) and the experimental surface tern- 

» 

perature (B INCR 4). These four main programs use the same subroutines: MOTON to 
preserve monotomicity, and TRIST to set up and solve tridiagonal systems. In all four 
programs, the solidus and liquidus curves are input in the table manner as discussed in 
BISCL. Also, the mesh points are set up in the same way as in BISCL. 


•i 
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The results in Tables 3-1 through 3-3 were obtained using BISCL for the analytic 
solution of the solidification of 12% Sn in Ni, with an instantaneous change of surface tem- 
perature from 1800°K to 1400°K at time t = 0. The Ni-Sn phase diagram is from Hansen 
(1958), while material property data was from Smithells (1968). The coefficient of 
3(t) ® P + at + p t + p t 3 and T.(t) = T. + T t + T t 3 + T.„t 3 are given. The values 

0 14 o 1 10 11 14 lo 

of interface point, interface derivative, and interface temperature are given for times 
t = . 03125 sec to 1 sec, at intervals of . 03125 sec. The values of T(x,t) and C(x,t) are 
given for t = . 03126 sec, . 0625 sec, and 1 sec. 
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TABLE 3-1 COMPUTED SOLIDIFICATION OF 12% Sn IN Ni, 

ACCORDING TO THE ANALYTIC SOLUTION (Sheet 1 of 2) 


0. 540109E-01 

T. 

10 


0. 157333E-04 


TIME-t 


0.312500E-01 


0 , 


-0. 888091E-04 


T. 


il 


-0. 168538E-00 


y(t) 


0. 190948E-01 


TIME-t 


0. 6250000E-01 


0 , 


0. 897718E-04 


i2 


0.134766E-00 


y(t) 


0« 305487E-00 


0.0 

0. 140000E 04 

0.333330E-01 

0. 100000E-01 

0. 149360E 04 

0. 598702E-01 

0. 200000E-01 

0.157956E 04 

0. 1958 08 E 00 

0.400000E-01 

0.169378E 04 

0. 139609E 00 

0. 800000E-01 

0. 178779E 04 

0. 120246E 00 

0.160000E 00 

0. 179999E 04 

0.120000E 00 

0. 320000E 00 

0. 180000E 04 

0.120000E 00 

0. 640000E 00 

0.180000E 04 

0.120000E 00 

0.128000E 01 

0.180000E 04 

0.120000E 00 

0. 256000E 01 

0. 180000E 04 

0.120000E 00 

0.512000E 01 

0.180000E 04 

0.120000E 00 

0.102400E 02 

0. 180000E 04 

0.120000E 00 

0. 204800E 02 

0. 180000E 04 

0.120000E 00 

0.409600E 02 

0. 180000E 04 

0. 120000E 00 


y(t) 


0. 270028E-01 


0, 


-0.425002E-04 


i3 


-0. 559896E-01 


T.(t) 


0. 157332E-Q4 


y(t) 


0. 215984E 00 


X 

T 

C 

0.0 

0. 140000E 04 

0. 333330E-01 

0.100000E-01 

0. 146657E 04 

0. 597497E-01 

0. 200000E-01 

0.153083E 04 

0.598738E-01 

0.400000E-01 

0. 163269E 04 

0. 166143E 00 

0. 800000E-01 

0. 175121E 04 

0.124121E 00 

0.160000E 00 

0. 1799 05 E 04 

0.120001E 00 

0. 320000E 00 

0.18 000 OE 04 

0.120000E 00 

0. 640000E 00 

0. 180000E 04 

0. 120000E 00 

0. 128000E 01 

0. 180000E 04 

0. 120000E 00 

0. 256000E 01 

0.180000E 04 

0.120000E 00 

0. 512000E 01 

0.180000E 04 

0. 120000E 00 

0.102400E 02 

0. 180000E 04 

0.120000E 00 

0. 2048 00E 02 

0. 180000E 04 

0. 120000E 00 

0.409600E 02 

0. 180000E 04 

0. 120000E 00 


T.(t) 


0.157332E 04 





TABLE 3-1 COMPUTED SOLIDIFICATION OF 12% Sn IN Ni, ACCORDING TO THE 

ANALYTICAL SOLUTION (Continued) (Sheet 2 of 21 

TIME-t y(t) y(t) T.(t) 

0. 937500E-01 0.330701E-01 0.176329E 00 0.157331E 04 

0. 125000E 00 0.381845E-01 0.152690E 00 0.157331E 04 

0. 156250E 00 0.426900E-01 0.136557E 00 0.157330E 04 

0. 187500E 00 0.467628E-01 0.124649E 00 0.157330E 04 

0. 218750E 00 0.505080E-01 0.115895E 00 0.157329E 04 

0* 250000E 00 0. 539936E-01 0.107935E 00 0.157329E 04 

0. 281250E 00 0. 572672E-01 0.101757E 00 0.157329E 04 

0. 312500E 00 0. 603633E-01 0.965308E-01 0.157328E 04 

0. 343750E 00 0.633079E-01 0.920349E-01 0.157328E 04 

0. 375000E 00 0.661214E-01 0.881137E-01 0.157328E 04 

G.4G62GGE GO 0.66SI3SE-Q1 u. 84bO‘fctJ,h-ul 0.157328E 04 

0.437500E 00 0. 714163E-01 0.815727E-01 0.157327E04 

0.468750E 00 0.739214E-01 0. 788048E-01 0.157327E 04 

0. 500000E 00 0. 763443E-01 0.763006E-01 0.157327E 04 

0.531250E 00 0.786926E-01 0. 740207E-01 0.157327E 04 

0.562500E 00 0.809726E-01 0. 719334E-01 0.157326E 04 

0. 593750E 00 0.831901E-01 0.700130E-01 0.157326E 04 

0. 625000E 00 0.853499E-01 0.682382E-01 0.15732GE 04 

0.656250E 00 0.874563E-01 0.665914E-01 0.157326E 04 

0.687500E 00 0.895131E-01 0.650579E-01 0.157325E 04 

0. 718750E 00 0.915235E-01 0.636250E-01 0.157325E 04 

0. 750000E 00 0. 934905E-01 0.622823E-01 0.157325E 04 

0. 781250E 00 0.954170E-01 0. 61 '/203E-01 0.157325E 04 

0. 812500E 00 0. 973051E-01 0.598312E-01 0.157325E 04 

0. 843750E 00 0.991571E-01 0.587082E-01 0.157325E 04 

0. 875000E 00 0.100975E 00 0.576451E-01 0.157324E 04 

0. 906250E 00 0.102760E 00 0.566368E-01 0.157324E 04 

0.937500E 00 0.104515E 00 0.556785E-01 0.157324E 04 

0. 968750E 00 0.106241E 00 0.547660E-01 0.157324E 04 

0. 100000E 01 0. 107939E 00 0.538958E-01 0.157324E 04 



0.0 

0. 140000E 04 

0. 333330E-01 

0.100000E-01 

0.141674E 04 

0.471641E-01 

0. 200000E-01 

0. 14 3344 E 04 

0. 557259E-01 

0.400000E-01 

0.146658E 04 

0. 597815E-01 

0. 800000E-01 

0. 153085E 04 

0. 599058E-01 

0. 1G0000E 00 

0. 163269E 04 

0. 166167E 00 

0. 320000E 00 

0.175121E 04 

0.124234E 00 

0. 640000E 00 

0.179905E 04 

0. 120001E 00 

0.128000E 01 

0. 180000E 04 

0. 120000E 00 

0.256000E 01 

0. 180000E 04 

0. 120000E 00 

0.512000E 01 

0. 180000E 04 

0.120000E 00 

0.102400E 02 

0. 180000E 04 

0. 120000E 0 r 

0. 204800E 02 

0.180000E 04 

0.120000E 00 

0.409600E 02 

0. 180000E 04 

0. 120000E 00 
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These results in Table 3-2 and Table 3-3 ar*> to be compared with those obtained 
by the numerical solution, B INCR 3, described above. The times are from . 0625 sec 
to 1 sec in intervals of . 0625 sec. The values of T(x,t) and C(x, t) are given at 
t = . 00097653 sec, . 0625 sec and 1 sec. 

TABLE 3-2 COMPUTED SOLIDIFICATION OF 12% SnIN Ni, 

ACCORDING TO THE NUMERICAL PROGRAM 


TIME, t 


0. 976563E-03 

X 0.00 

T 0. 140000E 04 

C 0. 333333E-01 

TIME, t 
0. 625000E-01 

X 0.0 

T 0.140000E 04 

C 0. 333333E-01 

TIME , t 

0. 127307E 00 
0. 187500E 00 
0. 250000E 00 
0.312500E 00 
0.375000E 00 
0.437500E 00 
0. 500000E 00 
0.562500E 00 
0. 625000E 00 
0.687500E 00 
0.812500E 00 
0.875 000E 00 
0. 937500E 00 
0.100000E 01 

X 0.0 

T 0. 140000E 04 

C 0. 333333E-01 


0. 262523E-02 


0. 13441 2E 01 


0. 100000E-01 

0. 200000E-01 

0.400000E-01 

0.176418E 04 

0.179763E 04 

0.179996E 04 

0.128821E 00 

0.120274E 00 

0.120002E 00 


y(t) 

0. 272345E-01 


y(t) 

0. 23335.1 E 00 


0.100000E-01 

0. 200000E-01 

0.400000E-01 

0. 146771 E 04 

0.153287E 04 

0.1 63831 E 04 

0.553996E-01 

C. 578464E-01 

0.160132E 00 


0.400000E-01 
0.476896E-01 
0.552284E-01 
0. 623012E-01 
0.687872E-01 
0. 746791E-01 
0. 800672E-01 
0.841858E-01 
0. 881141E-01 
0.919827E-01 
0. 996864E-01 
0.103489E 00 
0.107246E 00 
0.110950E 00 


0. 16305 7E 00 
0. 123216E 00 
0. 117407E 00 
0. 108 731 E 00 
0. 989489E-01 
0. 89925 7E-01 
0. 690561E-01 
0. 63G305E-01 
0. 614704E-01 
0. 619815E-01 
0. 611697E-01 
0.605048E-01 
0.597120E-01 
0.588031E-01 


0.153064L 04 

1 0. 800000E-01 

4 0. 180000E 04 

0 0. 120000E 00 


0.157701E 04 

1 0. 800000E-01 

4 0.176099E 04 

0 0.124030E 00 


0.15 768 7E 04 
0. 158008E 04 
0. 157892E 04 
0. 157872E 04 
0. 157869E 04 
0. 157852E 04 
0. 157848E 04 
0.158206E 04 
0. 158280E 04 
0. 158148E 04 
0.158031E 04 
0. 157999E 04 
0.157973E 04 
0. 157957E 04 


0.100000E-01 
0.141392E 04 
0.461461E-01 


X 0.160000E 00 0.320000E 00 

T 0.163849E 04 0.176112E 04 

C 0.158612E 00 0.123817E 00 


0.200000E-01 
0.143380E 04 
0.541801E-01 

0.640000E 00 
0.1 7981 2E 04 
0.120118E 00 


0.400000E-01 
0.146732E 04 
0.572740E-01 

0.128000E 01 
0.179999E 04 
0.120000E 00 


0. 800000E-01 
0. 15324 7E 04 
_0^573063E-0^ 

0. 256000E 01 
0. 180000E 04 
0.120000E 00 
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TABLE 3-3 COMPUTED SOLIDIFICATION OF 12% Sn IN Ni, COMPARISON BETWEEN 
ANALYTIC SOLUTION AND NUMERICAL PROGRAM 


Time, 

Interface 


Velocity of 


Interface 

t, sec 

Position 


Solidification 


Temperature 


Ay(t) 

Any 

y a (t) 

Ay(t) 

kill) 

y» 

A 7*(t) 

ATj(t) 

T*(t) 

0.25 

1.24 x 10“ 3 


2.30% 

9.47 x 10 -3 


8. 78% 

5.63 

0. 56% 

0.50 

3.73 x 10" 3 


4.88% 

-7.24 x 10" 3 

- 

-9.51% 

5.21 

0. 33% 

0.75 

2.36 x 10" 3 


2.53% 

-6.20 x 10~ 4 


-1.00% 

7.60 

0.48% 

1.00 

3.01 x 10' 3 


2. 78% 

5.01 x 10~ 3 


9. 31% 

6.33 


Average 

2.58 x 10" 3 

3.12% 

0.26 x 10~ 3 

1.90% 

6.19 

0. 39% 


The average accuracy of the numerically computed results in this example 
(12% Sn in Ni) is only 3.12% off the analytic values on the interfacial position y(t); 
only 1.90% on the velocity of solidification y(t); and only 0.39% on the interfacial tem- 
perature Tj(t). 

Within the range of study, our numerical program is therefore considered to be 
adequate for some applications to the solidification problem. 
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PART IV 


SUGGESTIONS FOR FUTURE WORK 
Section 1. IMPROVEMENT OF MODEL OF SOLIDIFICATION 

The present analytic model assumes all the material parameters to be constant. 
Such parameters include: 

• Specific heats 

• Conductivities 

• Latent heat of fusion 

• Diffuse-ties 

• Densities. 

In any real situation, these parameters are varying either as a function of concentration., 
or temperature, or both. By refining the program to account for such variations, we 
move toward a more accurate model. Even more significant, however, is the fact that 
we will have a tool for accurately weighing the role and significance of each of the para- 
meters in the solidification process, both in one gravity and in zero gravity. 

The second Improvement in the model Involves introduction of the evaporation 
effects which were so significant in work conducted during the first part of the contract. 
We have already briefly discussed evaporation and the role it plays in solidification. 

Only by accounting for this effect can we clearly define the Initial conditions of the melt 
before solidification in order to ascertain zero gravity effects during solidification. 

Finally, with the help of several fluids specialists in the Research Department at 
Grumman, we have begun to analyze the role of liquid metal fluid dynamics on solidifi- 
cation. Even in zero gravity, there is considerable fluid motion brought about by forces 
such as surface tension, field effects, etc. , possibly resulting, a^ain, from evaporative 
segregation. We want to account for, in some fashion, the effects of such fluid motion 
on solidification. 
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Section 2. EXPERIMENTAL CORRELATION 


Analytical programs are a valuable part of the investigation and, for this reason, 
we want to use the programs at all stages of their development to aid in understanding 
solidification and segregation in zero gravity. Effects of initial and of boundary condi- 
tions, concentration range, etc. , will also be systematically determined. 

In addition, we want to correlate the program with a very real semi-conductor 
material. Certain m-V semiconductor compounds, such as GaAs and GaP, form ideal 
candidates for single crystal growth in space because of their technical Importance, 
stringent demand for perfection, and high cost. GaAs, for example, is a high-mobility, 
temperature-stable material that has been widely used in lasers, light emitters, Gunn 
microwave oscillators, and the Hkti. While these devices are already performing 
numerous unique and Important [unctions, their lives, performances, and reliabilities 
are usually material-limited. That Is, the GaAs starting materials invariably contain 
too many physical and chemical defects in the form of point, line, surface, massive, 
or spiral types. The origin and formative conditions of these defects have not been 
fully explained. A clear understanding of these defects is, therefore, an urgent require- 
ment for material improvement. 

Today, in manufacture, exact control of a compound's stoichiometric composition 
is required, but it may have a phase diagram which is distinctly nonlinear. In addition, 
elemental As and P are so evaporative relative to Ga as to make the required composi- 
tional control extremely difficult and nonreproducible. Yet, most conventional studies 
consider neither the differential evaporative behavior of the constituents nor the curva- 
tures of the relevent phase boundaries. For this reason, w<* have placed major emphasis 
precisely on these two important areas. 

GaAs is, therefore, an ideal candidate material; it is susceptible to itrn. • vement 
through use of analytical techniques. Our programs can begin to explain the formation 
of defects In such a material and, therefore, be used to ascertain the effectiveness of 
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the space environment in improving the material. Finally, the same programs will 
allow us to devise improved GaAf ' pace growth procedures. 

Continuing a relationship wnich \ 'as initiated on the current contract, Grumman 
will work in cooperation with :*;9 State University of New York at Stony Brook in an 
effort to correlate our analytical efforts with experimental results. Single crystals 
will be grown at SUNY-SB under highly controlled conditions and carefully prepared 
statistical plans. These materials will then be analyzed and results used to confirm 
findings obtained on the computer. 
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ABSTRACT 


The phenomenon of evaporative segregation in binary alloys ha 
been investigated through a study of some experimental evaporation 
data relating to the Ni-Cr and Ni-Fe systems. In nornal evap- 
oration it is assumed that a) the evaporating alloy is always 
homogeneous, b) the vapor is instantly removed, and c) the alloy 
follows Raoult's law. The solutions of the evaporation equations 
for the two most important cases are presented and experimental 
data are analyzed with these equations. The difference between ob 
served and calculated values of evaporation constants lies within 
one order of magnitude. This is surprising because of the major 
assumptions stated above. Experimental results have shown that 
the evaporation time t and final solute concentration m are 
logarithmically related, further supporting our evaporation equa- 
tions. It is further shown that neglecting the nonlogarithmic 
term in these evaporation equations may introduce considerable 
errors in the analysis. 
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DISCUSSION 


The phenomenon of evaporative segregation in binary alloys has 
been quantitatively treated by Li (Refs. 1,2). By evaporative 
segregation, we mean solute enrichment or depletion in the evapo* 
rating alloy as a result of surface evaporation. 

The formalism of evaporative segregation in alloys works 
surprisingly well for the several systems so far studied, 
notwithstanding the several major assumptions involved. In the 
present memorandum, we shall study only some experimental evapora- 
tion data relating to che Ni“Cr and Ni*Fe systems, as reported 
by N. D. Obradovic et al.(Ref. 3). Specifically, we shall compute 
the degree and rate of solute enrichment or depletion in the evapo* 
rating source at different temperatures and evaporation times. 

In the treatment of the simple "normal" evaporation, it is 
assumed that a) the evaporating alloy is always homogeneous, 
b) the vapor is instantly removed, and c) the alloy follows Raoult's 
law. Such conditions exist or are nearly approached in an induction 
melt in vacuum or liquid in space. The solutions of the relevant 
differential equation for the following two important cases are 
(Ref. 2): 

Case I: The solute is much 

i.e., U » V. Then the time to 
is 

[ (1 - ra )m 

h * V" 


more evaporative than solvent, 
reach a final concentration m 


m - m 


(1 - m Q ) (1 - m) 


( 1 ) 
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where 


G 1 " " N o (1 ’ m 0 )/ AU < sec > 

2 

U ■ evaporation rate of the solute In mol/cm /sec 

2 

V * evaporation rate of the solvent in mol /cm ‘/sec 

m ■ Initial molar fraction of solute 
o 

m ■ final molar fraction of solute after time 
t^ (in sec) 

N * total initial number of moles of both solvent 
J and solute 

2 

A • evaporating area in cm 


For nearly equal m and m (i.e., m ^ m ) , or for small 

o o 

m and m Q (i.e., m ~ n> 0 ~ 0) > Eq. (1) reduces to the following 
logarithmic expression: 


t^ - G^ In 


(1 


m )m 

„qL. 


m (1 “ m) 
o 


(la) 


Case II: The solvent is much more evaporative than solute, 

i.e. , V » U, then. 


t 


2 


f (1 * m )ra m - m I 
G ? ln 77 8““ + a 

2| (1 - m ) m 0 mni o j 


( 2 ) 


where Gj ■ m o N o /AV. Again, for nearly equal m and m Q (i.e., 
m ~ n» 0 ) , Eq. (2) reduces to a logarithmic expression very similar 
to (la) : 


t 


2 


ln 


(1 - m Q )m 
m Q (l - m) 


(2a) 
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G's are the time constants of evaporation. For the case 
when U » V, solute depletion takes place, i.e., m < m Q and, 
therefore, the logarithmic term in the parentheses is negative 
(the second term in the parentheses is often negligible by comparl 
son with the logarithmic term). Since the evaporation time, t, 
must be positive, is negative, is always positive. 

2 

The evaporation rate in mol /cm /sec for solute and solvent 
is given by (Ref. 4) 


u - k( 10 A “ B “ /T ) X v ATt 

(3) 

, A *B /T 


V-k(iO V V )/ / HyT 

(4) 


where 

K • 5.833 x 10 5 , a constant for metals 

M .M v * molecular weights of solute and solvent, 
u respectively 

T ■ evaporating temperature in °K 
A u ,B u ■ evaporating constants for ths solute 
A v ,B v * evaporating constants for the solvent. 

Obradovic et al. (Ref. 3) give the final solute concentrations 
for Ni-Fe and Ni"Cr alloys after various evaporating times at 
1600°C under various ambient pressures, as shown in Tables 1 and 2. 
These data have been analyzed by means of Eqs. (1) and (2), above. 
The solute and solvent evaporating constants, A u , B u , A , and B v , 
are obtained from Ref. 4. 
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Relevant data for 

a) 60X Fe-40I Ni 

U M . « 6.386 x 
Ni 

V pe - 1.081 x 

b) 801 Ni-20* Cr 

U Cr - 3.'w9 x 

V Ni « 6.386 x 

It can be seen from Table 1 that for the Fe-Ni alloy, the values 
of 62 *s are constant within one order of magnitude. The difference 
between observed and calculated values of G2 Is also within one 
order of magnitude. This is very surprising considering the fact 
that elemental evaporation rates generally differ by many orders of 
magnitude. Also, the above analysis is based on the assumptions of 
a) the evaporating alloy is always homogeneous; b) ideal solution 
behavior (and the solutions here are highly concentrated!); c) the 
vapor is instantly removed. The experimental data of Obradovic et al. 
(Ref. 3) on the Fe”Ni system shows no systematic dependence of 
evaporation kinetics on pressure as one would expect, and suggests 
large experimental errors, errors in measurements, or both. 

Table 2 also shows the calculated and observed values of G^'s 
for the Ni-Cr alloy. Here, also, the values of G^’s are almost 
constant, and the differences between the observed and calculated 
values of G^’s are smaii. The data on the Ni-Cr system show 
some positive correlation between ambient pressure and the weight 
percent of chromium in the melt after 30 minutes. 


computation of the G's are given below: 
alloy: 

10* 6 mol/cm 2 /se;; N - 0.7478 mol; T * 1600' C; 

o 

10 ^ mol/cm 2 /sec; A * 10.0 cm 2 ; Case II: U « V. 
alloy: 

10 ^ mol/cm 2 /sec; N q = 0.7582 mol; T * 1600 °C; 

10 mol/cm /sec; A * 10.0 cm ; Case I: U » V. 
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The values of computed G's are constant within one order of 

magnitude. In particular, for the Fe-Ni system, at ambient pres* 

3 

sures of 1, 30, 50, 200 (x 10 ) torr, these same values are ac- 
curate to within 25 percent. Therefore, in some cases it is de- 
sirable to know these values to within 20 percent of their real 
values. 

In Tables 1 and 2 the time constants G's have been calculated 
in two ways. First, G's were calculated by using only the loga- 
rithmic term. Next, the complete equation (1) or (2) involving 
both the logarithmic and nonlogarithmic terms was used. It is in- 
teresting to note that the ratio of G's calculated in two cases 
is constant and has an average value of 1.67 0.01 for Fe-Ni 

alloy and 1.26 t 0.005 for the Ni-Cr alloy. Mathematically, 
this means that for the case U » V (in Ni-Cr alloy) 


In 

(1 - m Q )m 

1 m - m 

1 . . o 

m (1 - m) 
Q. ' 

(1 ‘ m 0 )(l “ m ) 

ln| 

(1 - m Q )m 

m (1 - m) I 

l o v ' J 


(5) 


or 


m 


m 


(1 - m Q ) (1 - m) 


In 


(1 - m o )m 


m Q (l - m) 


0.26 , 


( 6 ) 
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and for the case V » U, i.e., in Fe-Ni alloy 


m - m 


m m 


|d - tn)m 
In 2— 


jm o (l - in) 


0.67 


(7) 


Experimental results have shown that the evaporation time t 

and the final solute concentration m are logarithmically related 

(Refs. 3,5), supporting our Eqs. (1) and (2). However, these same 

equations indicate that such relations are true only for vastly 

different solvent and solute evaporating rates, and, in addition, 

for small initial and final concentrations, m and m. These 

o 

same equations also allow one to compute, e.g., the effective evap- 
orating area. A, the solute or solvent evaporating rates U or 
V, or the evaporating temperature T. Further, neglecting the 
non logarithmic second terms in Eqs. (1) and (2) may introduce 
errors in T, A, U, or V by 26 percent in the 80% Ni-207. Cr 
case and 67 percent in the 60% Fe-40% Ni case. 

Ttoe equations of normal evaporation can be made more realistic 
by introducing the concept of solvent-solute interaction. If one 
considers only first -o-der correction, then it can be shown that 
the ideal evaporation i *s of solutes and solvents will be changed 
by a new factor related to the activity coefficienc of the element 
in question. In other words. 


U' 


V 


V' 



IB) 
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where y u , y v are activity coefficients of soluce and solvent, 
respectively, assumed constant in the vicinity of m Q ; and U*, V' 
are effective evaporating rates of solutes and solvents. Using 
the activity data of Fe-Ni alloy (Ref. 6) at 1600C, we have 
recalculated the values of Gj* Now 


where 


C 2* 


m n 
o o 

AV ; 



(9) 


7 V ■ Yp g “ 0.977 for 60 mol % Fe alloy. 

G^ (cal) ■ 0.4471 x 10* 4 sec. 

On applying similar first-order correction for Ni~Cr system 
(Ref. 7), the revised value of G^ becomes, 

Gj (cal) » 0.2138 x 10* 4 sec . 

These corrections, while not very large, do make a little better 
agreement between calculated and experimentally observed G values. 

In conclusion, the simple normal evaporation approach seems 
to give fairly accurate and useful predictions, even for the highly 
concentrated Ni - Cr and Fe-Ni alloys. 
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ABSTRACT 


The kinetics of evaporation in binary alloys have been quanti- 
tatively treated. The formalism so developed works surprisingly 
well for several systems so far studied, notwithstanding the sev- 
eral major assumptions involved. In this memorandum, we have 
studied the kinetics of purification of beryllium through evapora- 
tion data actually acquired during vacuum induction melting. This 
study shows that normal evaporation equations are generally valid 
and useful for understanding the kinetics of beryllium purification. 
The normal evaporation analysis has been extended to cover cases 
of limited liquid diffusion. It has been shown that under steady- 
state evaporation, the solute concentration near the surface may 
be up to six orders of magnitude different from the bulk concentra- 
tion. Corrections for limited liquid diffusion are definitely 
needed for the highly evaporative solute elements, such as Zn, Mg, 
and Na, for which the computed evaporation times are improved by 
five orders of magnitude. The commonly observed logarithmic rela- 
tion between evaporation time and final concentration further sup- 
ports the validity of our normal evaporation equations. 
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INTRODUCTION 

The kinetics of evaporation in binary alloys have been quanti- 
tatively treated by C. H. Li (Refs. 1,2). The formalism so devel- 
oped works surprisingly well for the several systems so far studied, 
notwithstanding the several major assumptions involved. Specifi- 
cally, in this normal evaporation approach, it is assumed that 

a) the evaporating alloy is always homogeneous in composition, 

b) the alloy follows Raoult's law, and c) the vapor is instantly 
removed. In this work, we have studied the kinetics of puri- 
fication of beryllium through evaporation data actually acquired 
during vacuum induction melting by R. F. Bunshah and R. S. Juntz 
(Ref. 3). 
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EXPERIMENTAL DATA 

The details of the vacuum induction melting procedure are as 

follows. Beryllium (SR grade flakes) in the form of a cylindrical 

compact was cruciLle-free induction melted under an ambient pres- 
*6 

sure of 10 torr. The cylindrical compact was 1-3/4 inches in 
diameter and 4 inches high, and weighed about 100-130 g. Since 
the actual temperature of the melt was not known, for the purpose 
of our computations, the temperature was assumed, for the sake of 
comparison, to be 1500°C, a temperature commonly used in Be 
melting (Ref. 3). 
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EQUATION OF NORMAL EVAPORATION 


The solutions of the kinetic equations fo* - two important cases 
are given in papers by Li and Mukherjee et al. (Refs. 2,4). In 
Case I, the solute evaporating rate, U, is much greater than the 
solvent evaporating rate, V, i.e., U » V. The equation re- 
lating the evaporating time t^ to final concentration, m, is, 
for this case, 


Gi 


In 


(1 ~ m p ) tn 
ra o (l - m) 


m * m 


(1 - 


m o )(l * m) 


( 1 ) 


and in Case II when U « V, the evaporation time is given as 
follows: 


where 


In 


(1 - m )m 
o 

m (1 " m) 
o ' 


m 


il 


mm 


( 2 ) 


Cj_ - - N 0 (l - m Q )/AU 
G 0 ■ + N m /AV, and 

4 O O 

m,m Q ■ mole fractions of solute at t » t and t - 0, 
respectively. 


G’s are the time constants of evaporation. For Case I when 
U » V, solute depletion takes place (or m < m Q ) and, therefore, 
the logarithmic term in the parentheses is negative (the second 
term in the parentheses is, for small values of m and m Q , neg- 
ligible by comparison with the logarithmic term) . Since the evap- 
orating time t L must be positive, Gj^ must be negative. G 2 > on 
the other hand, is always positive. 
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In Table 1 the times tj,, to reach the final concentrations 

m from initial concentration m Q (through evaporation from area A 

and with an initial total of N moles of both solute and solvent) 

o 

for each element in the melt have been compiled using Eqs. (1) and 
(2) above. In the table, the ideal evaporation rates at temperature 
T°K of each element have been computed using the formula given in 
Ref. 5, viz . , 

for the solute: 

A -B /T 

U » K 10 u u / v/lTT 

and ( 3 ) 

for the solvent: 

A "B / T 

v - k 10 v v / /-fyr , 

where K • 5.833 x 10 ^ for metals; A , B , A , and B are 
tlemental evaporation constants given in Ref. 5; and the M's are 
molecular weights. 
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TYPES OF IMPURITIES 

Three types of impurities can be distinguished: 

Type I. - The impurities that have much higher evaporating 
rates U than the solvent Be. For these impurities, such as Zn, 

Na, and Mg, evaporation £q. (1) should be used. 

Type II. ~ The impurities that have much lower evaporating 
rates U than the solvent Be. For these impurities, such as Fe, 

Cr, Ni, Cu, Ti, and Si, Eq. (2) should be used. 

Type III. - The impurity that evaporates at about the same 
rate as the solvent Be, i.e., U ~ V. Only a single element, 

Al, is in this group. For this impurity, a different evaporation 

equation in infinite series form should be used (Ref. 2). 
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PRELIMINARY DATA ANALYSIS 


fi 

-» s 



* 

T 


i r -■ 

‘ t ‘ 

\ i 

* * 



Since the exact temperature of the melt was not known, the 

exact evaporation rates could not be computed. However, since all 

the solute and solvent elements were evaporating from the same 

liquid-gas interface of a fixed area A for the same length of tirre, 

we can compute, from the measured m 's and m's, the values of 

o 

P, which is the product of the evaporating time and evaporation 
rate for a given solute with the help of Eqs. (1) and (2). The 
value of P should be constant for all the solutes. 


ideal 


'N (1 - m ) 
o o' 


AU. 


In 


(1 - m Q )m 
m Q (l - m) 


m 


m 


(1 “ m) (1 - m o ) j 


(4a) 


or 


^ ^ “N (1 “ 

t l U - t ldeal U 2— 

C l U i C 1 U i A 


V 


(1 - I r ) 

In — y . — ~ + 
m Q (l -m) 


m 


m 


(1 " m) (1 “ m o ) | 


* P 


(4) 


where t^ is the estimated evaporation time computed separately 
for each element at the assumed temperature of 1500 C. Similarly, 
from Eq. (2) and the ideal evaporation rate of the solvent V. at 
the same temperature, we have: 


ideal 

t 2 


^ideal 


P . 


In Table 1 are listed the value of P's for each element. P 
is defined, as shown, such that 


P ■ t. x U. when V. < U. 
ii i i 


P * t^ x when . 
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It can be seen from Table 1 that P is almost constant for the low 
vapor pressure (Type II) impurities except for copper. 

Because of our three assumptions of normal evaporation, and 
because of the assumed evaporating temperature of 1500 C, the 
computed values of t^ vary from element to element, even though 
they should be identically the same. The variations in t^, how- 
ever, are not very large so that the equations presented here are 
still useful. 

The high vapor pressure (Type I) impurities (Mg, Na, Zn) 
a 1 so have uniform P values, which are roughly one order of magni- 
tude higher than those for the low vapor pressure (Type II) ele- 
ments. The negative value of P for aluminum shows that aluminum, 
cannot be treated by Eqs. (1) and (2) because for this element 
U ~ V and a different evaporation equation in infinite series form 
should be used. 
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STEADY-STATE EVAPORATION 



To explain the large variation in P values among different 
types of impurity elements, one should consider limited liquid dif- 
fusion. The concept of effective evaporation rates, U e and V £ 
upon reaching a steady-state condition is now introduced. 

P = t V = constant (la) 

e e 

or 

c.V. = P * constant * t V , (2a) 

11 e e 

where the subscrip.vs i refer to ideal or normal evaporation val- 
ues, while the subscripts e refer to effective or nonnormal val- 
ues. 

In the following paragraphs, the concept of U e will be de- 
veloped. 

In normal evaporation, we assume that the evaporating alloy 

is always homogeneous, i.e., the liquid diffusion constant D., = *>. 

But because of the large (orders of magnitude) differences in 

evaporation rates of elements, the surface concentration of the 

solute m in the melt is not the same as the bulk concentration 
s 

m Q due to the limited liquid diffusion. Figure 1 depicts the 
situation (at t = 0) (i.e., initial), and after reaching steady 

states, for the case U » V. 

For U » V, the evaporating surface region soon becomes de- 
pleted of solute, and the surface concentration of the solute, m g , 
will be different from the bulk concentration, m Q , specifically, 
m g « m Q . After the initial transients, a steady-state for the 
alloy system and evaporating conditions will be reached. The steady 






Fig. la Initial State, Fig. lb Steady State, 

t - o t » t s 


state concentration profile is shown in Fig. lb. The equation of 
the profile, i.e., m versus x, can be derived according to a 
procedure used by Tiller et al. (Ref. 6) in their study of freez- 
ing with limited liquid diffusion, i.e., 

m - m 0 1 1 + — “ ex p',i & x ) t • ,5 > 

o e 

where the receding rate R of the evaporating surface can be com- 
puted from U^, V^, m g , M u , M v> and material densities. 

In steady-state, 


and 


m U. 
s i 


m U 
o e 


fl - m a )V l - ft - n, 0 )U e , 


or 


Vi 


(1 - m s ) Vi 



(6a) 


(6b) 


(6c) 
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Solving for m g , 


m 

s 


m V. 
.. A..A. 


(1 - m )U. + ra V. 
v o' 1 o i 


(7) 


Since 


m U. = m U 
si o e 


( 8 ) 


where U g is the effective evaporation rate of the solute due to 

the change in surface concentration. Substituting the m from 

s 

Eq. (7) in Eq. (8) gives 


m U.V. 

e m u i (1 - m )U. + m V. 
o o' 1 o 1 

Similarly, for the case V » U, 


(1 - m )V. = (1 - m ) V 
s i o 


e 


(9) 


( 10 ) 


or 


V 


1 - m V.U. 

§. y c: L_i _ IJ 

1 - m i (1 - m )U. + m V. e 
o o' 1 o 1 


(ID 


Thus, it can be seen mathematically that V « U . This means 

j e e 

that in both cases, U » V, or V « U, the effective value of 
U or V is the same. This equality means that in the steadysta te , 
the surface concentration m g does not change with time, as ex~ 
pec ted. 

Now consider Eq. (7) ii some detail. Several cases can be 
distinguished; these are: 
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X. — U » V and m ~ 0 


m 

s 


(1 


m V. 
o 1 


- m )U 
o' 


i 


m V . 
U. 

l 


or 


(Ha) 


m 

s 


m 

o 



II. — V » U 
a) 


U. 

m V » U. or « m 

Oil V . o 


(lib) 


then 


b) 


then 


m = 
s 


1 


m V. « U. or 

o 1 1 



» m 

o 


m 

s 


m V. 
o 1 

(1 - m o )U 1 


m V. 
-Q- ± 


(11c) 


when 


tn ~ 0 . 
o ~ 
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c) 


mV. U. 
o 1 ” 1 


or 


u i 

v:* % 


then 



(lid) 


when 


m ^ 0 . 
o — 


The effective evaporation rate of the solutes, U g , can be 
computed by Eq. (9). The effective evaporation rates have been 
variably raised or lowered to equal the ideal evaporation rate of 
the solvent (Be) . In the case of high vapor pressure impurities 
like Mg, Zn, Na, this means a lowering of evaporating rates by 
about five orders of magnitude. 

The surface concentration, m , and the ratio of surface to 

s 

bulk concentration have been computed in Table 1 by Eq. (7). These 
ratios must also equal V^/lb . Now we have 


The effective time t , to reach the final concentrations m 

e 

have been computed for all the impurities in Table 1; t g should 
be a constant. It can be seen that for the low vapor pressure im- 
purities, t g is almost constant. Again, the t fi for Cu is one 
order of magnitude less than others in the low vapor pressure im- 
purity group. Among the high vapor pressure impurities, the t 
is almost constant. It can be seen that t for the high vapor 
pressure impurities is one order of magnitude larger than that of 
low vapor pressure impurities. 
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NONIDEAL SOLUTIONS 


The normal evaporation equations can be further modified by 
introducing the concept of solvent “solute interactions. If one 
considers deviation from the Raouit’ s law, it can be shown that the 
evaporation rates of solvent and solute wJl be changed by a new 
factor, the activity coefficient. For very dilute solutions, 
which we are considering, evaporation rates for solute and solvent 
can be given by U /=! y U and V' * V , respectively, where 
is the activity coefficient of the solute (Henry constant, in this 
case) . 

In general, such consideration should not affect the estimated 
time, t g , for the type I impurities (Fe, Ni, Cr, Ti, etc.), ex- 
cept in the case of copper. From the available phase diagram, it 
is evident that Be and Cu atoms have significant interactions 
with each other and form many compounds (Refs. 7,8). Thus the 
excess free energy for the system should be negative, resulting in 
a lower value of y_ compared to the ideal solution case 
\y^ u * l.Oj. In such a situation, is not very much higher 

than U' , and this leads to a greater effective time, t 

LU 6 « vll 

to reach the final concentration. The exact value of v is 
difficult to evaluate because of the lack of data in this composi- 
tion range. 

In the cases of Z n and Mg, the solubility in beryllium is 
negligible (Ref. 8) and so > Zn and should be greater than 

one. This would lower the estimated time, t e< meaning a better 
agreement with the other elements. 

It is important to note that exact calculation for the devia- 
tion from the ideal solution case is difficult because of the lack 
of data, and such estimations are purely qualitative. 
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1. Normal evaporation Eqs. (1) and (2) are generally valid 
and useful for the study of beryllium purification. This agrees 
with a previous study on the results of Ni~Fe and Ni”Cr sys- 
tems (Ref. 4). 

2. The computed values of P * t^U^ or t^V^ are essen- 
tially constant for the same type of solute elements, although dif- 
fering by one order of magnitude between different types of ele- 
ments . 

3. The normal evaporation analysis has been extended to cover 

cases of limited liquid diffusion. It is shown that under steady- 

state conditions, the solute concentration near the surface, i.e., 

m , may be up to five orders of magnitude different from the bulk 
s 

concentration, m . 

o 

4. Corrections for limited liquid diffusion are definitely 
needed for the highly evaporative solute elements Mg, Zn, and 
Na, for which the computed evaporation times, t e> are improved 
by three or four orders of magnitude. 

5. After these corrections, the computed times are comparable, 
even between different types of solute elements. 

6. The agreement among the computed times can be further im- 
proved by considering the deviation from the ideal solution condi- 
tion. Because of lack of data, it is difficult to give any quanti- 
tative calculation. 

7. Estimations of solute evaporation behavior are now well 
within one order of magnitude. 



8. The commonly observed logarithmic relation between evapo - 
ration time t and final concentration m further supports the 
validity of our normal evaporation equations. Such relations are, 
however, true only for widely different solute and solvent evapo - 
rating rates and for dilute solutions, as can be seen from Eqs. (1) 
and (2). 
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APPENDIX C 


MOMENTUM GAIN BY EVAPORATING SURFACES 


by 

James T. Yen 

Grumman Aerospace Corporation 
Research Department 
Bethpage, New York 11714 


In this appendix, we evaluate the rate of momentum acquired by a liquid surface 
when its molecules are evaporating into a vacuum. 


First, let us evaluate the evaporating rate and evaporating pressure by considering 
the equilibrium case. When the surface is at equilibrium with its surroundings, the rate 
of molecules leaving the surface by evaporation will equilibrate the rate of condensat-'on, 
resulting in zero gain o? less by the surface. These rates will be called the equilibrium 
rates of evaporation and condensation, and an evaluation of one will yield the value for 
the other. The condensing rate is readily given by 


* 



B. - 

1 

(2n RT) 2 


( 1 ) 


At equilibrium, the rate of molecules impinging on a unit surface area is J nv. 
Here, n is the number density, o the mass density, T the absolute temperature, and R 
the gas constant of the molecules in equilibrium with the liquid. The equilibrium 
evaporating rate also given by Eq (1); this is the maximum rate at which molecules 
can escape from a liquid surface wMcn is kept at a temperature T. Thus, the evapor- 
ating rate into a vacuum from this liquid surface is also given by Eq (!'. 


The preseure p in Eq tl) is the vapor pressure at saturation. Since pressure is 
a scalar quantity, it represents the sum ihe pressures acting rn the sur f ace by the 
evaporating and the condensing molecules. At equilibrium, thes»i two contributions 
equilibrate with each other. Hence, the equilibrium evaporating pressure is equal to 


C-\ 



half of the saturation vapor pressure. When the surface is surrounded by a perfect 
vacuum, the condensing pressure vanishes and the pressure due to the evaporating 
molecules is just given by the equilibrium evaporating pressure, or half of the 
saturation vapor pressure. Consequently, we can expect that the rate of momentum 
acquired by a liquid surface when its molecules are evaporating into a vacuum is 
equal to half of the saturation vapor pressure. This is verified analytically below. 


The molecules evaporating from a liquid surface kept at temperature T have a 

~ 2 -v^/2RT 

Maxwellian distribution function n(2nRT) e . Each beam of molecules 

emanating with a velocity v from a unit surface area and at an angle 9 inclined from 

the normal to the area carries a flux rate of nvcose and a momentum rate of 
2 2 2 

mnv co s 9 = ov cos 6. Hence the total momentum rate in question is given by 


o(2ttRT) 


W 


2 2 « 
v cos 


_v 
2RT 2 


v sinddvde. Integrating v from 0 to » , 9 from 0 _o 


y- , and ® frcm 0 to 2n readily yields | oRT = |p and completes the verification. If 
the surrounding vacuum is imperfect, then the rate is given by 

|or o with 


a can be determined either from the given degree of imperfection or experimentally. 
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ABSTRACT 


Perfect liquid mixing but no solid diffusion is assumed in 
normal freezing. In addition, the molar compositions of the freez- 
ing solid and remaining liquid, respectively, follow the solidus 
and liquidus curves of the constitutional diagram. For the linear 
case, in which both the liquidus and solidus are perfectly straight 
lines, the normal freezing equation giving the fraction solidified 
at each melt temperature and the solute concentration profile in 
the frozen solid was determined as early as 1902, and has since 
been repeatedly published. Corresponding equations for quadratic, 
cubic, or higher-degree liquidus and solidus lines have also been 
obtained. This paper gives the equation of normal freezing for 
ideal ternary liquid solutions solidified into ideal solid solu- 
tions of the pseudobinary type. Sample computations with the use 
of this new equation were made and are given for the Ga-Al-As 
system. 
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INTRODUCTION 


Perfect liquid mixing but no solid diffusion is assumed in 
normal freezing. In addition, the molar compositions of the 
freezing solid and remaining liquid, respectively, follow the 
solidus and liquidus curves of the constitutional diagram. 

For the linear case, in which both the liquidus and solidus 
are perfectly straight lines, the normal freezing equation giving 
the fraction solidified at each melt temperature and the solute 
concentration profile in the frozen solid was determined as early 
as 1902 (Ref. 1), and has since been repeatedly published (Refs. 2-9). 
Corresponding equations for quadratic, cubic, or higher-degree 
liquidus and solidus lines have also been obtained (Refs. 10-12). 

A recent publication presents the normal freezing equation for 
ideal dilute solutions or alloys in which the liquidus and solidus 
are exponential functions of the reciprocal absolute melt tempera- 
ture (Ref. 13). The logarithm of the mole fraction segregation 

coefficient k for the dilute component is, therefore, also a 
m 

linear function of the reciprocal of the absolute temperature. 

Further, in such ideal solutions, the ratio of partial pressure 
of each component, p^, to that of the pure component, p°, 
equals the atomic fraction x^, i.e., x A ■ p A /p°. This last 

equation permits the computation of the solidification behavior 
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of such solutions or alloys and the resultant solid concentration 
profiles from basic thermodynamic quantities without the use of 
phase diagrams. 
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THE PSEUDOBINARY SYSTEM 


This memorandum gives the equation of normal freezing for 
ideal ternary liquid solutions of A, B, and C atoms (Fig. 1) 
solidified into ideal solid solutions of AB-AC compounds with 
A atoms on one sublattice and B or C atoms on the other. 
Such systems can be treated as pseudobinary systems. Account is 
taken of the entropy of mixing of the three atomic species to 
form the ideal liquid solutions along the pseudobinary. 



Fig. 1 Constitutional Diagram for a Ternary 
System of the Pseudobinary Type 
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For the compound AB, we have 


AB(s) - A(l) + B(l) . 


The equilibrium constant for the reaction is 


K - 



( 1 ) 


( 2 ) 


and the standard free energy AF°g and the free energy of forma- 
tion AF^ are related to K at temperature T by 

nD 

- tf AB + AF AB * RT ln K ' <3 > 

For the compound AB, the free energy of formation, AF^g, 
is related to the free energy of melting, and the heat of 

melting, AH™g, at the melting point T™ fi by 

* AH AB (1 - T/I AB> * _aF AB ' RT ln 4 • < 4) 


It has been assumed that the difference in heat capacity be- 
tween the liquid and solid is negligible. The -RT • ln 4 term 
is the ideal entropy of mixing to form the stoichiometric liquid 
phase . Hence , 

x 1 • x 1 

-RT ln 4 * A B - AH^g (1 - T/T® B > . (5) 

X AB 

Similarly, for the compound AC, we have 
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RT In 4 • 


( 6 ) 


XjL T^ - <c < l - t ' t ac> * 

X AC 

where AH^, is Che heat of formation of the compound AC, assumed 
constant; and T™ c is the melting point of AC. 

For the tern^y liquid phase of the pseudobinary, 

X A ' i ’ 


1 . 1 l 

C B + X C " * 


(7) 


Hence , 


and 


-RT In 4 


A • x 1 /x S 
? B' AB 




t/t Ib> 


( 8 ) 


-RT In 4 • i • (i - xj)/(l - x® B ) - AH® C (1 - T/T^ c ) . (9) 

Equations (8) and (9) are essentially the same as those previously 
given (Refs. 14,15). From these equations, 


n 1 S 

2 ' *8 * X AB ' eXp 


AH 


m 


.m 


tM 1/T * 1/t ac 


S ; 

x AD * exp(a + b/T) ■ x 

Ao i 


AB 


U 


( 10 ) 


i * 2 • x i - 


(! - x AB ) • exp 


AH 


m 

AC 

R 


1/T - 1/T 


m 1 
AC J 


( 11 ) 


(1 - x Afi ) • exp(f + g/T) - (1 - x AB ) • V 


where 
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b " -AH AB /R » a "‘ b/T AB * g "' AH AC /R • and f *" g/T AC (12) 

are constants for a given ternary system of the type shown in 
Fig. 1, while the temperature-dependent variables 

U - exp(a -f b/T) and V ■ exp(f + g/T) . (13) 

Equations (10) and (11) yield the solidus equation for the 
compound AB 

1 " X AB ' U “ ( X * X AB^ ‘ V * ( 14 > 

Hence , 

X AB " (1 " V)/(U ' V) * (15) 

s 

The value of x Ag can be as small as zero when, according to 
Eqs. (15) and (13), V - 1 » exp(f + g/T), i.e., when T - -g/f - 
T®£ and x^ • x Ag /2 ■ 0 [Eq. (11)), as expected. The value of 

g 

x AB can also be as large as 1 when U * 1 ■ exp(a + b/T) in 
Ao 

accordance with Eqs. (15) and (13), i.e., when T ■ -b/a » T Ag 
X s 

and x n ■ x a _/2 - \ [Eq. (10)], also as expected. 

d Ad 

Notice that there are no AB or AC compounds in the melt, 
while in the solid there are only compounds AB and/or AC, but 
no individual elements A, B, and C. Hence, the thermodynamic 
molar concentrations x®g, x® c , x^, x^, x£. 

Considering the ternary liquid or solid solutions as a whole, 
and, without regard to crystallogiaphic lattice arrangements, we 
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can define the chemical molar concentrations m^, m^, and m^. 


Evidently, 


s i i 
m A ' m A ■ 4 • 


”B + m c‘ n B +m c‘ i ’ 


m B “ X B * 3nd 


( 16 ) 


s i 
m B ' 4 X 


AB ' 
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i v- •* 


NORMAL FREEZING EQUATION 


r.: 


We can now balance Che amount of component B in the remain- 

1 s 

Ing liquid by its chemical molar concentrations m^ and m^ 

then obtain the usual differential equation for normal freezing 
(Refs. 9, 11-13) as follows 


dm* , 

2 — . — 2 £_ 

1 s 1 - p U/J 

m B * “B 

where p is the fraction solidified. 

Substituting into the normal freezing differential equation 
results in 


; > i 


' f i 

i I ’ 


where 


- 1 ”< 1 - p) - : gfe ia/T) 


L_1 

LU - 


U) T (1 - V) T (1 


+ n-rr- d(i) 


q - 1 


1 - S 


1 - r 


D-ll 


and 


q - exp(f - ga/b) ' 

ag - bf 
r " exp -J“TtT 

s - exp kC- ~. a & > 

U Q - exp(f + g/T o ) 

V o - exp(a + b/T o ) 


( 19 ) 


r 1 - u ] 
0 

a 

fl - V 

o 

e 

[1 - v /u 
0 0 

1 - U 


t l - V 


1 - v/u 


( 20 ) 


where T q is the initial freezing temperature and can be deter- 
mined from the pseudoliquidus line on substituting for x^ g the 
initial melt composition, x* ■ 2 • x^ 
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THE Ga-Al-As SYSTEM 


This ternary system can be considered as a GaAs-AlAs pseudo* 
binary. In the liquid melt, we have 


x. 


As 


m 


As 


0.5 

1 


1 , 1 1 , 

X B + X C " X Ga + X A1 


1 . 1 
1 B + m C 


1 . 1 
m Ga + m Al 


0.5 


In the GaAs-AlAs mixed crystal, we have 


GaAs 


AlAs 


2m 


Ga 


2m 


Al 


s s 

X GaAs X AlAs 


2(m Ga + “Il> 


1.0 


To determine the phase diagram [i.e., liquidus and solidus, 
Eqs. (8) - (15) J and, in addition, to compute the normal freezing 
behavior according to Eqs. (18)-(20), the following four thermo* 
dynamic constants are required (Ref. 14 gives the values of these 
constants) : 


Melting point of GaAs 
Melting point of AlAs 
Heat of formation of GaAs 
Heat of formation of AlAs 


1511 K 
2013 K 

16.64 Kcal/mol 
22.8 Kcal/mol 
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Hence, a ■ 5.700, b - -11,470, f • 5.542, g ■ -8,374, q ■ 3.983, 
r » 0.1505, s ■ 166.5, a » -1.335, 0 * 0.1772, and y * 1.006. 
Table 1 summarizes the solidification behavior of the Ga-Al-As 


system, treated as a ^seudob inary . The first three columns describe 

the phase diagram, i.e., the concentration of aluminum in the liquid 

1 s 

melt and that of AlAs in the solid crystal x ^^ As at tem- 

peratures from 1720°C down to 1238°C at 20°C intervals. The 

next two columns, obtained from the first three, show that the 

s 1 

segregation coefficient k for Al, i.e., = x A ]_ As /2x 

g 

always exceeds unity while that of Ga, i.e., k^ a = (1 - x A ^ As )/ 

(1 - 2x^) , never does so except at 1238°C. Both of these segre- 
gation coefficients, however, are widely variable, increasing with 
decreasing melt temperature, by factors of over 6.3 in the case 
of k.^ and of over 3.7 in the case of kg fl . Thus, for this 
system, as well as for many others, the linear normal freezing equa- 
tions (Refs. 1-9) involving constant segregation coefficients are 
generally not applicable. 

These phase diagram data can, of course, be obtained also by 
means of the equations of Ilegems and Pearson (Ref. 14), or those 
of Steininger (Ref. 15). The direct computation of the normal 
freezing data from basic thermodynamic constants, however, appears 
to have never been attempted or published. By means of Eqs. (18) 
and (19) of this memorandum, these normal freezing data are cal- 


culated and given as the last nine columns. Here, the fraction 
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solidified, p, is given as a function of both the melt tempera- 

ture, T, and the A1 concentration in the original liquid melt. 

The concentration profile can be obtained by plotting the 

concentration of AlAs in the solid crystal, x ^]^ s > against the 

fraction solidified, p, for the appropriate initial A1 concen- 

1 


tration in the melt, x 


Al,o* 
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OTHER SYSTEMS 


Many other ternary systems can also be treated as pseudo- 
binaries. The necessary thermodynamic constants are usually 
available in the literature. Steininger (Ref. 15), for example, 
has examined a number of important systems and given the neces- 
sary thermodynamic constants required to compute the solidifica- 
tion behaviors and the phase diagrams. These systems include: 
InAs-GaAs, InSb-InAs, InSb-GaSb, GaSb-AlSb, InSb-AlSb, CdTe-ZnTe, 
CdTe-CdSe, HgTe-CdTe, GeTe-MnTe, PbBr 2 -PbCl 2> InAs-InP, SnTe-PbTe, 
PbTe-PbSe, HgTe-HgSe, ZnTe-ZnSe, and HgTe-ZnTe. For the first 
13 systems, his calculated solidus curves are in good to excellent 
agreement with the published experimental values. These systems 
can, therefore, be considered ideal. One can thus compute with 
some confidence their normal freezing behaviors according to 
Eqs. (18) -(20) of this memorandum. 
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ABSTRACT 


In the study of normal evaporation, it is assumed that the 
evaporating alloy is homogeneous, that the vapor is instantly 
removed, and that the alloy follows Raoult's law. The differen- 
tial equation of normal evaporation relating the evaporating time 
to the final solute concentration is given and solved for several 
important special cases. Uses of the derived equations are ex- 
emplified with a Ni“Al alloy and some binary iron alloys. The 
accuracy of the predicted results are checked by analyses of 
actual experimental data on Fe-Ni and Ni-Cr alloys evaporated 
at 1600°C, and also on the vacuum purification of beryllium. 
These analyses suggest that the normal evaporation equations pre- 
sented here give satisfactory results that are accurate to within 
an order of magnitude of the correct values, even for some highly 
concentrated solutions. Limited diffusion and the resultant sur- 
face solute depletion or enrichment appear important in the ex- 
tension of this normal evaporation approach. 
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INTRODUCTION 


In normal evaporation, we assume that the evaporating alloy 
is always homogeneous in composition, that the vapor is instantly 
removed, and that the alloy follows the Raoult's law (Ref. 1). 

Such conditions exist or are approached in an induction-stirred, 
melt-in-vacuum or liquid drop- in- space. 

This memorandum deals with the normal evaporation of binary 
alloys. In particular, we study the evaporative segregation pat- 
terns, i.e., thi. :ype and degree of enrichment or depletion of 
solute in the evaporating source, at different evaporation tempera- 
tures and times. 


DIFFERENTIAL EQUATION OF NORMAL EVAPORATION 


The appendix of Ref. 2 gives the exact solution for the normal 
evaporation of binary alloys. This equation relates the concentra- 
tion of the solute, m, at a time, t, when the mole fraction of 
the alloy remaining is F, as follows: 


F 




_JL 

U-V.l - 


1 - 


m 


_LL 

TJ “V 


Sl\ 


) 


( 1 ) 


where N and N q are, respectively, I s . .lumber of moles of both 
solvent and solute at evaporating times t * t and t * 0; m is 
the initial molar concentration of the alloy; and U and V are, 
respectively, the evaporation rates of pure solute and solvent. 

For pure solute and solvent, respectively, these evaporation 
2 

rates in mol/cm /sec are 


/4 D / 1 I 

U - K 10 U U (M u T)~ (2) 

and 

A -B /T i 

V - K 10 V V (M y T) * (3) 

• 5 

where K ■ 5.833 x 10 a, a ■ 1 for most metals (Ref. 1); M u 
and M v are, respectively, the molecular weights of the solute 
and solvent atoms; T is the evaporating temperature in degrees 
Kelvin; and A^ and B u or A^ and B y are evaporating con- 
stants for the solute or solvent given, for example, by Ref. 1. 
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Differentiating Eq. (1) yields 


but 


dm 


QL“. vM i _ , m ) dN 

[mU + (1 - m)V]N 


(A) 


dN - - [mU + (1 - m) V ]A • dt (5) 

where A is the evaporating area of the alloy, assumed coi.^tant 
here. 

Substituting Eq. (5) into Eq. (4) results in 

dt ■ Gm a (l •• m)^dm (6) 

where 


G 


„ . a+2 

N (1 * m ) 

-fi a 

A(U - V)m“ +1 


(") 


a - (2V - U)/(U - V) (8) 

3 - (V - 2U)/(U - V) (9) 
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SPECIAL CASE SOLUTIONS OF THE 


NORMAL EVAPORATION EQUATION 

Equations (6) -(9) allow us to determine the evaporating time, 
t f for an alloy to reach a specific solute concentration, m. Un- 
fortunately, these equations are not exactly solvable in the gen- 
eral case. All of the following special and important cases, (ex- 
cept Case V), however, are solvable in closed forms. 

Case I: The solute is much more evaporative than the solvent, 

i.e., U » V; or a » -1 and 0 ■ -2. In this case 

G.dm 

dt l 1 2 < 10 > 

1 m(l - m 

where 


G 


1 


V 1 ~ -o> 

AU 


and 


in 


(1 ' 


m )m 
JL 


m - m 


m (1 
o' 


- >») (1 - 


m o> (1 


- m) 


(ID 


When m ~ m Q , the second term in the bracket is nearly 0. 

The evaporation time for this first case, t^, is then a logarith- 
mic function of m, as has been experimentally observed (see e.g.. 
Ref. 3). 


Also, when m ~ m Q ~ 0, 


N m 


( 12 ) 
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Case II: The solvent is much more evaporative than the 
solute, i.e., V » U; or a • -2 and 0 - -I. In this case 


dt 


G 2 dm 


^ (1 - m)m^ 


( 13 ) 


where 


and 


m N 


AV 


in 


(1 ■ m )m m - m 


m (1 - m) m ie 
o v ' o 


(14) 


For dilute solutions (i.e., m ~ m Q ~ 0) , also becomes a 

logarithmic function of m, as has been observed. 

Case III: a » (2V - U)/(U - V) - 0, i.e., U « 2V and 0 - -3 


t 3 - C 3 [(l - m)' 2 - (1 - m o )‘ 2 j 


(15) 


where 


G 


3 



2 AVm 
o 


Case IV: 0 - (V - 2U)/(U - V) - 0, i.e., V - 2U and a - -3. 

In this case 




( 16 ) 


where 


N m 2 

r _ o o 

G 4 " 2AU(1 - m o ) 

Equations (15) and (16) show that under some conditions, the 

evaporating time, t, is more adequately represented by linear 

-2 -2 

functions of (1 * m) or m , rather than by logarithmic func- 
tions of m. 

Case V: For relatively dilute alloys, i.e., m and £3m C'' 1, 

the following solution by series expansion can be obtained from 
Eq. (6): 


dt r * Gm 


a 


(Pm) ^ ^ (Pm) *- 

1 - pm + 2 j + * * * + (-1) 


(17) 


and 


t 5 -G 


a+1 a+1 

m - m 
Q 


a+1 


a+2 a+2 a+i+1 a+i+1 

m - m m - m 

2 ~ P+ * * * + ( _1 )( a + i + ?)i! 


a + 2 


(18) 


For computer calculations, it is desirable to know the ratio 
f~Vi f-u 

of the l term to the (i - 1) term, thus 


T. P(a + l)(m a+i+1 - m a+i+1 ) 


T. , * . , . n / a+i a+i\ 

i*l i(a + i + l)ym - m J 


(19) 


which is generally less than Pm/i or em 0 /l. 


E-10 


r 



Because 0m « 1 and 1 constantly increases with each addi- 
tional term, this series converges rapidly unless 0 is very large, 
i.e., unless U * V, which leads to the following interesting case: 

Case VI; The solute and solvent are evaporating a t equal rates, 
i.e., U * V. In this case, we would expect 

m « m for all t . (20) 

o 

There is, then, no evaporative segregation, that is, there is 
neither solute enrichment nor depletion in the evaporating source. 

For any pair of solvent and solute, there is a unique tempera" 

ture, T , at which U * V and, hence, the alloy concentration 
s 

remains stable or constant. Equations (2) and (3) give: 


B - B 
_U iL 


s \ ' A v ‘ °- 5 lo *<W 


( 21 ) 


Case VII: With extremely dilute alloys, i.e., m ^ m Q % 0, 

we have (Ref. 2) 


F - N/N q = (m/m Q ) 


a+1 


( 22 ) 


Hence, 


Hat 


a+1 


*7 " AVI 1 ' ( " /m o> J 


(23) 
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EXAMPLES OF COMP UT AT ION 


• The Ni-Al System: 

As an example of the use of the various derived equations, 

the evaporation behavior o£ an alloy containing 8 percent by weight 

of A1 in Ni at the melting point of pure Ni (i.e., 1453 C) 

is computed. Here, the solute element (Al) is compare clvely 

highly evaporative relative to the solvent (Ni) . Equation (11), 

therefore, applies, and the time, t^, to reach a final solute 

concentration m from a specified initial concentration, m^. is 

directly proportional to N q A/U 'in the constant) . Table 1 

gives the times to reach various final Al concentrations for one 

mole (53.66g) of the 8 percent Al in Ni alloy (m = 0.159) 

° 2 

evaporating at 1453 °C from its (supposedly constant) 10 cm 
surface . 


Table 1 

NORMAL EVAPORATION OF 8% BY WEIGHT ALUMINUM IN NICKEL AT 1453 C 


Final Solute Evaporation 

Concentration Alloy Remaining Time (sec) 


Weight 

Fraction 

Mole 

Fraction 

Weight (g) 

Mole 


0.080 

0.1590 

53.66 

1.0000 

0.00 

0.079 

0.1573 

53.60 

0.9977 

14.08 

0.078 

0.1555 

53.53 

0.9955 

28.29 

0.077 

0.1536 

53.47 

0.9932 

42.64 

0.076 

0.1518 

53.41 

0.9910 

57.12 

0.075 

0.1500 

53.34 

0.9887 

71.75 

0.074 

0.1481 

53.28 

0.9865 

86.54 

0.073 

0.1463 

53.22 

0.9842 

101.46 

0.072 

0 . 1444 

53.15 

0.9820 

116.55 

0.071 

0.1426 

53.09 

0.9798 

131.80 

0.070 

0.1407 

53.03 

0.9775 

147.22 
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• Iron Alloys: 

The evaporation behavior of binary iron alloys containing 20 
different solute elements has also been studied. Table 2, listing 

Table 2 

EQUI-EVAPORATIVE TEMPERATURES FOR 
TWENTY BINARY IRON ALLOYS 



U/V at 




U/V at 

T ' C 

Solute 

1600 

°C 


V c 

Solute 

1600 c C 


S 

Cd 

3.765 

X 

io 6 

13,540 

W 

4.412 x 

io‘ 12 

-69,590 

Zn 

2.065 

X 

10 6 

15,610 

C 

8.210 x 

10" 8 

5,263 

Mg 

9.074 

X 

10 5 

19,610 

Mo 

1.879 x 

io" 7 

-12,130 




4 




-6 


Ca 

5.503 

X 

10* 

9,335 

Zr 

1.950 x 

10 

“47,790 

Pb 

3.333 

X 

10 3 

4,976 

B 

6.832 x 

io" 5 

9,510 

Mn 

1.070 

X 

10 3 

20,740 

V 

3.804 x 

io -3 

8,574 

A1 

45.65 



7,918 

Ti 

2.254 x 

10' 2 

34,650 

Cu 

12.25 



5,613 

Co 

4.362 x 

io" 1 

4,318 

Sn 

9.978 



2,686 

Si 

5.238 x 

io' 1 

2,825 

Cr 

3.158 



“215 

Ni 

5.897 x 

io" 1 

3,036 


the equi “evaporative temperatures for these 20 different alloy sys- 
tems, also gives the ratios of the solute evaporating rate, U, at 
1600°C, to that of the solvent iron, V. At 1600°C, 10 of these 

solute elements evaporate faster than the solvent (the three left 

columns) and 10 slower (the three right columns) . Moreover, these 

-12 

ratios vary widely over 18 decades, from 4.41 x 10 for the 
slowest evaporating, W, to 3.76 x 10 8 for the fastest evap- 
orating, Cd. Because of this wide variation in evaporating rates, 
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and because of the extreme sensitivities of the evaporating surface 
to unsuspected contaminants, predicted or experimental evaporating 
results cannot generally be very accurate. 

The equi -evaporative temperatures in iron alloys also vary 
widely. Binary iron alloys containing Cr, Zr, Mb, and W have 
no practical equi -evaporative temperatures. One can, therefore, 
always expect these alloys to change compositions continuously 
with the evaporating time. 

Table 3 shows the effect of evaporating temperature on the 
U/V ratios for four different solute elements Mg, Ca, Mn, and 
Al. In the range of 1500 °C to 1900 °C and beyond, increasing 

Table 3 


EFFECT 

OF SOLUTE 

ELEMENTS AND 

EVAPORATING TEMPERATURES 

ON 

THE VALUES 

OF U/V (x 

1000) 

IN IRON 

ALLOYS 



Solutes 


Temp 

(°C) 

Mg 

Ca 


Mn 

Al 

1500 

2132.0 

118.2 


1.649 

0.06037 

1600 

907.4 

55.03 


1.070 

0.04565 

1700 

421.2 

27.68 


0.7261 

0.03552 

1800 

210.5 

14.88 


0.5110 

0.02831 

1900 

112.1 

8.466 


0.3719 

0.02304 


the evaporating temperatures always decreases the U/V ratios. 
This can also be seen from Table 2, as the equi-evaporative temper 
atures shown for these four binary iron alloys are higher than 
1600 °C, at which temperature the four solute elements evaporate 
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much faster than the solvent Iron. Table 3 also shows that for the 
same temperature variation, the more evaporative the solute element, 
the more percentage variation the U/V ratio. In the case of Mg, 
e.g., the U/V ratio decreases by about 20 times from 1500°C tc 
1900 °C, whereas for Al, the same ratio decreases only by less 
than three times over the same temperature interval. 

Figure 1 shows the effect of solute elements and evaporating 

temperature, T, on the evaporating time t, for a given set of 

initial and final solute concentrations (i.e., m * 0.01 and 

o 

m ■ 1 ppm) . For given m Q and m, the log t versus 1/T curves 
for these highly evaporative solute elements are approximately lin- 
ear and have positive slopes. This can be expected from Eq. (11) 
since, for given m Q and m, log t is linearly related to 

B/T “ 0.5 log T, and since 0.5 log T is small relative to B/T 

within the evaporating temperature range studied. Thus, one can 
determine, from Fig. 1, the value of the elemental evaporating con- 
stant B, or heat of evaporation AH * 4.574 B (Ref. 1), for the 
solute elements by plotting log t versus 1/T and measuring the 
slope of the resultant, nearly straight lines. 

In Fig. 2, the log t versus 1/T relationships also appear 
nearly linear for all the four different initial concentrations 
(i.e., m Q - 10 10 10 \ and 10 ^) of Al in Fe. This 

can also be seen from Eqs. (11) and (12). All these nearly straight 
lines have identical slopes, from which the heat of evaporation of 
pure Al can be evaluated. 

Figures 3 and 4 display the effect on the evaporating time, 
t, of final solute concentration, m, and either initial concen- 
tration m Q (in Fig. 3) or solute elements (in Fig. 4). 

Figure 5 indicates that Mg, Ca, and Mn in Fe alloys are 
so evaporative at 1600 °C that practically all of these elements 
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1900 1800 1700 1600 1500 


EVAPT, °C 


SOLVENT: Fe 

m =0.01 

o 

m = 1 ppm 


Fig. 1 Evaporating Times for Different Dilute Iron Alloys 

to Change Concentration from itu, = 0.01 to m = 1 pj 





1900 1700 1500 


EVAPT, °C 


Fig. 2 


Evaporating Times for Dilute Iron Alloys Containing 
A1 and Starting at Four Different Initial 
Concentrations, to Reach a Final Concentration of 
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• SOLVENT: Fe 

• SOLUTE: Al 

• TEMP * 1600 °C 


EVAP t, 
SEC 



10 ” 10 v 10 
FINAL CONC, m 


Pig. 3 Effect of Initial and Final Solute Concentrations on 

the Evaporating Times in Dilute Iron Alloys Containing 
Al at 1600 °C 
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Fig. A Effect of Solute Elements and Final Concentration 

°" t ??jY aporatin S Ttmes in Dilute Iron Alloys 
at 1600 C 








FINAL CONC, m 


Fig. 5 Evaporation Loss of Dilute Iron Alloys at 1600 c C 
As a Function of Solute Elements and Final Solute 
Concentrations 
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are removed by evaporation, without much evaporative loss of the 

solvent Fe atoms. On the other hand, A1 is comparatively less 

evaporative so that much of the solvent Fe atoms are evaporated 

off together with Al. To achieve a purification factor of 10^ 

— 9 

(i.e., to m ■ 10 ) from m Q ■ 0.01, for example, the initial 

evaporating alloy must lose over 30 percent of its material. 
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ACCURACY OF PREDICTED RESULTS 


To check the validity of our derived equations for normal 
evaporation of binary alloys, the current literature has been 
searched. Several sets of alloy evaporation data have been found 
that are amenable to normal evaporation analysis. These sets in- 
clude the following: 

• The Fe-Ni and Ni“Cr Systems: 

An analysis has been made of the data by Obradovi , et al. 

(Ref. 4) for Fe-Ni and Ni~Cr alloys evaporated at 1600 C 

for different times under various ambient pressures. In the 

80% Ni-20% Cr case, evaporation started with N = 0.7582 moles 

° 2 

of the alloy having an evaporating area A = 10.0 cm . The solute 

and solvent evaporating rates are, respectively, U = 3.479 x 10 

“6 2 

and V... = 6.386 x 10 mol/cm /sec. This is a case where the 
Nl 

solute, Cr, is comparatively highly evaporative relative to the 
solvent, Ni. Hence, the evaporating time t^ given in Eq. (11) 
applies, i.e., 


f (1 - m )m m - m I 

h G lp m o (l - m) + (1 - m o )(l - m) j 

where G. * ~N (1 - m )/AU, and m and m are the initial and 
1 o o o 

final molar solute concentrations in the alloy. 

Least square fits of the Obradovic data give nearly constant 

A 

values of observed G, : -1.785, -5.064, and -5.064 (x 10* sec) 

1 -3 

for ambient pressures of 1, 100, and 500 (x 10 torr) . The cal- 
culated values of G^, though also nearly constant, are, however, 
one order of magnitude smaller, indicating surface solute depletion. 
Similar analysis for the 60% Fe-40% Ni alloy evaporated with 
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N * 0.7478 moles and A • 10.0 cm^ gives U... * 6.386 x 10 ^ 
o - 5 2 Nl 

and Vp e * 1.081 x 10 mol/cm /sec. Here, the solvent Fe is 

much more evaporative than the solute Ni, and Eq. (14) aprlies. 

The values of observed evaporating coefficients G 0 * m N /AV are 

A ^ O O 

2.469, 1.369, and 8.386 (x 10 sec), for ambient pressures of 

-3 

1, 100, and 500 (x 10 torr) , respectively. The calculated 
G 2 *s are again nearly constant but also an order of magnitude 
smaller than the observed G 2 's, again indicating solvent deple- 
tion (or solute enrichment) at the surface. Details of these 
analyses are given in Ref. 5. 

• Beryllium Purification 

The kinetics of normal evaporation for beryllium have been 
quantitatively checked with actual results of beryllium purifica- 
tion during vacuum induction melting. Details are given in Ref. 6. 

In these tests, beryllium was crucible- free, induct ion -melted under 

“6 

an ambient pressure of 10 torr. The actual temperature of the 
melt was not known, and the exact evaporating rates for the beryl- 
lium solvent and various solute elements cannot be computed. How- 
ever, because all solute elements (i.e., Fe, Cr, Mn, Ni, Si, Al, 
Mg, Cu, Zn, and Na) and the solvent. Be, were evaporating 
from the same or common liquid-gas interface of a fixed area for 
the same length of time, we can compute the values of P, defined 
as the product of the evaporating time, t, and solvent evap- 
orating rate, V (for V » U) , or solute evaporating rate, U 
(for- U » V) ; thus: - P * tV cr tU. -Table A., gives the -initial 
concentrations and final concentrations (In ppm) of the various 
solutes, together with the value of P, actual surface concentra- 
tions, ratio of actual surface to bulk concentrations, and effec- 
tive times to reach the final concentrations under the assumption 
that the evaporating temperature was 1250 c C. The following con- 
clusions can be drawn from Table 4: 
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1. The solute elements can be divided into three groups: 

Fe, Cr, Ni. fi, and Cu evaporate much more slowly 
than Be; R.-,, Zn, and Na evaporate much more rapidly 
than Be; A1 evaporates at about the same rate as Be. 

2. The computed times to reach the final concentrations 

under normal evaporation conditions (complete liquid 

mixing) are fairly constant for the solute elements 

that evaporate much slower than Be (i.e., Fe, Cr, 

4 

Ni, Si, and Cu) , being about 4 x 10 . 

3. The computed times to reach the final concentrations 
for the highly evaporative Mg, Zn, and Na are also 
fairly constant, but about four orders of magnitude 
smaller. 

4. After correction for limited liquid mixing (Ref. 7), 
the effective times to reach the final concentrations 
are much more constant, even between the groups of 
solute elements. In particular, the highly evapora- 
tive elements Mg> Zn, and Na have their results 
improved by several orders of magnitude. 

5. Limited liquid diffusion must, therefore, be con- 
sidered, particularly for the highly evaporative 
elements. 

6. The surface solute concentrations can be changed by 
up to six orders of magnitude, so that the effective 
evaporating rates can be similarly changed. 
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DISCUSSION AND CONCLUSION 


These and other similar analyses suggest that the normal evap- 
oration equations presented in this Memorandum work surprisingly 
well for the several alloy systems studied so far. These equations 
generally give results accurately to within an order of magnitude 
of the correct values, perhaps even for highly concentrated solu- 
tions (e.g., 20% Cr in Ni or 40% Ni in Fe) . Such accuracies 

are usually sufficient for many evaporation studies. 

In other more critical studies, however, refined evaporation 
analyses may be needed. In particular, the effect of limited dif- 
fusion and, hence, surface depletion or enrichment of the solute, 
must often be accounted for. This is true especially for the 
highly evaporative solutes such as, e.g., Hg, Zn, and Na in 
the solvent Be. We already have results that confirm this con- 
clusion. Details of these results will be published in Ref. 7 and 
elsewhere . 
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A NEW MECHANISM FOR THE FORMATION OF GROWTH SPIRAL 
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ABSTRACT 


A new mechanism is proposed to explain the formation of 
growth spirals in controlled melt-grown, crystalline samples of 
Ba^SrNb^O^Q , Several distinguishing features of these spirals 
cannot be explained on the basis of the classical, core-first 
growth theory. These features are examined in detail, especially 
in comparison with the inherent or implied features of the classi- 
cal spiral model. Characteristics of the spirals are predicted 
based on solidification theory, and are found to be in agreement 
with the observed results. 
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J. Applied Physics 

* Research Department, Grumman Aerospace Corporation, Bethpage, 
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Introduction 


The spiral growth mechanism was first proposed by Frank 
(Ref. 1) in 1949 and has been developed into the classical BCF 
model (Refs. 1*3). Spiral markings had, however, long been ob- 
served on silicon carbide crystals (Refs. 4-5). Since 1949, many 
more spirals have been observed during evaporation-condensation 
(Refs. 6-7) and solidification (Refs. 8-10). Some of these spirals 
are rounded (Refs. 9) while others are segmented (Refs. 8-9). 

The BCF model has been able to explain many of the observed 
effects and growth features. However, in our controlled crystal 
growth experiments on Barium Strontium Niobates (BSN) , we observed 
some new phenomena that cannot be easily explained by the classi- 
cal BCF model. In fact, some of these observed phenomena contra- 
dict the same model. This paper briefly describes the implications 
of the classical BCF model and some experimental results recently 
observed on these niobates in our laboratory. From these results 
we propose a new spiral growth model which seems to explain better 
the observed effects on the niobate materials. 

Existing Theory 

According to the BCF model of spiral formation, a nucleus is 
formed at a favorable spot on the solid substrate surface, and ad- 
ditional atoms then solidify or condense from the liquid or vapor 
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phase to spiral around this nucleus, i.e., to form the spiral. 

The distinguishing features of this theory are: 

1. The spiral starts at its center or core region 
in the form of a screw dislocation and spirals 
out. 

2. Adjacent or successive turns of the spiral should 
be separated by a thin boundary no thicker than a 
subgrain boundary, i.e., angstroms thick, since 
even the grain boundary in a polycrystalline mate- 
rial is generally believed to be only 10A or 
less thick (Ref. 11). 

3. The primary temperature gradient is perpendicular 
to the plane of the spiral (Ref. 10). 

4. The step height in spirals, at least according to 
the BCF model, should be only one atomic layer 
thick. 

Experimental Procedures 

By using standard sintering techniques, polycrystalline 
samples of hypoeutectic composition (ABal^Og : Srltt^O^) were 
prepared (Ref. 12). The raw materials used in preparing these 
ceramic samples were: barium carbonate, from J. T. Baker, re- 

agent grade; strontium carbonate, from J. T. Baker, reagent grade 
and niobium pentoxide, from A. D. Mackay, 99.9 percent. The 
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sintering was done at 1350°C for more than 80 hours in air and 
the sintered material was checked for the completion of chemical 
reaction by using an IRD-Guinier camera. 

Directional solidification of the sintered compacts was car- 
ried out in a 4 in. long, £ in. diameter alumina crucible heated 
in a vertical tube-type platinum furnace. The temperature gradient 
was varied from 20-40°C/cm by changing the tapping of the re- 
sistance wire of the furnace, and the temperature was monitored by 
using three thermocouples placed at three different places. After 
heating the material to about 30-40°C above its melting tempera- 
ture, the alumina crucible was withdrawn from the furnace at a 
constant rate of from 1-40 ^m/sec. The longitudinal and trans- 
verse sections on the resultant, elongated samples were then ex- 
amined on the solidified specimens. 

New Experimental Results 

The following observations have been made on these melt- 
grown, crystalline niobate samples prepared at 40 ixm/sec: 

1. The spirals consist of linear oriented segments, 
typically 150 n,m or shorter. 

2. Many spirals have only 1, l£, 2, 3, ..., outer 
linear segments (see Fig. 1 at a) . 

3. The linear segments have been identified by elec- 
tron microprobing as proeutectic crystals. 
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4. The spacings (Figs. 1-3 at b) between adjacent 
linear segments of the spirals are as wide as, or 
wider than, the thickness of the proeutectic cry- 
stals, i.e., tens or hundreds of \im rather than 
angstroms thick. 

5. By probing, the materials in these spacings are 
identified as being of eutectic composition 
(Figs. 1-3 at b). 

6. Both single spiral (Fig. 2, top) and paired spirals 
(Fig. 3, right) occur. 

7. Clockwise and counterclockwise paired spirals exist, 
even only 70 pm apart on the same sample cross 
section (Fig. 3 at c) . 

8. The two adjacent, clockwise and counterclockwise 
spirals (Fig. 3) share the same single crystalline, 
outermost linear segment (Fig. 3 at d) . 

9. During spiral formation, there was a longitudinal 
temperature gradient of about 40°C/cm perpendicu- 
lar to the axis of growth spirals. 

10. The growth directions of the linear segments were 

often independent of the externally applied, longi- 
tudinal temperature gradient (see Figs . 2-3 at a) . 
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Discussion 


Many spirals are near ideal helices, each containing a single 
smooth spiraling segment with no discontinuities in the change of 
its slope (Ref. 9). Many other spirals are segmented (Refs. 8-9), 
similar to those noted in Observation 1. These segments apparently 
follow the preferred crystallographic orientations which probably 
coincide with the directions of fastest crystal growth or high 
thermal conductivities. 

The spirals with only 1, l£, 2, 3, ..., outer linear seg- 
ments (Fig. 1) appear to represent different stages of growth, 
from the outermost points toward the center or core regions. The 
hypothesis of complete spiral formation first with subsequent core 
dissolution or homogenization is untenable, because the core would 
be solid, according tc the classical theory, and, hence, have 
limited material diffusivities . Solid diffusion coefficients at 

near the melting point or crystal growth temperatures is typically 
-8 2 

only about 10 cm /sec, and the time to grow a 100 p.m seg- 
ment at a rate of 40 iim/sec is only 2.5 seconds. The diffusion 

length in this particular case is of the order of - Vl0~ 8 x 2.5 

* 

L6p.m, that is two orders of magnitude smaller than the segmental 
length . 

The fact that linear segments are proeutectic crystals (Obs. 3) 
is in accordance with the relevant phase diagram (Ref. 12). 
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Observation 4, i.e., the vide spacings between adjacent seg- 
ments, is in sharp contrast with prediction or implication of the 
conventional theory. 

Observation 5, i.e., the materials in the spacing between 
spiral segments are eutectics, indicates that the linear segments 
froze first, protruding into the melt (somewhat like dendrites) 
and rejecting high-solute materials ahead of the liquid-solid 
interface and also on one or both sides of the protruding seg- 
ments, to be subsequently frozen as eutectics. 

The interesting fact that adjacent clockwise and counter- 
clockwise spirals share the same single crystalline, outermost 
linear segment (Obs . 8) also appears to defy satisfactory explana- 
tions according to the classical, core-first BCF growth model. 

The statistical probability is almost nil for two independent, 
separated nuclei in the melt to form the two segmented, clockwise 
and counterclockwise spirals that grew to reach perfect alignment 
and matching widths in the last and common, outermost linear seg- 
ment. We would have to hypothesize that one nucleus must be trans- 
mitting complicated signals on growth characteristics, such as 
growth rates, growth directions, ... through the intervening 1' 'i d. 
This liquid is probably in a disordered state, so that at best, the 
signals are transmitted through with great distortion and absorp- 
tion. The other nucleus must simultaneously receive and decode the 
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signal and, furthermore, take instantaneous, exact corrective 
action to achieve perfectly matched growth conditions. Even if 
we assume that somehow the two nuclei giew from opposite direc- 
tions in an identically and inf .ranizaneously matching manner, we 
still expect the two approaching crystal segments to be separated 
by a solute-enriched layer, due to the limited liquid diffusion 
and solute pile-up in front of the growth interface. A loss or 
epitaxial relations must occur, i.e., two crystalline segments 
must form, one on each side of a distinct, solute-rich laye 
Such was not observed. 

It is also interesting to note that the external or primary 
temperature gradient is transverse to the axis of the spiral (Obs . 9) , 
rather than longitudinal thereto, in disagreement with the normal, 
spiral growth predictions. Some segments must, therefore, grow in 
the direction of this primary temperature gradient (at fast growth 
rates), some must grow against the temperature gradient (at 40°C/cm, 
or over 0.6°C over 150 pm), and some others must grow essentially 
isothermally. 

A New Mechanism 

It is, therefore, proposed that at least on the samples re- 
ported here, the spirals were formed from outside toward the center 
or core region. Such a mechanism has been suggested (Refs. 13-14). 
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All the above observations and discussions , including the 
formation of eutectics between the spiral spacings, can now be 
readily explained by the solidification theory (Ref. 15). Even 
the growth against a temperature of over 0.6°C can be ascribed 
to solute pile-up in front of the liquid-solid interface, and the 
resultant constitutional supercooling phenomenon (Ref. 16). 

The two adjacent clockwise and counterclockwise spirals, such 
as those in Fig. 3, most likely shared the same nucleus, which was 
probably located in the middle of the outermost segment and grew 
in opposite directions to form the common outermost linear segment 
for the two spirals. 

According to our core-last spiral growth mechanism, the com- 
bined effect of temperature and solute concentration fields makes 
it less likely for nucleation to occur at the center or core re- 
gions. Only under a special combination of phase diagram, mate- 
rial thermophysical properties, and external temperature fields 
(such as when the material in the core region has the highest melt- 
ing point) can core nucleation become possible. 

When nucleation does occur at an outer corner, the growth ap- 
pears to follow the outline of some kind of local concentration 
cell, often only slightly influenced by the externally applied tem- 
perature gradient (Figs. 2-3). This can be seen from the fact that 
in our crystal growth experiment, the external temperature gradient 
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was 40°C/cm, or a maximum of 0.6°C over the 150 p.m outer seg- 
ments. According to the liquidus of the Srlft^Og - pseu- 

dobinary phase diagram (Ref. 12) , such a temperature difference 
corresponds to a change in liquid-solute concentration of only 
about 0.3 mole percent. On the other hand, the concentration ef- 
fect from solute segregation and pile-up ahead of the solid-liquid 
interface is, according to the separation between liquidus and 
solidus on the same phase diagram, over 30 times greater. 

The proeutectic crystals or linear segments, thus, appear to 
grow mainly along the isoconcencration lines generally following 
the cell boundaries, with rejection of solute (segregation coef- 
ficient being less than unity here) both longitudinally ahead of 
the grown segments and laterally at least toward the center of the 
spirals. We could thus predict that for these core-last grown 
spirals, the solute concentration of spiral segments should gen- 
erally increase along the segments toward the last-freezing core 
region. Notice, also, that the differently oriented segments should 
have different growth rates and concentration profiles. The centers 
of the spirals must, in particular, have much higher solute concen- 
trations approaching the eutectics. 

Hence, at least in the BSN samples, the spiral growth seems 
to start from the outside inward. 
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Fig. 1 
Fig. 2 
Fig. 3 


Captions 


Broken Spiral Segments on Ba^Sr NbloC>3o Melt-Grown 
Sample Formed at 40 iim/sec under Temperature 
Gradient of 40°C/cm. Magnification 220X 

Segmented Spirals on Ba 4 Sr Nbio030 Me It -Grown at 
40 p.m/sec under Temperature Gradient of 40°C/cm. 
Magnification 220X 

Segmented Paired Spiral on Ba 4 Sr Nb^^Q Melt-Grown 
at 40 p.m/sec under Temperature Gradient of 40°C/cm. 
Magnification 220X 
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STEADY-STATE SOLIDIFICATION IN THE Ni-Sn SYSTEM 
INTRODUCTION 


It is well known that solute redistributes as an alloy solidifies from its molten 
state. Pfann (1962) and several other authors have expressed the instantaneous solute 
concentration c as a function of the fraction solidified, g, i.e. , 

c = kC o (l-g) k 1 , (1) 

where C is the initial concentration of solute in the melt, k - c /c, is the segregation 
o si 

constant, and c and c are the solute concentrations in the solid and liquid, respectively. 
si 

Equation (1) is based on the following assumptions: 

1. Diffusion in the solid is negligible. 

2. Th° mixing in the liquid is complete. 

3. k is constant. 

Frenkel (1946) pointed out that assumption (2) is hardly true because the diffusion 
coefficient is too small to assume completely mixing. In addition to assumptions one 
and three, Chalmers (1961) assumed that convection in the liquid is negligible and thai 
solute is transported by diffusion instead of being completely mixed. Then the con- 
centration of solute in the liquid, c^, is, under steady state conditions, a function of 
distance from the solid-liquid interface, x: 



where D is the diffusion coefficient in the liquid and R is the constant growth rate. 
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There is always an equilibrium temperature which corresponds to a certain 

solute concentration in the liquid phase. Chalmers (1964) pointed out that supercooling 

can occur if the temperature gradient, G, is controlled in such a way that the actual 

temperature in the liquid phase is lower than the equilibrium temperature. Figure G-l 

shows that no supercooling can occur if G is larger than the critical temperature 

gradient of supercooling G . If G is smaller than G , however, a distance of super- 

c c 

cooling, x, exists. 

The constant value of k in Eq (1) and (2) is only a simplified assumption, since 
we know that the solidus and liquidus curves in most phase diagrams are not linear. 

In this appendix we will investigate the errors involved in calculating the critical tem- 
perature gradient G and distance of supercooling x when a constant k is assumed. For 
convenience in comparison, three models are assumed. The first is a linear model in 
which both the solidus and liquidus curves are linear functions of the melt temperature. 
This model corresponds to constant k. noth the solidus and liquidus curves should, of 
course, pass through the point representing the temperature of solid-liquid equilibrium 
of the pure solvent, and the slopes should be obtained by a least-square fitting of the 
experimental data over the entire solidification range. The second model is a pseudo- 
linear model. This model is the same as the first one, except that the slopes of solidus 
and liquidus curves are the slopes, or tangents, at the melting temperature of pure 
solvent. The third model is a quadratic model for both the solidus ~.nd liquidus curves, 
with the same terminal restrictions that they precisely pass through the solid-liquid 
equilibrium points for the pure solvent and eutectic. The shapes of the curves are also 
determined by least-square fitting of the experimental data. 

THEORETICAL BACKGROUND 

To express the liquidus and solidus curves, C^ and C g , in terms of the melt 

temperature, we choose T , the equilibrium temperature of pure solvent, as the origin 

h 

of the coordinate, so that the constant terms in both liquidus and solidus drop out 

automatically. Then the relation between the new coordinate or modified temperature 

T and old coordinate T' is T = T_ - T' . 

E 
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Fig. G-l The Effect of R/D on x 
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For the linear model, the liquidus curve and the solidus curves can be expressed 


as 


and 


For the quadratic model; 


and 


c i = A i T 


C - B T 

s 1 


C S B 1 

k = — = — — = constant 

1 Cj a 2 


c r A i2 T + A 22 T 


C s = B 12 T f B 22 T 


. C s B 12 - B 22 T 

3_ V A 12^ A 22 T 


(3) 

(4) 

(5) 


( 6 ) 

(7) 

( 8 ) 


which is a function of the temperature, T. If the temperature T is small, i.e. , if the 

concentrations of solute c and c„ are small, then the higher-degree terms of ttm- 

s 1 

perature in both C and C are negligible, we thus obtain a new constant for k, i.e. , 
j l s 


B 


12 


2 A 


12 


(9) 


which is the pseudo-linear model. 

G c can be found by the chain rule relationship; 
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. u 

r*x 


.. ^ *C l >c 

lim r>T ill 


c x-.o Sx Sx / 


/ 


Both ^ C l I ~)x and * C jt j ?>T can be found from Eq (2), (3), and (6). Thus: 


• For the linear case, 

• For the pseudo-linear case, G 

• For the quadratic case, G 


RT 

cl 55 ~D 2 ’ (1 ' k l l 
RT 

« = (1-kJ 

c2 D ' 2 


R’i 


o A 12 + A 22 T o 

c3 ' D ‘ ' S> A 12 * 2A 22 T c 


RT 

Tr "- k 2 » 


t- 


( 10 ) 

( 11 ) 

( 12 ) 


A 12 (A 22 B 22 ) ~ 2 A 22 B 12 
A 12 ( A 12 _ B 12 ) 


•J 


where T 0 is the equilibrium temperature of solid-liquid interface expressed in the new 

coordinate system. The real temperature is T* = T - T . 

£ o 

The distance of supercooling can be found by introducing the equation for the actual 
temperature in the liquid phase, after transformation into the new coordinate system: 


T = T - Gx (13) 

o 

Thus: 

• From Eq (3), C = A, (T -Gx) (14) 

l 1 o 

■ From Eq (6), C = A 12 (T o~ G x * (15) 

• From Eq (6), C = A i 2 ( V Gx) + A 22 (T o“ Gx)2 (16) 
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Then x can be found by equating Eq (14), (15), and (16) to Eq (2) for any model, and any 
given R/D and T q . 


COMPUTATION SCHEME 

To find the values of A^, B , A^, A 22> B^, and B 22> 
programs were used which minimized the total error y = £ [ 


least-square computer 

T T 2 "1 2 

1_A 12 C " A 22 C I 
l l J 


Since we want both end points to be precise, we solve the following simultaneous 
equations for A and A : 


and 


-2 

aA 12 


1-A JX_ A -X . 2 

12 C , 22 C c< C ? A 12 


A Tu 2 

1.-^— * A22 ^ 
i U Jtu 


(17) 


where C^ u and T^ are the liquid concentration and temperature at the eutectic point 
respectively. The same procedure can be used to compute A , B. , and B. , B . and 

X X ££ 

the required data sets obtained from the phase diagrams. 

To find k's and Ge's. Once we have obtained those coefficients, we can compute 
k's and Gc's from Eq (5) and (8), and (10), (11), and (12), respectively. 

To compute the distance of supercooling, x. From equations (7), (16), and (2), 
we solve for x, for the quadratic case, as follows: 


_R_ 

A 12<V g *> ♦ A 22 ( V GX) 2 - (B 12 T o+ B 22 T 2 , a* - k - e 

0 


) =0 


( 18 ) 


Since the above equation is nonlinear, Newtons method can be used to find x by 


minimizing y, where 



! 


a 


2 1_k "d X 

y = A 12 (T o -Gx) ♦ A 22 ,T o -Gx) 2 - (B 12 T o *B 22 T o ) , 1 +T ^ e , 


For the linear case: 


R. 

I * l-k " D *’] 

1 + “C e J 


( 19 ) 


( 20 ) 


i . 


i! 


and for the pseudo-linear case: 


R 


l-k D 

y - A 12 < V Gx, - B 12 T o ( 1 + -j^ e , 


( 21 ) 


where k is the segregation coefficient at a given solid-liquid equilibrium temperature, 
o 


! ; 



COMPUTATION RESULT 

A Ni-Sn System, in which both the solidus and liquidus cur are parabolas, was 
chosen for the study. Sn is the solute and Ni is the solvent in this system. Ten equili- 
brium data points, as tabulated below, were obtained from the phase equilibrium diagram 
given by Hansen (Ref. 4): 

TABLE G-l THE Ni-Sn PHASE DIAGRAM 


T, K 

C 

s 

C 

l 

273.1 

0.0 

0.0 

283.1 

0.2000E 01 

0.2000E 01 

298.1 

0.8500E 00 

0.4000E 01 

323.1 

0.1800E 01 

0.6900E 01 

373.1 

0.3700E 01 

0.1220E 02 

423.1 

0.5900E 01 

0.1700E 02 

473.1 

0.8600E 01 

0.2150E 02 

523.1 

0.1200E 02 

0.2600E 02 

573.1 

0.1700E 02 

0.3000E 02 

604.1 

0.2000E 02 

0.3230E 02 


I! 


% 




From these data, the required coefficients can be obtained by a least-square method: 

A12 A22 B12 B22 A1 731 

0.1495 -1. 568x1 0~ 2 0.02337 1.120xl0~ 2 0.1157 0.03364 

The segregation constants k's for the three different cases, i.e., k(l), k(2), and 
k(3) corresponding to the linear, pseudo-linear, anc quadratic models, respectively, 
were computed for a given T and are listed in Table G-2. 

TABLE G-2 VARIATION OF THE SEGREGATION COEFFICIENTS 


WITH TEMPERATURE FOR THE Ni-Sn SYSTEM 


H 

m 

_l 

K(l) 

K(2) 

K(3) 

1728 

0.2908E 00 

0.1563E 00 

0.1563E 00 

1718 

0.2908E 00 

0.1563E 00 

0.1655E 00 

1703 

0.2908E 00 

0.1563E 00 

0.1797E 00 

1678 

0.2908E 00 

0.1563E 00 

0.2045E 00 

1628 

0.2908E 00 

01563E 00 

0.2583E 00 

1578 

0.2908E 00 

0.1563E 0G 

0. 316bE 00 

1528 

0.2908E 00 

0.1563E 00 

0. 38 7 3E 00 

1478 

0.2908E 00 

0.1563E 00 

0.4656E 00 

1428 

0.2908E 00 

0.1563E 00 

0.5559E 00 

1397 

0.2908E 00 

0.1563E 00 

0.6192E 00 


The deviations of segregation constant k for the linear models from that of the 
quadratic model are given in Table G-3. 


TABLE G-3 DEVIATIONS OF THE SEGREGATION COEFFICIENTS 



[N THE Ni-Sn SYSTEM 

T, K 

(K(L)-K(3))/K(3) 

(K(2)-K(3))/K(3) 

273.1 

0.8604E 00 

0.0 

283.1 

0.75S7E 00 

-0. 5574E-01 

298.1 

0. 61 78 E 00 

-0.1304E 00 

323.1 

0.4221E 00 

-0.2356E 00 

373.1 

0.1258E 00 

-0.3949E 00 

423.1 

-0.8790E -01 

-0.5097E 00 

473.1 

-0.2493E 00 

-0.5965E 00 

523.1 

-0.3755E 00 

-0.6643E 00 

573.1 

-0.4770E 00 

-0.7189E 00 

SHU 

-0.5304 E 00 

-0. 7476E 00 


The critical temperature gradients, G , were computed by Eq (10), (11), 

c 

and (12) for any given R/D and T and are listed in Table G-4. Computed distances 
of supercooling, x, are given in Table 5. 
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TABLE 0-4 COMPUTED CRITICAL TEMPERATURE GRADIENTS 
IN THE Ni-Se SYSTEM FOR THE THREE MODELS 


T. K 

R/D 

G c (l) 

G c (2) 

C c (3) 

273.1 

0.1000E n 

0.0 

0.0 

0.0 

283.1 

0.1000E 01 

0.7082E 01 

0.8437E 01 

0.8434E 01 

288.1 

0.1000E 01 

0.1773E 02 

0.210SE 02 

0.2107E 02 

323.1 

O.tOOOE 01 

0.3546E 02 

0.4218E 02 

0.4211E 02 

373.1 

O.IOOOE 01 

0.7092E 02 

0.8437E 02 

0.8402E 02 

423.1 

O.IOOOE 01 

0.1064 E 03 

0.126CE 03 

0.1256L 03 

473.1 

0.1000E 01 

0.141NE 03 

0.1687E 03 

0.1668E 03 

$23.1 

O.IOOOE $1 

0.1773E 03 

0.2108E C3 

0.2073E 03 

$73.1 

O.IOOOE 01 

0.212ME 03 

0. 2531 E 03 

0.2464E 03 

604.1 

O.IOOOE 01 

0.2348E 03 

0.2793E 03 

0.2K83E 03 

273.1 

0.2500E 01 

0.0 

0.0 

0.0 

283.1 

0.2500E 01 

0.1773E 02 

0.2109E 02 

0.21 (WE 02 

288.1 

0.2300E 01 

0.4433E 02 

0.S273E 02 

0.5260E 02 

323.1 

0.2S00E 01 

0. 8H65E 02 

0.I055E 03 

0.1053E 03 

373.1 

0.2500E 01 

0.1773E 03 

0.2109E 03 

0.2I00E 03 

423.1 

0.2500E 01 

0.2660E IS 

0.3164E 03 

0.3141E 03 

473.1 

0.2500E 01 

0.3546E 03 

0.4218S 03 

0.4171E 03 

$23.1 

0.2500E 01 

0.4433E 03 

0.5273E 03 

0.51M2E 03 

573.1 

0.2580E 01 

0.5319E 03 

0.6328E 03 

0.6158E 03 

604.1 

0.2S0OE 01 

0.5H69E 03 

0.68H2E 03 

0.6733E 03 

273.1 

0.5000E 01 

0.0 

0.0 

0.0 

283.1 

0.5000E 01 

0. 3546E 02 

0.421SE 02 

0.4217E 02 

288.1 

0.5000E 01 

0.8865E 02 

0.1065E 03 

0.1054 E 03 

323.1 

0.5000E 01 

0.1773E 03 

0.2108E 03 

0.2105E 03 

373.1 

0.5000E 01 

0.3546E 03 

0.J218E 03 

0.4201E 03 

423.1 

0.5000E 01 

0.5319E 03 

0.U328E 03 

0.G2S2E 03 

473.1 

0.5000E 01 

0.70S2E 03 

0.8437E 03 

0.K341E 03 

523.1 

0.5000E 01 

0.H8C5E 03 

0.1055E 01 

0. 1036E 04 

573.1 

0.5000E 01 

0.1064E 04 

0.1266E 04 

0.1232E 04 

601.1 

O.SOOOE 01 

0.1174E 04 

0 . 1 3 JOE 01 

0.1347K 04 

273.1 

O.IOOOE 03 

0.0 

0.0 

0.0 

283.1 

O.IOOOE 03 

0.70B2E 03 

0.8437E 03 

O.H434E 03 

288.1 

O.IOOOE 03 

0.1773E 04 

0.2108E 04 

0.2107K 04 

323.1 

O.IOOOE 03 

0.3546E 04 

0.4219E 04 

0.421 IE 04 

373.1 

O.IOOOE 03 

0.70S2E 04 

0.N437E 04 

0.K402E HI 

423.1 

O.IOOOE 03 

0.1064E 05 

0.126GE 05 

0.1256E 05 

473.1 

O.IOOOE 93 

0.1418E 05 

0.1687E 05 

0.1668E or. 

523.1 

O.IOOOE U 

0.1773E 0$ 

0.2108E 05 

0.2073E OS 

$73.1 

O.IOOOE 03 

0.2128E 00 

0.25311. OS 

0.2464 E 0$ 

604.1 

o.ioooe or 
— — . 

0.2348E 05 

0.2793E 05 

0.2603E 05 


The «!»»!«>•*• oC supercooling, x, were con cd by Eq (1«), (IS), and (20) for given 

values of R/l’ f T . 

a 





TABLE 0-S COMPUTED DISTANCES OF SUPERCOOLING IN THE Ni-S» SYSTEM 
FOR THE THREE MODELS (SIMM 1 of S> 


G 

R/D 

K -275.1 

K -253.1 

K -298.1 

K -323.1 

K -373.1 

K -423.1 

K -473.1 

K -023.1 

K -5 73.1 

K - SOt.l 

S.S 

0.1 

0.0 

0.0 

0.0 

0.0 

7.4519 

-0.0000 

-0.0000 

9.0000 

-0.0000 

-0.0000 

5.0 

0.5 

0.0 

0.0 

2.5003 

6.8629 

14.1726 

21.2763 

0.0000 

0.0000 

-0.0000 

-0.0000 

s.o 

1.0 

0.0 

0.7452 

3.4314 

7.0863 

14.1545 

21.2767 

28.3690 

35.4612 

42.5535 

40.9507 

5.0 

2.5 

0.0 

1.3725 

3.5450 

7. 0922 

14.1845 

21.2767 

28.3690 

35.4612 

42.5535 

40.9507 

5.0 

5.0 

0.0 

1.4173 

3.5461 

7.0922 

14.1845 

21.2767 

25.3690 

35.4612 

42.5535 

40.9507 

5.0 

10.0 

0.0 

1.4154 

3.5401 

7.0922 

14.1545 

21.2767 

29.3690 

35.4612 

42.5 S 3 S 

46.9507 

5.0 

25.0 

0.0 

1.4154 

3.5461 

7.0922 

14.1845 

21.2767 

29.3690 

35.4612 

-.2.5535 

46.9607 

5.0 

50.0 

0.0 

1.4154 

3. 5461 

7.0922 

14.1545 

21.2767 

2*. 3630 

35.4612 

42.5535 

46.9507 

5.0 

100.0 

0.0 

1.4154 

3.5461 

7.0922 

14.1545 

21.2767 

2*. 3600 

35.4612 

42.5535 

46.9607 

5.0 

250.0 

0.0 

1.4154 

3.5461 

7.0922 

14.1845 

21.2767 

25.3694 

35.4612 

42.5535 

46.9507 

10.0 

0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

1.2507 

7.4519 

- o.oooo 

<0.0000 

0.0000 

10.0 

0.5 

o.e 

0.0 

©.© 

2.5603 

6.5629 

10.5549 

14.1726 

17.7251 

21.2762 

23.4751 

10.0 

1.0 

0.0 

0.0 

1.2502 

3.4314 

7.0563 

10.6351 

14.1545 

17.7306 

21.2767 

23.4753 

10.0 

2.5 

0.0 

0.5121 

1.7508 

3.5456 

7.0922 

10.6354 

14.1845 

17.7306 

21.2767 

23.4753 

10.0 

5.0 

0.0 

0.6503 

1.7728 

3.5461 

7.0922 

10.6354 

14.1845 

17.7306 

21.2767 

23.4753 

10.0 

10.0 

0.0 

0.7056 

1.7731 

3.5461 

7.0922 

10. tt 3 K 4 

14.1845 

17.7306 

21 . 2767 

23.4753 

10.0 

25.0 

0.0 

0.7002 

1.7731 

3.5461 

7.0922 

10.63*4 

•14.1845 

17.7306 

21.2767 

23.4753 

10.0 

50.0 

0.0 

0.7002 

1.7731 

3.5461 

7.0922 

10.6354 

14.1545 

17.7306 

21.2767 

23.4753 

10.0 

100.0 

0.0 

0.7002 

1.7731 

3.5401 

7.0922 

10.6354 

14.1845 

17.7306 

21.2767 

23.4753 

10.0 

250.0 

0.0 

0.7002 

1.7731 

3.5401 

7.0922 

10.6354 

14.1845 

17.7306 

21.2767 

23.4753 

15.0 

0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

25.0 

0.5 

o.o 

0.0 

0.0 

0.0 

1.4904 

3.5251 

5. sail 

G . MG 29 

*.3*19 

9.3004 

25.0 

1.0 

0.0 

0.0 

0.0 

0.7452 

2.6330 

4.1910 

5.6539 

7.0563 

*.509© 

9.3893 

25.0 

2.5 

0.0 

0.0 

0.5121 

1.3726 

2.5345 

4.2552 

5.6738 

7.0922 

*.5107 

9. 3901 

25.0 

5.0 

0.0 

0.1400 
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5.6738 

7.0922 

8.5107 
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0.2535 
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1.4154 

2.5369 

4.2553 

5.6735 
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0.0 
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0.0 
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0.0 
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TABLE 0-5 COMPUTED DISTANCES OF SUPERCOOLING IN THE Ni-Sn SYSTEM 
FOR THE THREE MODELS (Sheet l of 3) 


0 

R/D 

K-tTS.t 

K-283.1 

K-293.1 

K-3S3.1 

K -3 73.1 

K-423.1 

K=473.1 

K =523.1 

K-573.1 

K -604.1 

S.0 

0.1 

0.0 

0.0 

0.0 

0.0 

11.5665 

22.6949 

32.4305 

41.5216 

50.2909 

55.6391 

5.0 

0.5 

0.0 

0.0 

3.47TO 

8.3043 

16.8704 

25.3110 

33. 7481 

42.1852 

50.6122 

55.8532 

5.0 

1.0 

0.0 

1.1507 

4.1522 

8.4352 

16.8741 

25.3111 

33. 7481 

42.1852 

50.6222 

55.8532 

5.0 

2.5 

0.0 

1.6600 

4.2184 

8.4370 

16.8741 

25.3111 

33.7481 

42.1852 

50.6222 

55.8532 

5.0 

5.0 

0.0 

1.6070 

4.2185 

8.4370 

16.8741 

25.3111 

33.7481 

42.1852 

50.6222 

55.8532 

5.0 

10. 0 

0.0 

1.6874 

4.2185 

8.4370 

16.8741 

25.3111 

33.7481 

42.1852 

50.6222 

55.8532 

5.0 

25.0 

0.0 

1.6874 

4.2185 

8.4370 

16.8741 

25.3111 

33.7481 

42.1852 

50.6222 

55.8532 

5.0 

50.0 

0.0 

1.6874 

4.2185 

8.4370 

16.8741 

25.3111 

33.7481 

42.1852 

50.6222 

55.8532 

S.0 

100.0 

0.0 

1.6874 

4.2185 

8.4370 

16.8741 

25.3111 

33. 7481 

42.1852 

50.6232 

55. 8532 

5.0 

250.0 

0.0 

1.6874 

4.2185 

8.4370 

16.8741 

25.3111 

33. 7481 

42.1852 

50.6222 

55.8532 

10.0 

0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

4.9109 

11.5665 

17.3847 

22. 6948 

25.8133 

10.0 

0.3 

0.0 

0.0 

0.2269 

3.4770 

8.3043 

12.632? 

16.8704 

21.0920 

25.3110 

27.9266 

10.0 

1.0 

0.0 

0.0 

1.7385 

4.1522 

8.4352 

12.6555 
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21.0926 

25.3111 

27.9266 

10.0 

2.5 

0.0 

0.6954 

2.0981 

4.2184 
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12.6556 
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21.0926 
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10.0 
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2.1092 
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10.0 
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0.4937 

2. 2596 
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0.0 

0.0 

0.0 

2.3133 

4.5390 
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8.3043 
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11.1278 

25.0 
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0.0 
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0.0 

1.1567 

3.2430 

5.0291 

6.7417 

8.4352 
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11.1705 
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2.5 

0.0 

0.0 

0. 0054 

1.6609 

3. 3741 

5.0622 

6. 7496 

S.437© 

10.1244 

11.1706 

25.0 

5.0 

0.0 

0.2314 
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1.6870 

3.3748 
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0.8437 

1.6874 
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11.1706 

25.0 

50.0 

0.0 
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0.8437 

1.6874 

3.3748 

5.0622 

6.7496 

s.4370 

10.1244 

11.1706 

25.0 

100.0 

0.0 
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0.8437 

1.6874 

3.3748 

5.0622 

0.7490 
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10.1244 
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0. 0 
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3. 1770 

4.5390 

5.1627 

50.0 

1.0 

0.0 

0.0 

0.0 

0.0 

1.1567 
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3.2430 
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5. 0291 

5. 5639 

50.0 

2.5 

0.0 

0.0 

0.0454 

0.6954 

1.6609 

2.5265 
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3. 0622 
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50. 0 

5.0 

0.0 

0.0 

0.3477 
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3. 0622 

5.5853 
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10.0 

0.0 
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0.4152 

0.8435 

1.6874 

2.5311 
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4.2185 

5. 0622 

5. 5*53 

50.0 

25.0 

0.0 

0.1661 

0.4218 

0.8437 

1.6874 

2.5311 

3.3748 

4.2185 

5. 0622 

5.5853 

50.0 

50. 0 

0.0 

0.1687 

0.4219 

0.8437 

1.6874 

2.5311 

3.3748 

4.2185 

5. 0622 

5.5*53 

50.0 

100.0 

0.0 

0.1687 

0.4219 

0.8437 

1.6874 

2.5311 

3. 3748 

4.21 *5 

5.0622 

5.5853 

50.0 
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0.4219 

0.8437 

1.6874 

2.5311 

3.3748 
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5. 0622 
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0.0 

0.0 

0.0 
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0.0 
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0.0 

0.0 
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0.0 
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1.1567 

1.7385 
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2.5813 
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0.0 

0.0 

0.0454 

0.6954 
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1.6609 

2.0981 
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TABLE G-5 COMPUTED DISTANCES OF SUL HCOOUNG IN THE Ni-Sn SYSTEM 
FOR THE THREE MODELS (Sheet 3 of 3> 



3. KIKM 
- 0. 0000 
10 . 7 * 3 H 
10. 7*40 
10. 7*40 
10. 7*40 
10. 7*40 
10. 7*40 
10. 7*40 
10. 7*40 


- 0.0000 
- 0. 0000 
13. 12*9 

33.4240 
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13.4249 
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- 0 . 0000 
15. 4501 
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The deviations of critical temperature gradients for the linear model, from the 
quadratic model for R/D « 100, are listed in Table G-6. 

TABLE G-6 DEVIATIONS OF THE CRITICAL TEMPERATURE 
GRADIENTS IN THE Ni-Sn SYSTEM 


T, K 

GC(1)/GC(3)-1 

GC(2)/GC(3)-1 

273.1 

0.0 

0.0 

283.1 

-0.1591E 00 

0. 3376E-03 

298.1 

-0.1587E 00 

0. 8688E-03 

323.1 

-0.1578E 00 

0.1842E-02 

373.1 

-0.1559E 00 

0.4183E-02 

423.1 

-0.1533E 00 

0. 7259E-02 

473.1 

-0.1497E 00 

0. 1148E-01 

523.1 

-0.1446E 00 

0.1762E-01 

573.1 

-0.1364E 00 

0. 2738E-01 

604.1 

-0.1283E 00 

0. 3699E-01 


The deviations of the distances of supercooling for the linear models from the 
quadratic model are given in Tables G-7 and G-8. 
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TABLE G-T DEVIATIONS OF THE DISTANCES OF THE SUPERCOOLING FOR MODEL 1 

FROM THOSE OF MODEL 3 (X -x,)/x, IN THE Nl-Sn SYSTEM 
13 3 


G 

R/D 

10.0 

0.05 

10.0 

0.10 

10.0 

0.50 

10.0 

1.00 

10.0 

2.50 

10.0 

5.00 

10.0 

10.00 

10.0 

25.00 

10.0 

50.00 

10.0 

100.00 

25.0 

0.05 

25.0 

0.10 

25.0 

0.50 

25.0 

1.00 

25.0 

2.50 

25. 0 

5.00 

25.0 

10.00 

25.0 

25.00 

25.0 

50.00 

25.0 

100.00 

50.0 

0.05 

50.0 

0.10 

50.0 

0.50 

50.0 

1.00 

50.0 

2.50 

50.0 

5.00 

50.0 

10.00 

50.0 

25.00 

50.0 

50.00 

50.0 

100.00 

100.0 

0.05 

100.0 

0.10 

100.0 

0.50 

100.0 

1.00 

100.0 

2.50 

100.0 

5.00 

100.0 

10.00 

100.0 

25.00 

100.0 

50.00 

100.0 

100.00 

250.0 

0.05 

250.0 

0.10 

251.0 

0.50 

250.0 

1.00 

250.0 

2. 50 

250.0 

5.00 

250.0 

10.00 

250.0 

25.00 

250.0 

50.00 

250.0 

100.00 


K-283.l| K-298.1 K=323.l| K =373.1 K =423.1 K=473. 1 K=523. 1 K-573.ll K =604.1 





















































































DISCUSSION 


The segregation constant k. The solidus and the liquidus for the Ni-Sn system 
display moderate curvatures (Hansen 1958). The errors involved for linear model 
assumption are illustrated in Fig. G-2 through G-6. 

Figures G-2 through G-6 show that, for the pseudo-linear model, the errors of 

k, G and x are always approximately proportional to T , which in turn is proportional 
c o 

to the concentration of solute and zero for pure solvent. For the segregation coefficient, 

k, the error is about 0.3%/degree K, with the rate of change slightly decreasing at 

higher T . For G , the error is about 0. 08%/degree K, with the rate of change slightly 
o c 

increasing at higher T q . For x, the error is about 0.11%/degree K, with the rate also 

slightly increasing at higher T q . For the linear model, the error in k varies from 

+85% to -55%, the rate of change being about 0.4%/degree K. The error in G is always 

o 

negative at about -16% to -13%, the average rate of increase being 0. 01%/degree. The 

error of x varies from -25% to +38% for T = 0 to T * 331K, the average rate of change 

o o 

being about +0. 19%/degree. 


As expected. Table G-2 shows that k^ and k^ are independent of the melt tem- 
perature, while k increases by over four times (from 0.156 to 0.6192) as the modified 

u 

temperature T increases from 0 K to 331K. Table G-3 shows that, compared to k_, 
o o 

kj is overestimated by 86% at small temperature difference and underestimated by 53% 
at T = 331° (eutectic point); and k is always underestimated, the amount being pro- 

0 Vs 

portional to the modified temperature, T q . A relationship between — and 

k - k 3 

2 3 


k„ 


with T q is shown in Fig. G-2. 


The critical temperature gradient. From Eq (11), (12), and (13), we see that 

G , G , and G in Table G-4 should be linearly proportional to R/D and T , and 
C 1 C 2 C S 0 

G and G should also be proportional of (1-k). Table G-6 and Fig. G-3 show how 
C 1 C 2 


/ 


G-16 














Fig. G-4 Fractional Deviation in x vs To, for Model 1 

G-19 









Fig. G-6 Distance of Supercooling vs To, for Various G and R/D, 
Model 1 


G-21 
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much G and G >un ocviate Irom G • .',ioci lie ally, G ’s arc always underestimated 
C 1 1 2 c 3 C 1 

by about 13 ; but G 's arc always slightly overestimated from zero up to 3', . 


the distance of supercooling. When K I) x '3 in Eq »1HI, (ID), and ( 20 ), then 

-II / I) 

e x >■ 1 1 ■ Equations (18), (19), and (20) can then be solved for x. 

T 


Linear model: 


• Quadratic model: 


*1 G ,1 - k l» 


0 

x 2 : ir (1 

V 
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0 
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tM 
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< ; 
1 z 

1 Ai2 

n !2 

B 2 S 1 

X 3 3 ~G~ 

L 

2 A 22 T „ 

1 o A 22 

q 2 J 


The values of x’s in Eq (21), (22), and (23) correspond to the BC section in Fig. ( -<>. 

In this section, x is proportional to T but inversely proportional to G. The portion of 
Fig. G -6 for R D x< 3 corresponds to the AB section. If R D is so small that G ^G^, 
which corresponds to curve (1) in Fig. G-l, then there is no supercooling. This condition 
of no supercooling continues as R/D increases, until G - G^_. If R D increases again, 
then x continues to increase, which corresponds to curves (3) and (4) in Fig. G-l. 

Table G-7 and Fig. G-4 show that the fractional deviation of x , defined as ( x 1 " x 3 > * 3 , 
is underestimated bv about 3o ( at small values 01 modified me’t temperatures T, but 
overestimated by about 40^ at large values of T. Table G -8 and Fig. G-3 show that the 
fractional deviation of x , defined as (x -x )/x , is overestimated from 0 to 40G as 

2 2 0 o 

the modified temperature increases from 0 to 331K, i.e., as the melt temperature 
decreases from the melting point of pure Ni (1"26K or 1455 C) to the Ni-Sn eutectic 
temperature (1397K or 1124°C). 
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